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The Bending Stability of Thin-Walled 
Unstiffened Circular Cylinders 
Including the Effects of 
Internal Pressure’ 


HERBERT S. SuER,* LEONARD A. Harris,* WILLIAM T. SKENE,** AND ROLAND J. BENJAMIN*** 


North American Aviation, Ine. 


SUMMARY 


In a recent paper, the authors presented a statistical, semi 
empirical design procedure for the determination of the buckling 
strength of unpressurized and pressurized cylinders under axial 
compression. This procedure has been extended in the present 
paper to the bending of unpressurized and pressurized cylindrical 
shells and allows the calculation of the critical bending stress with 
a knowledge of the cylinder geometry and the internal pressure 
only 

Because no published data could be found, an extensive series 
of bending tests of pressurized cylinders has been performed 
These new data for pressurized cylinders are treated semiempir 
ically together with all of the other known test data for unpres 
surized cylinders. Best-fit curves are presented using applicable 
theoretical parameters Design curves for determining the 
critical buckli:g stress for unpressurized and pressurized cylinders 
in bending are then developed as 90 per cent probability curves 
from the test data 


SYMBOLS 


bending buckling stress coefficient in the equation 
ocr = Cp)E(t/r) for unpressurized cylinders 
= bending buckling stress coefficient in the equation 
Ocrp = ChpE(t/r) for pressurized cylinders 
increase in the buckling stress coefficient, Cy, caused 
by internal pressure (AChp) = Crp — CO 
= modulus of elasticity 
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moment of inertia of cvlinder 
length of cylinder 

unpressurized critical buckling moment 
pressurized critical buckling moment 

internal pressure 

radius of cylinder 

wall thickness of cylinder 

unpressurized critical buckling stress (¢ © M.,r/I 
pressurized critical buckling stress (a¢;p) = Merp r/1 
increase in the critical buckling stress caused by in 


ternal pressure (Ager = Ger; a 
INTRODUCTION 


Statement of Problem 


_ the last few years. the increased use of pure 
monocoque construction in aircraft and missile 
structures has forced a closer inspection of the stability 
analysis of unstiffened circular cylinders. It is well 
known that in many cases unsatisfactory agreement 
exists between the values of the critical buckling stress 
determined experimentally and those predicted by the 
classical theory. Experimental values of the critical 
stress are almost always lower than those calculated 
theoretically. For some loading conditions, large- 
deflection theories which include an initial imperfection 
parameter have been developed'~* to explain the dis- 
crepancy between test and linear theory. However, 
use of these more refined theoretical treatments is not 
satisfactory from a design standpoint because an estt 
mate must first be made of the probable initial imper- 
fection. To circumvent this design difficulty for the 
case of axial compression, the authors® have recently 
presented statistical, semiempirical design curves for 
the determination of the critical compressive stress of 








AERONAUTICAL 











282 JOURNAL OF THE 
PRESSURE /-— LOAD STRUT 
GAGE fL ‘ 
S| __ LOAD STRUT 
& V4 / PRESSURE INLET 
& ULE IP pressure InveT 
CLAMP gui | [HEAD J} 
____ GRCULAR 
SHELL 
qb 
COMPRESSION | oe’) 
LOAD CELL J 
©) | 
UI “HEAD ~===—CLAMP 





C 





CABLE ~ 
. : 

. ex (- --PULLeY 

LOAD STRUT——2 ‘%—BENDING LOAD CELL 


Fic. | Schematic drawing of test jig 


unpressurized and pressurized circular cylinders. 

In general, the behavior of an unstiffened cylinder 
under bending is similar to that of a cylinder under 
compression—that is, the buckling mode shapes are 
similar, snap buckling may occur, and the buckling 
stress is sensitive to initial imperfections. In addition, 
unpressurized cylinders under bending and under com- 
pression share the distinction of exhibiting equally poor 
agreement between small-deflection theory and test. 
More precise theoretical treatments which include the 
effects of large deflections and initial irregularities have 
not been applied to cylinders in bending, as far as is 
known, although it is probable that better correlation 
could be achieved between test and theory by their use. 
For pressurized cylinders in bending, no theoretical 
analyses have been published. 

Because of the difficulty of achieving a theoretical 
treatment of the problem, and because it was desired 
to compare the critical stresses under compression and 
bending, it is proposed in the present paper to treat the 
behavior of cylinders under bending semiempirically, 
in a manner similar to that previously used for cylinders 
under compression. Use of the statistical method of 
the analysis of cylinder stability data may be con- 
sidered to amount to an empirical evaluation of the 
effects of initial irregularities. The basic assumption 
in this evaluation is that the initial imperfections of the 
test cylinders were representative of those likely to be 
encountered in actual structures. Where 


sufficient data are available, such a procedure should 


aircraft 
provide reliable, statistically defined design curves. 


Scope of Investigation 

Because test data in bending were not available for 
short cylinders, the present paper has been restricted 
to long cylinders unaffected by end fixity conditions. 
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As a consequence, the theoretical parameters f1 
classical form of the buckling equation for long eyjj 
ders, ¢., = C,E(t r), have been used to develop a de 
sign curve for unpressurized cylinders. To aluat 
the test data, the experimental values of the 
buckling coefficient, C,, have been plotted aS a TUNCtio; 
of the r ¢ ratio. A best-fit curve was placed 
the test data and a 99 per cent probability curve was 


statistically defined. A thorough literature  seare} 
indicated that considerable test data were utlab] 
for unpressurized cylinders under bending and the r 


sults of these tests, with all but a few exceptions 
used in the study. However, as will be seen below, th 
data were scanty for values of r ¢ greater than 1,500 

In the present treatment of pressurized evlinders j 
bending, the parameters from the theoretical develo 
ment by Lo, Crate, and Schwartz‘ for pressurized cy] 
inders under compression have been used. The use 
these parameters would appear justified because of th 
similarity between the behavior of circular cylinders 
under bending and under compression. 

For pressurized cylinders under bending, apparenth 
no previous tests had been performed. A total of 5s 
tests were performed in the Structural Test Laboraton 
and these tests provide the only data known to be avail 


able for the bending stability of pressurized cylinders 





ABLE | 
Results of North American Aviati res 
1S-S Half-Hard Stainless Steel, E = 27 X I psi, ? 8.7 
21.54 

Spec b, Ser, Spec 

No." rit psi psi No ps sit 
Unpressurized 
210 1,006 0) 9,176 72(2 2,734 ( 1 218 
211 9,797 wi 62 
212 7,442 «(2 62 
213 1,006 7,838 78 2 
‘2 2,734 0 1 ,299 78(2 2,734 
Pressurized—Anxially Pretensioned 

5 1,006 1.5 11,290 89(2 2,734 } 
282 1.8 10,987 S6 8 16,9 
285 18.4 21,245 ] 8.5 23,5 
286 1,006 18.4 18,967 1] 14 l 
283 2,734 a.2 8,481 9 18 f 
2s4 2.2 9,019 lt 18 .( t 
85 t 4 9,742 13 IS 8,7: 
87 2,734 $.4 10,070 3 2,734 18 
Pressurized— Axial Pretensioning Balanced by Compress I 
975 1.006 1.0 9,70; S84 9 734 5,02 
276 1.0 11,649 73 8.4 li 
2770 i .&@ %2,437 71X(2 781 
278 1.8 10,842 76 8 
280 3.5 11,360 76(2 5 
281 3.5 12,425 77(3 7,112 
287 S 4 13,SOS8 ai(A4 6,98 
288 8.4 14,415 78(3 6,60 
208 16.4 14,82¢ 78 S 

9 1,006 16.4. 15,200 116 2 7,3 
272 2,754 ‘a 1, S41 17 2 7 oe 
274 b.3 3.475 4 Ss 
270 2.2 +, 965 79(2 
273 2.2 +344 Si! 7, (94 
274(2 3.5 1,758 SO (2 lt 8,11 
83 2,734 + 4 5,683 8] 2,734 2 8,1 

* Numbers in parentheses refer to the number of t suc 
ceeding test for those specimens tested more than on Eat 
succeeding test was made at an internal pressure equal t 


greater than any previous test of the same cylinder 


To = i? 


NDING STABILITY OF UNSTI 


ng pe ! For sign purposes, it appeared advantageous to 
lop a de. | present data for pressurized cylinders in two parts 
ies, i that] the critical stress which could be attained 
bendin if the cvlinder were unpressurized and that part by 
functi * | which the unpressurized critical stress is increased by 
through internal pressure. As a consequence, the pressurized 
UrVe wac data have been analyzed and a design curve developed 
Sear using the parameter (Ac., -) (ry t) as a function of the 
silabh pressure parameter (p £) (7 ¢)*, where Ao-., is the in 
| the a crease in the critical bending stress associated with the 
ns, wer internal pressure, p. This presentation has the ad 
low, th vantage of giving the designer a direct indication of the 
500 effect that internal pressure has in increasing the allow 
nders jr able critical stress. 
le vel [ 
ed cyl EXPERIMENTAL INVESTIGATIONS 
C use 
e of the North American Aviation Investigation 


vlinders Sa 
i Prior to the present investigation, there appears to 

have been no comprehensive testing of pressurized 

irently , : T , . . 
vlinders in bending. Therefore, an extensive series of 


l of BN F ° : . . » 
bending tests of model circular cylinders were performed 
rat TV or. 2 ° . 

“Y ~ at the Structural Test Laboratory of the Missile Devel 
opment Division to obtain test data within the param 


lers : saps ; 
eter range likely to be encountered in missile or aircraft 


A total of 5S tests was performed, 10 of unpressurized 
( 


vlinders an Iwo series 


| structures 
| 


{SN of pressurized cylinders. 
som, | f tests were conducted of pressurized evlinders in 
bending. In one series, the longitudinal component 
: i the internal pressure was balanced by an axial com 
pressive load before the bending moment was applied, 
18 while, in the second series, the longitudinal component 
bad of the internal pressure was not balanced. 
fi24 Che test evlinders were fabricated from 0.0032-in. 
” and 0.0OS7-in. thick, half-hard, 1S-S stainless steel foil 
The above sheet thicknesses were the averages of a num 
ber of micrometer measurements. The cylinders were 
3, 00 prepared with a length of 21.5 in. and a radius of 8.75 
in, and had r ¢ ratios of approximately 1,000 and 
1.6 730. The internal pressure was varied from 0 to 24 
psi. Most of the cylinders were fabricated by wrapping 
i hail the sheet and adhesive-bonding a single * ,-in. wide 
), 521 longitudinal seam with Epon VI adhesive. Some of 
the later cylinders were fabricated by seam-welding, 
008 which proved to be a faster and simpler technique. 
: No appreciable difference could be noted between the 


982 buckling stresses from the adhesive-bonded and from 
poe the seam-welded evlinders. The test cylinders were 
307 oriented with the vertical seam at the neutral axis. 

A schematic drawing of the test jig with a evlinder in 
The test evlin- 


' ' 


position for testing is shown in Fig. 1. 
119 ders snugly fitted the lower and upper heads of the test 

ng and were clamped to the heads by metal straps. 
uc- rhe straps prevented slippage and leakage between the 


icl h ] ° 
heads and the evlinder. 


Internal pressure was pneu 
matically supplied, but for reasons of safety most of the 
cylinders were partially filled with water for those tests 


1 


in which the internal pressure exceeded 2 psi. A 
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ABLE 2 
Summary of Previous Experiment Investigations 
No 
Spec ( 1 il Parameters 
Reference Material mens 
Donnel Brass and Steel 0. 90-28 8 5K 9 » 
al 17ST Aluminum 11 2.5-5 0 }-322 2.4-4.8 
17ST Aluminum 7 7 2 50 
Aluminum Alloy gs 0 IRs 2 0 
17ST Aluminum 5 rt ’ 80-749 0 25 





All specimens unpressurized 


manometer was used to measure the internal pressure 
In computing the pressure parameter from the test 
results, an internal pressure equivalent to a head of 
water of approximately half the specimen length was 
added to the applied pneumatic pressure reading for 
those specimens which were filled with water. Bending 
and axial compressive loads were applied by hydraulic 
struts and measured by SR-4 load cells. 

The geometrical parameters and the buckling stresses 
of the test cylinders are recorded in Table | Phe 
buckling loads of the test specimens were determined 
visually. For most specimens, a visual determination 
of the buckling load was relatively simple because the 
buckles often snapped into position. Occasionally 
buckles appeared to grow out of relatively large initial 
irregularities. Incipient buckling of a_ pressurized 
evlinder was often indicated by the formation of small 
circumferential ripples. Although initially very shal 
low, these ripples (similar in appearance to the shape 
of an axisymmetrical compression buckle) could easily 
be seen in the reflecting surface of the cylinder 

Photographs of typical buckle patterns in unpres- 
surized and pressurized cylinders are shown in Figs. 
2, 3, and 4. These diamond-shaped buckle patterns 
are typical for the buckling of unstiffened cylinders 
loaded in bending or in axial compression. Internal 
pressure changed the buckle pattern by decreasing 
the axial and circumferential wavelengths, although 
the longitudinal wave length was decreased more than 


the circumferential wave length. 


Summary of Other Investigations 


An important aspect of the present semiempirical 
treatment of unpressurized cylinders has been the ac 
cumulation of test data from past investigations. A 
large number of investigations were carefully reviewed 
to obtain test data, but only the experimental results 
summarized in Table 2 
Other ex 


from the investigations,* 
have been used in developing the analysis. 
perimental data were not included in the study for one 
or more of the following reasons 1) insufficient data 
were given; (2) test specimens were fabricated from non 
metallic materials such as paper, celluloid, or rubber 
which frequently exhibit variations in properties and 
thickness; (3) critical stresses were above the vield 
point; (4) critical stresses were above the proportional 


limit, and insufficient data were given to make a plas- 
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Fic, 2, Typical buckle pattern for unpressurized cylinder. 





Typical buckle pattern for cylinder with low internal 
pressure, 

















Fic. 4. Typical buckle pattern for cylinder with high internal 
pressure 


8 





NX a | 
3* THEORETICAL COMPRE E 
ee KLIN EFFICIE? 
' 
° DONNELL 
oT T >, wT ee — IMPERIAL 
° ©LUNDQUIST 
Py | | | * MOSSMAN & ROBERTSON 
2a rf @NAA 
5 — ——— ¢ PETERSON 
3 oe 
tae. || | 
a | 
ob °9 } 
4/42 $,———__+——__+_-_+ + +— 
Qe ' “a “a 
= = 44 of € 
E 4, } co > - 
a o ° 
= geo ; a a | ae | | _ 1 — 
° 
4+ go 
] ° a 
dq > 
= - 8 40 | oe + + 4 - 
o a 
° o,p 
av a 
| o & io} e 
q § 
me 4 + + 
} 
we 50C OO 1500 By Tay, 35, 40 
' 


Unpressurized bending buckling stress coefficients as a 
function of r/f 
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BENPING STABILITY OF % 


ticity correction; (5) the r f ratio was less than 100; or 


6) the L? rt ratio was less than 40. 


ANALYSIS OF UNPRESSURIZED CYLINDER DATA 


The buckling stress for a long unpressurized cylinder 


in bending can be expressed as 


where, for a particular case, the theoretical value of 
C, has been found by Fliigge'* to be about 30 per cent 


129 expert- 


higher than the corresponding classical compressive 
I In Fig. 5, 


wickling coefficient of 0.605. 
mental values of the buckling coefficient from the tests 
of six investigations are plotted as a function of the ry f¢ 
ratio. As is the case with compressive buckling, there 
is considerable scatter in the bending test data and the 
experimental values of the buckling coefficient are con- 
siderably lower than the theoretically predicted values. 
This large scatter and the discrepancy between theory 
and test has often been attributed to the effects of initial 
imperfections or defects in the test specimens. In 
general, experimental values of the buckling coefficient 
decrease with an increase in the y ¢ ratio, indicating an 
apparent increase in the magnitude and or the effect 
of the initial imperfections with increasing yr ¢. The 
relationship between C, and r ¢ will be established by 
empirically fitting a curve to the test data and deriving 
design curves on the basis of probability theory. 

In Fig. 6, the logarithm of the buckling coefficient, 
C;, has been plotted as a function of yf. In addition 
to more clearly defining the relationship between C, 
andr ¢, the semilog plot has the advantage of producing 
amore uniform scatter band in the range in which data 
A “best-fit” curve is shown on the graph 
The shape of this curve 


are available. 
of log C, as a function of r ¢. 
was selected by the investigators as being representative 
of the trend of the experimental data. The vertical 
location of the curve with respect to the test data was 
achieved by satisfying the equation 


\ 
> (y: — ve) = 0 (2) 
i=] 

where y; = log C,, and vy, = log C, for points from the 


best-fit curve at 7 ¢ of test point 7. If Eq. (2) is 
satisfied, it can be shown that the following relationship 


is also satisfied : 
> &; — ¥.)? = a minimum (3) 


Thus, by satisfying Eq. (2). the assumed curve was in 
the least squares position. Note that the best-fit curve 
may not be reliable for values of the 7 ¢ ratio greater 
than 1,500 because of the scarcity of test data. 

A statistical analysis of the data was performed after 
a best-fit curve had been established with the sum of 
IBM 


‘01 digital computer was used to solve the equations of 


the squares of the deviations a minimum. An 


the statistical analysis in the following form: 


NSTIFFENED 
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e PROB, CURVE FOR ‘ 
OMPRESSIVE BUCKLING (REF 6 \ 
90% PROBABLITY CURVE 
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500 1000 500 x SOO 3000 5 4000 
FIG. ¢ Unpressurized bending buckling stress coeflicients as 
function of + semilogarithmic grid 
y, v Kos }) 
where 
y, = log C, for a specified probability 
\ 
S = \ > (vi - ¥ V-1 
i=] 
° ») 
standard deviation 
N= number of test points 
K, = small sample probability coefficient 


The solution of Eq. (5) for the test points gave a value 
The 


1.47 was selected from a small sample 


of 0.130 for the standard deviation of the log C,. 
value of A, = 
statistical table’ 
The 90 per cent probability design curve is 


for a 90 per cent probability design 
curve. 
This curve agrees well with and gen- 
The 


99 per cent probability curve is also shown on this 


drawn in Fig. 6. 
erally parallels the lower bound of the test data. 
graph. The relatively large difference between the 90 
per cent and the 99 per cent probability curves is typical 
of a statistical analysis of test data with a large scatter. 

In Fig. 6, the 90 per cent probability bending design 
curve is compared with the 90 per cent probability de- 
sign curve for the compressive buckling of thin-walled 
evlinders from reference 6. It may be observed that 
there is not a uniform ratio between the buckling co 
efficients for the two loading conditions over the com- 
plete range of ry ¢f. At low values of r ¢, the bending 
critical stress is about 60 per cent greater than that 
for compression. This difference decreases with in 
creasing r ¢ to a value of about 25 per cent for large r f. 

In the procedure just described, it has been necessary 
for the investigators to select a best-fit curve which 
indicates the general trend of the test data. Because 
the data is scanty at large values of 7 ¢, the reliability 
of the selected shape is necessarily somewhat question 
able in this region. The similarity between the be 
havior of a cylinder under axial compression and under 
bending suggests that the variation of the buckling 
coefficient with the , ¢ ratio might be of the same shape 
for both cases. In compression, considerably more 
data were available, and it is of interest, therefore, to 


perform a statistical analysis of the bending data using 
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function of ¢/t (semilogarithmic grid 


the previously selected shape of curve® to represent the 
variation of the compressive buckling coefficient with 
Fit. 

In Fig. 7, the curve with the compression shape is 
shown in the least squares position for the bending data. 
A statistical analysis has been performed as before, and 
the 90 per cent probability curve is drawn in Fig. 7. A 
comparison of the 90 per cent probability design curve 
for bending based on the compression curve shape and 
the 90 per cent probability design curve for compression 
from reference 6 indicates that the bending coefficients 
can be obtained by increasing the axial compression 
coeflicients by about 35 per cent. However, the basic 
shape of the curve from the compression data does not 
fit the bending data well for all values of r ¢. There- 
fore, the constant ratio between the buckling coefficients 
determined in the above manner may not be reliable. 

Because the final design curve of this study is de 
fined in terms of a statistical probability of failure, 
it is incumbent on the user of this data to decide what 
probability is required by the design criteria demands 
of his particular structure. If the 99 per cent prob- 
ability design curve is used, there is a slim | per cent 
If the 90 per 


cent probability curve is selected, there is a 10 per cent 


probability of failure at ultimate load. 
probability of failure—associated with a_ significant 
weight-saving in the shell. 


ANALYSIS OF PRESSURIZED CYLINDER DATA 


The critical buckling stress of pressurized cylinders 
in bending may be expressed by a relationship analogous 
to that for unpressurized cylinders, 


Scrp = Cop (t/r) (6) 


/ 


where experimental results indicate that C,, is a function 
of the internal pressure as well as the geometry of the 
cylinder. An accurate theoretical treatment would be 
difficult to achieve and, as is the case with unpressur- 
ized cylinders, a prior knowledge of the initial imper- 
fections would probably be required for an accurate 
For this reason, 
the 


prediction of the buckling stress. 


an empirical treatment has been made _ using 
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parameter derived theoretically’ for the con 
The 
the behavior of cylinders in compression and in bendin 


| ITessive 


1 


buckling of pressurized cylinders. similarity jy 
suggests the validity of this approach. 

In Fig. 8, the values of the buckling coefficient ¢ 
are plotted as a function of the parameter (p /:) (r { 
The triangular symbols are used to indicate those tes 
2t by 


the cylinder heads. The 


specimens axially pretensioned to a stress of p) 
the 
circular symbols indicate those specimens to which ay 


internal pressure on 
axial compressive load, sufficient to balance the pre 
tension caused by the internal pressure, was applied 
prior to the application of the bending moment. |t 
should be noted that in computing the value of C,,,, th. 
critical bending stress o,,, has been obtained by diyid- 
ing the critical moment by the section modulus of the 
cylinder. For pretensioned cylinders, the actual mayi- 
mum compressive stress is found by subtracting pr 2 
from o 

A line is drawn in Fig. S to indicate the magnitude 
of the axial pretension stress pr 2¢in terms of C,,. The 
pretensioned specimens must always fail at a bending 
stress greater than o., = pr 2t, which represents the 
compressive bending stress required to balance the 
longitudinal tensile stress produced by the internal 
that, 
buckling coefficient for most of the axially stressed 
2t to the 


critical stress for the specimens with the axial pretensior 


pressure. It can be observed in general, the 


specimens can be obtained by adding pr 


balanced. At low values of the pressure parameter 
the axial stress has a negligible effect on the buckling 
stress coefficient. At large values of the pressure 
parameter, the effect of the axial component of the 
stress induced by the internal pressure becomes pre- 
dominant. 

In Fig. S a line is also drawn to indicate the general 
region of test data for the compressive buckling of pres 
surized cylinders. The experimental buckling coeffi 
cients for cylinders in bending are significantly greater 
than those for cylinders in axial compression. 

In order to better demonstrate the effect which in- 


ternal pressure has in increasing the buckling stress 
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coefficient, the experimental data are replotted in Fig. 0 a Sa cosa 
4 jn terms of the increment AC,, to the buckling co- gf en kk Ae 
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The Effect of an Internal Compressible 
Fluid Column on the Breathing Vibrations 
of a Thin Pressurized Cylindrical Shell 


J. G. BERRY* ann E. REISSNER** 
The Ramo-Wooldridge Corporation and Massachusetts Institule of Technology 


SUMMARY 
The free oscillations of a thin pressurized cylindrical shell con- 
taining a compressible fluid are studied here. The use of an 
approximate set of shell equations (shallow shell theory ) leads to 
a relatively simple formula for the natural frequencies of the 
coupled fluid-cylinder system. The results of some computa 


tions are presented 


(1) INTRODUCTION 


I THE FOLLOWING, we consider the determination of 
the inertia and compressibility effects of an internal 
liquid or gas column on the breathing vibrations of thin 
pressurized circular cylindrical shells. By breathing 
vibrations we mean those vibrations of the points of 
the cylinder walls which correspond to radial deforma- 
tions proportional to cos né or sin 0, where the integer 
n has at least the value 2. We do not consider rota- 
tionally symmetric vibrations corresponding to n = O 
or beam bending vibrations corresponding to m = 1. 

It has been shown in earlier reports’ * that it is pos- 
sible to analyze very accurately the breathing vibrations 
of thin pressurized cylindrical shells which do not 
contain an internal fluid by means of a system of ap- 
proximate shell equations. This system, which is con- 
siderably simpler than the classical equations as given 
by Love, consists of the equations of the theory of 
shallow shells, with an additional simplification due to 
the omission of certain of the inertia terms.* ' The 
nature of this theory is such that it is also applicable 
when the inertia and compressibility of an internal fluid 
are taken into account. It is this application which is 
carried out here. 

The present paper consists of three main sections. 
In Section (2) the equations for transverse vibrations 
of pressurized shallow cylindrical shells are discussed. 
The section concludes with a formula for the frequen- 
cies of free vibrations. This formula includes a term 
for an apparent inertia of the shell walls due to the 
presence of an internal fluid. 

In Section (4) the equations for the small motion of 
compressible circular cylindrical fluid 
The surface of the fluid column 


an inviscid 
column are discussed. 
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is assumed to be acted upon by a pressure which js q 


function of the axial and circumferential space coord; 
nates and time. In particular, we determine the fluc 
tuating pressure over the cylindrical surface of the fluid 
column in terms of the radial acceleration of points on 
the cylindrical surface. A specific expression is ob 
tained for the apparent inertia term appearing in the 
frequency equation of Section (2). 

Section (+) contains a discussion of the natural fre 
quencies of the elastic cylinder — internal fluid system 
A qualitative discussion is given concerning the effects 
of fluid compressibility and mass on the frequencies of 
the coupled system. The results of some numerical 
calculations are presented. 


(2) TRANSVERSE VIBRATIONS OF SHALLOW 
CYLINDRICAL SHELLS 


General Theory 
As discussed in reference 4, it is possible to reduce the 
system of differential equations which govern finite 
transverse vibrations of shallow shells* with middle 
surface equation 
2 = 2(x, y) (2.1 
to the following form: 


DV*w + m(O0*w/dt?) = pz + 
[(07z/Ox?) + (0?w/dOx?)] (0?7F/dOy?) 4 
[(0?s/Ov?) + (0?w/Ovy?)] (0? F/Ox?) — 
2[(07s/Oxdv) + (0?w/Oxdy)] (O?F/OxvOy) (2.2 
ViF = —Eh[(0*s/Ox?) (0?w/Oy?) + (07s/Ov") X 
(0?w/Ox?) — 2(072/Oxdy) (0?w/Oxdy) + 
(0?w/Ox) (0°w/Ov") — (O?w/OxNdY)*| (2.5 
In these equations, w is the deflection of points in the 
middle surface of the shell in the direction of the z-axis; 
F is Airy’s stress function which serves to describe 
middle surface stresses in the shell; m is the mass of the 
shell per unit surface area; p, is the intensity of surface 


loading on the shell in the direction of the z-axis; and 


D is the bending stiffness factor, given by D = Eh’ + 
12(1 — »?), where v is Poisson's ratio, E is Young's 


* A shell is said to be shallow when the ratio of rise to span ts 
less than 1, say 1/8. In the present case, it turns out that shal 
low shell theory gives good results? when the rise to span ratio Is 


as large as 1/2. 
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ind hf is the constant wall thickness of the 


modulus 
shell. The operator V+ is the iterated Laplace operator 
9!/dx* + 204/Ox*0y? + 04/Oy%. 


We assume that the load intensity p, is composed 
f two parts, p- po + pi, where fo is a constant in- 
ternal pressure, and /, is a fluctuating part which may 
be the interaction force between the oscillating shell 
and the oscillating internal fluid. The constant load, 
is resisted by membrane stress resultants .V,0, N,0, 
_ which depend on an Airy stress function, Fo, 


and N, 
as follows 
N 0?Fy/Oy?, Ny 0° Fy Ox?, 
Nryo = —(O?Fo/Oxdy) (2.4) 
From Eq. (2.2), the relationship between Fy and fo 1s 
py + (0°s/Ox*) (0? Fo/ Oy") + (O°%/Oy") X 
(02Fy),/Ox2) — 2(07s/Oxdy) (07? Fo/OxO0y) = O (2.5) 


The complete stress function F, corresponding to the 
loads po + fi, will be taken in the form 


F=Fh,+ Ff, (2.6) 


We now assume that the actual states of stress and 
displacement are small deviations from the state de- 
scribed by Fy and w = 0.* The substitution of Eq. 
2.6) into Eqs. (2.2) and (2.3) yields a set of equations 
which are then linearized with respect to the incre- 
mental functions F; and w. In this way, the system 
(2.2) and (2.3) reduces to 
m(O*w, Ot”) pi + (07?Fo/Oy?) X 

(O?w Ox?) + (0°Fo/ Ox") (0°, OV") — 


2(0°Fy, OxOy) (O7w/Oxdy) + (07s /Ox*) (0? F)/ Oy?) + 


DV‘w + 


(0°2/Oy?) (0? F,/Ox?) — 2(0%s/Oxdy) (0? F,/Oxdy) (2.7) 
VF, —Eh|(07s/Ox?) (0?w/Oy?) + (072/Oy") X 
(O?w/Ox?) — 2(07s/Oxdy) (O?w/AxdY)}| (2.8) 


In arriving at Eqs. (2.7) and (2.8), we have made ex- 
plicit use of Eq. (2.5). 
Specialization to the Circular Cylindrical Shell 
Eqs. (2.7) and (2.8) are now specialized to the case of 
a circular cylindrical shell by setting 
2/H = 
= H — (y?/2a) 


1 — (y?/b?) (2.9) 


or Z (2.10) 
where // is the rise of the shell, its half width, and a 
is the meridional radius of curvature. A comparison 
of Eqs. (2.9) and (2.10) shows that the rise //, width 


26, and radius a are related as 


6? = 2a (2.11) 

We now imagine that the complete cylindrical shell, 
of which the shallow portion of projected width 20 is a 
part, deforms according to the pattern cos @. Fur- 


thermore, we stipulate that the region —7/2 < n@ < 


a : 
w = 0 must be replaced by a more general expression unless 
VF, = 0 


In what follows, we assume that V#F) = 0 


x/2 corresponds to the region —b < y < band there- 
fore, in view of the shallowness of this portion, we have 


y ag 


and thus b ma/2n 


Eqs. (2.7) and (2.8), with z given by Eq. (2.10), now 
assume the following form: 
DV*w + m(0*w/Ot*) pi + (O0°Fo/Oy?) X 
(O°w/ Ox") + (0°Fo/Ox*) (O?w/Oy?) — 2 X 
(0°Fy/Oxdy) (07w/OxOy) — (1/a) (0°? F)/Ox" 


Ox") 


2.14) 


VF, (Eh a) (O*w 2.15) 


We finally stipulate that, for the class of problems under 
consideration here, 0°/)/OxOy 0; the load in 
tensity py does not give rise to shearing stresses .V 


1.€., 


In addition, anticipating a result from Section (3), we 
write for the fluctuating pressure component /), 
pi —m,(O07?w/ Ot*) 2.16) 


Eq. (2.16) states that the interaction pressure exerted 
by the internal fluid may be considered as an apparent 
mass contribution, m,, to the mass per unit area of the 
shell. In view of Eq. (2.16) and the fact that 0°F; 
Oxdy = O, Eqs. (2.14) and (2.15) become 


DV tw + (m + my) (0?w/OF?) N,0(O0*w/Ox*) + 


N,0(0?w/Oy?) — (1/a) (0?Fi\/Ox?) (2.17) 
VF, (:h/a) (0?w/Ox?) (2.18) 
We note that for a constant internal pressure fo, the 


quantities VV, o and N,9 are constants. In particular, 


for a shell with closed ends, 
Nix (2.19) 


Nw = pod, (1/2) poa 


The Frequency Equation 
We now assume that the motions of the shell wall < 

such that w and F, have the form 

(2. 


w = A sin (rx/L) cos (ny/a) sin wl 


F, = B sin (rxv/L) cos (ny/a) sin wt (39 


Substitution of Eqs. (2.20) and (2.21) into Eqs. (2.17) 
and (2.18) yields the following two simultaneous alge 
braic equations for the coefficients A and B: 
A\D[(x/L)? + (n/a)?]? + Nyo(2/L)? + Ny X 

(n/a)? — (m + my)w} — Bi(1/a) (x/L)?} 0 


A}\(1 a) (tr L)*} a 
B}(1/Eh) [(#/L)? + (n/a)?}?} 0 (2.23) 


99» 


The vanishing of the determinant of the system ( 
and (2.23) leads to the following frequency equation: 
Eh (r/L)* 


a? [(r/L)? + (n/a)?}? 


nN 
+ ( 
a 


(777 os Mw" = 


ol(z) 
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The first term on the right in Eq. (2.24) represents the 
effect of membrane action, the second term the effect 
of wall bending, and the third term the effect of the 
pressure /. 

Due to the simple form of Eq. (2.24), we can easily 
draw conclusions concerning the relative importance 
of the various terms. With regard to pressurized cyl- 
inders, the conclusion that is of most interest is as fol- 
lows. If V,o is of the same order as .V,o, then wall 
bending will have a negligible effect on the frequency, 
as compared to the effect of po, providing 


No[(3r/L)? + (n/a)?] > D[(r/L)? + (n/a)?] 


(2.25) 
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(2.26) is satisfied, then Eq. (2.24) reduces to 
Eh (3/L)# 


(m + My)wr - 
a? [(x/L)? + (h/a)?|? 


3) SMALL OSCILLATIONS OF A COMPRESSIBLE CIRCULAR 


CYLINDRICAL FLUID COLUMN 


General Theory 


Assuming small motions of an inviscid compressible 









































or, with No poa, D Eh*/12(1 — v?), fluid and linearized theory, we have as basic equations 
(E/ po) (h/a)® < 12(1 — v*)/[n? + (ra/L)*] (2.26) py(Ov,/ OF) —(Op Ov) 3.1) 
If we assume that the order-of-magnitude relation py (Ov, / OL) —(Op/Or) 3.2 
Table 1 Values of K waft Table 3 Values of K walt 
Po? ' ‘ Po* 0.200, C = 4.00, ~ 
os 0.200, C = 4.000, = 0 ri ; ; 500 
R 7/10 | n/4 | n/2 30/4 a 7 3a/2 2a AS w/i | n/4 w/2 | 3n/4 7 30/2 2 
| 10. 03750 0. 2026 0.6163 |0.8772 {1.018 1.284 iF 435 | j0. 1074 | 2471 0.5992 857¢ 1. 026 1. 306 i 
, [0.08047 | 0.5256 0.7234 1.077 1.231 1.478 1.78 a fen te b 8179 0. 906 1. O7¢ 1. 242 1.499 
| 0.7949 0. 8169 0.9035 {1.788 }1. 863 2. 050 | 2. 308 1. 687 696 1.731 1. 788 1.864 2. OF 
| |0.03932 |0.2182 0. 6230 an 1.162 1. 384 }1.48 1476 0. 263 0. 6437 0. 9668 1.18 1.4 
| 
0.01715 |0.1015 0. 3405 0.5938 [9.8250 1.128 r en 2171 0. 3989 645 8693 1.19 
| 1.077 1.092 1.146 1.231 }1. 342 11. 61( pen 1.078 }! 1.1 1.23 1. 344 1.61 
” 12.083 2.021 2.049 2. 096 12. 160 2. 333 2.554 2.01 2.021 2.049 2. 09¢ 2.16 2 
| |0.01771 |0.1047 0. 3514 0.6230 0.8541 162 1. 33 0.2134 | 0. 2398 0.4179 0. 668¢ 894 1.205 
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|. |1.879 1.888 1.918 1.968 | 2.037 |2. 224 | 2.457 PD eeete ssi 1.923 1.975 2. 044 2. 22 
2.941 2. 946 2.966 |, 998 3. 043 3.167 |3. 333 2. 941 shia 12. 966 2.998 43 ! 
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Table 2. Values of K = w aye Table 4. Values of K = wa ig 
Po@ P P,8 ; ; ‘ 
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| T 7 7 Ral a ] ; T ei a i 
fs | #/t0 | ee | 34 | 7 7 7 ESS Ba | 4 | #2 sa/4 | - | 3x2 20 
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po(Ov'g/ OF) —(1/r) (Op/08) (3.3) 
Co"(p — po) 3.4) 


Op/Ot) 7 po) (1 2) [O(rv,)/Or 1. 
1/7) (29/8) + (Qv,/dx)} = 0 (3.5) 


wher¢ 
ry. 6 cylindrical coordinates 

velocity components 
reference pressure 
reference density 
reference sound velocity, |(¢dp dp) Co? | 
time 

Eqs. (3.1) to (3.5) may be reduced to a single differ- 


ential equation for the pressure fluctuation g, defined as 


b— p 5.0) 


This equation is 


0°g/Or?) + (1/r) (Og/Or) + (1/r?) (0°g¢/08") + 
(0°q/ Ox") (1/co?) (0?g/Ot?) (3.7) 
We note that for an incompressible fluid, c while 
for typical compressible gases, 
Co A 3 X 104 in./sec. (3.8) 


Solutions of Eq. (3.7) which are suitable for our pur- 


poses are of the form 


g e™ cos né sin (rx/L)Q(r) 3.9) 


where LZ is the axial half-wavelength. Substitution of 
Eq. (3.9) into Eq. (3.7) yields an ordinary differential 
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Al,,[V (9?/L?) — (w?/co?)r 
O(r) 3.11) 


AJ,(V («?/co”) + (9?/L*)r 


where /, is the Bessel function of the first kind and /, 
is the modified Bessel function of the first kind. 


Expression for an Apparent Mass Factor 


For purposes of determining the interaction effects 
between the fluid and oscillating shell, we express the 
fluctuating pressure g at the radius r a of the fluid 
column in terms of the radial acceleration Ov,/Of at 


r a. In view of Eq. (3.6), Eq. (3.2) becomes 


Ov,/ Ol -(1/ po) (Og/ Or) 3.12) 
therefore, with 


g rit os J,{ cos nO sin (wx / Le 3.13) 


Eq. (3.12) becomes 


V (x?/L?) (w?/co")1 
—(A/ po) x 
V (w?/co?) — (9?/L*)J,,’ 


mx/L)e™ (3.14) 


Ov,/ Of 
cos né sin 


where a prime denotes differentiation with respect to 
the argument of J, or J,. Eqs. (3.13) and (3.14) may 
be combined to yield the following expression for g(a), 
the fluctuating pressure at r = a, in terms of @,(a), the 


radial acceleration at r = a: 














ain nl I,(V\?2 —k 
equation for V(r), <> 
) ) / , 
20) /dp? i 10/d Vd — RI,” (VA — RF 
de / —+- | { _ .) 
m/ar') 4 r) (dQ/dr) q(a) — ppd ,(a) 
F , , 2/9 ' 9/91 7 9 » 9 
[(a?/L?) — (w?/co?) + (n*/r7)|Q = 0 (3.10) Jn(V k? — x 
k>xX 
) 9 , ) 
‘ , —s Pi ; : . ‘ 2 ( ‘ 
Solutions of Eq. (3.10) which are regular at 7 O are ot Vk A? J Vi A 
- ) | 
the form 3.19) 
Table 5. Values of K = wa/~ Table 6. Values of K = wayf 
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where k = w/t (3.16) 
js a reduced frequency and A = ma/L. Eq. (5.15) 


mav be rewritten in the form 


g(a) = —praf,(r, k)w (3.17) 


where f,,(A, &) is defined by reference to Eq. (3.15) and 
w is the radial deflection at the surface of the fluid 
column—i.e., W = v,(a). 

The fluid column and the elastic shell are now coupled 
by imposing continuity requirements at the interface 
of the fluid and the surrounding shell. These require- 
ments are 
(3.18) 


Weshell = Wiuid = W 


Pi = q\a), 

In view of Eq. (2.16) we then have for the apparent 

mass contribution m, in Eq. (2.24), 

my = poaf,(X, k) (3.19) 

(4) FREQUENCIES OF THE COMBINED CYLINDER — FLUID 
COLUMN SYSTEM 


Some General Conclusions 


In order to facilitate a discussion of the effect that 
an internal fluid has on the frequencies of a cylindrical 
shell, we recast Eq. (2.24) as follows. Let 


Q2 = paa/Gh, K = wa V p./G = wa/C., 


E = C;/Co, k =< CK, (4.1) 
m ph, X= wa/L, Ny = pod, 
Nr = (1/2)poa, 1 = poa/psh 


+ y) is the shear modulus and p, 
Using the defini- 


where G = £/2(1 
is the density of the shell material. 
tions provided by Eqs. (3.19) and (4.1), we rewrite 


Eq. (2.24) in the form 


K*{1 ™ nt nA, CK) } = {5 3) [A4 (n? + ?)?] 4 


(1/4) (h/a)* (A? + n?)? + Q[n? + (1/2)A7] (4.2) 


where we have assumed that Poisson's ratio has the 
value 1/3. Wenow let 
K,? = (8/3) [A4/(n? + A2)2] + (1/4) (h/a)? X 

(A2 + m?)? + Q[n? + (1/2)A?] (4.3) 


or, if the order-of-magnitude relation (2.26) is satisfied, 


K,2 = (8/3) [A4/(m? + dA2)2] + Q[n? + (1/2)d2] (4.4) 
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K,, is the value of the frequency parameter A for a 
pressurized shell vibrating in vacuo. Finally, we write 
Eq. (4.2) in the form 


f,(\, CK) = (1/n) [((C?K,2/C?K?) — 1 1.5) 


or if, for the sake of brevity, we define 


g,(CK) = (1/n) [((C?K,?2/C?K2) — 1 (4.6) 


f(x, CK) = g,(CK) (4.7) 


then 


When both sides of Eq. (4.7) are plotted against CA, 
a system of curves similar to that shown in Fig. 1 is 
obtained. 

The solutions of Eq. (4.7) are those values of CA at 
which the curve g,(CA), see Fig. 1, intersects the 
various branches of the curve /,,(A, CA). 

Referring to Fig. 1, the following qualitative obser 
vations may be made. 

(1) The intersection of the g,-curve with the CK 
axis gives the value of the frequency parameter A in 
the absence of any inertia effect of the internal fluid. 

(2) The intersections of the f,-curve with the CK- 
axis give the values of the frequency parameter A for 
vibrations of the fluid column in the absence of the 
surrounding thin shell. The first intersection corre- 
sponds to a vibration with one nodal circle, the second 
The 


existence of these vibrations depends on finite com 


to a vibration with two nodal circles, and so on. 


pressibility of the fluid. 

(3) The values of CA at which f, becomes infinite 
give the values of the frequency parameter A for vibra- 
tions of the fluid column if the surrounding cylindrical 
shell is assumed to be rigid. The existence of these 
vibrations also depends on finite compressibility of the 
fluid. 

(4) The lowest frequency of the fluid vibration with 
rigid enclosing cylinder has a lower value than the cor- 
responding frequency for the fluid column with a free 
cylindrical boundary. 

(5) The effect of inertia and compressibility of an 
internal fluid is to /ower the frequency of the wall vibra- 
tions of the shell in comparison with the value of this 
frequency when inertia and compressibility of the in- 
ternal fluid are absent. In addition to this lowered 
frequency, there occur higher frequencies each of which 
is associated with a particular number of nodal circles 
in the oscillating fluid. 


An Approximate Frequency Formula for Heavy Liquids 


Inspection of Fig. 1 reveals that as the value of the 
mass ratio n increases, the lowest value of the frequency 
parameter decreases. When 7 is quite large, we may 
with good approximation replace the actual value of 
the function f, by its value for A = 0 and determine A 


from the relation 


K = K,,/Vnf,(d, 0) +1 (4.8) 


The second term inside the square root will be negli- 
gible, and we have in first approximation, when | < 7, 
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the following formula: 


K = K,/V nfl, (d)/AT, (A) ] (4.9) 
for that frequency of vibration which takes place when 
no nodal circles are present in the fluid. 

We also note that the higher frequencies, for vibra 
tions with nodal circles, are very nearly equal to the 
frequencies of vibration of the free-boundary column of 


liquid. 


Numerical Results for Some Special Cases 


The first three roots of Eq. (4.5) have been obtained 
for some specific values of 7 poa/psh, C = c;/co, and 
K,. The computations were carried out for the fol- 


lowing three pairs of values of n and C: 


n 0.200, C = 4.0; 
n> 0.464, C = 9.3; 
n = 300.0, C = 2.5 


and for values of K,,, as defined by Eq. (4.4), associated 
with the following values of V, A, and Q: 


i= 3 a. 4 ae: 
Q = poa/Gh = 0,1 ‘400, 1/200, 1/100; 
h = ra/L = w/10, w/4, 1/2, z, 32/2, 2 


The results of these computations are presented in 
Tables 1 through 12. The four rows for each value of 
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n in these tables have the following significance: the 


first three rows represent vibrations with zero, one. and 
two nodal circles in the fluid, respectively; the fourt, 
row contains the corresponding value of K,,, whic 


holds when inertia and compressibility of the interna] 
fluid column are neglected. A comparison of the ny 
merical values in the first and in the fourth rows shows 
that the effects which are determined in this report 
may be quite appreciable. . 
Tables 10, 11, and 12 contain only the values of K 
(4.5) and K 


The second and third roots of Eq. (4.5) are the same 


corresponding to the first root of Eq. 


to within four figures, for all values of 2 considered 
here—1i.e., the same as given in Table 9. 
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listing in the next edition of the Roster. 
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This is your last opportunity to be listed correctly in the 


that appeared for them in last year's Roster. 
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Random Gust 


s, Including a Sample Problem 


K. D. SAUNDERS* 
Douglas Aircraft Company, Inc. 


SUMMARY 


A general power-spectrum equation for stationary random 
gusts is derived on the basis of a mathematical abstraction that 
iodels the probability density distributions of observed mete 
ological data The 
n terms of the modes and moments of simple probability-density 
Sta 


50,000 


scale and intensity of turbulence are given 
and occurrence 
altitudes of 
and thunderstorm 


listributions for the magnitude, duration 


tistical parameters are tabulated for use to 


ft. Extreme spectra are given for clear-sky 
gusts 
A sample 


factors at the c.g. of a swept-wing fighter in flight through low 


problem is worked for expected incremental load 
level random gusts at Mach Number 0.86 


INTRODUCTION 


( * at power-spectrum equations are used in getting 
outputs—e.g., load factors, pitch 


attitudes, airframe stresses, autopilot outputs, etc. 


incremental 


of an airplane in flight through random gusts. The 
product of the gust spectrum and the modulus-squared 
of the transfer function is the output spectrum. The 


integral of the output spectrum over the frequency 
half-plane is representable as the variance of the out- 


put of a normal distribution with mean zero. The 
probability of an output exceeding a given value is ex- 


pressible in terms of the complimentary error function. 
The expected number of times per second an output 
passes through a given value is measured by the output 
for instance, references 1-9. 


spectrum—see, 


POWER SPECTRA 


Consider a model of rectangular pulses that satisfy 


probability-density functions for velocity, duration, 
and occurrence, p(v), p(/), and p(z). Assume the 
pulses are uncorrelated so that 
(v; — 9) (v; — 3) = 01+ F 
= 9 2 1 = 7 
Truncate the model in the interval 27, 7 > ¢ Let 


v, stand for pulse of duration / at time /;. Let 7 stand 


for a time shift overlapping v; with v,(¢ + 7), Fig. 1. 
Let R(r),; be the autocorrelation of this pulse. Then 
l eliti 
Riv). = — kn v,7dt 
2T r+ JStits 


= — lim (/ — |r/)o,7, |r] </ (2) 
27 


~ 7—-_> 
= U, |r| >i 
Received August 6, 1957 
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r takes on the absolute sign since R is even. The 


expected value of R(7),; is 


1 — 
R(r); = . lim p(l) (1 — |r) )plw)erdl dz 
27 | in vie 
Shifts of 7 greater than / add nothing to R(7r),, be 
cause of Eq. (1), so 
. . 
R(r),; = lim } | pil) X 
- /- ei rT « 
(1 — |r| )p(v)vrdl dv (3 
Now returning to the ensemble of pulses: the ex 
pected value of v,;° is the variance of v, ¢,*, p(n) gives 
the mean 7, and p(/) gives the mean /. The number 
of pulses like v; in an ensemble referenced to an airplane 
in flight at speed |’ for time 27° is 7127. So the auto 
correlation for the ensemble is 
R(r) = aVo,? | (1 — |r )p(l)dl (4) 


Let /,, be the mid-point of the least class interval of 
pl) the pil 
show that the probability-density 


2. These 


and £, be mode of Investigations 


functions for p(/) 
are shaped as shown in Fig. functions fit a 


(/ — |,) product exponential (/ — /,,) funetion with two 


undetermined coefficients. To get the function, take 
pt!) as 
pl) = 0, l<l, 
= a(/ — /,) exp, —d(1 — l,i, 1 >), (5) 
then impose constraints 
| p()jdl = |] 
a (6) 
dlp(/)| 
[f ) 0 
a Jiat 
From Eqs. (5) and (6) it follows that 
pl) = [1/(& 2] (7 —1,) & 
exp {—[(! — In) /(&r — Iw} (7) 


Substituting Eq. (7) in Eq. (4), performing the inte- 


gration, and simplifying gives the autocorrelation of 


v(t), 
R(r) = 


nVo,*; [|r 
. " yt [ . , 
[Im /(&1 — Im) ] + 2} exp [—|7| (Er — Ind] (8) 
Eq. (8) is the autocorrelation of the ensemble of pulses. 
Writing /,, = Vim, Vie. — ! 
transforming to obtain the power spectrum of v(f), gives 


and £;, — Im; = m), and 
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Now observe that the data-reduction procedures reduce a maximum frequency 71” /,,.. 
The sample trace in Fig. 3 demonstrates the filter operation in the reductio, 


crossing irregularities are not modeled. 
of the trace to an ensemble of rectangular gusts. 


The time trace is smoothed to a maximum frequency content 71” /,, 


ABRRONAUTICAL SCIENCES 


MAY, 1958 


Irregularities smaller tha; 


The maximum harmonic of the smoothed trace is shown dashed jy 


Accordingly, 


Vio) = W (xV/1,,) = 0 
A filter operation is required on Eq. (9). A natural filter to use is a one-minus-unit-step at w = rl” /, 
The final form of the power-spectrum equation is 
fa, 2(E i, — 'ms)? EXP tims (Ets — Ins) f., — |.) /, (E),, — Ins)? ]I 
W(w) = Ps - eee a Tr wo “- l + w P 4 
mV [lL + wrt — Ims)?/ V?2? | | E — Ie L V2 ( 


For turbulence with small eddies, /,,, — 0, and Eq. 


(10) becomes 


Bne 7E,2[1 + (we? /3 V7) | 


: 11 
rV[1 + (w%,,2/V2)]? sii 


W(w) 


In this form &,, is recognized as the scale of turbu- 
lence, usually denoted L. The spectrum at w = 0 


gives the intensity of turbulence 
W(0) = (3% o,” £152)/aV (12) 


The units for I1’(0) in Eq. (12) are (ft. see.*) (rad. ft.). 

In Eq. (11), where /,,, — 0 for turbulence with small 
eddies, it is readily seen that this equation is, except for 
a constant, the approximation for the Karman- 
Howarth" power spectrum for isentropic turbulence of 
fluid-velocity fluctuations normal to the direction of 


The spectrum is written 


_ o) [1 + (80°L?/V?)] 
7 Vl + (w*L? v2)} 


mean flow. 


o(w) (13) 
This equation is usually used to represent the power 
spectra of gusts, sometimes putting 


¢(0) = 780, L = 300 


Sometimes other values are used for ¢(0) and ZL. Re- 
gardless of the selection of these parameters, Eq. (13) 
does not lead to a prediction of outputs. The moments 
of the output spectra are divergent. 
Momentarily turning to p(m): the 
density distribution for m closes on n = 0. 
tuting » for / in Eq. (7) and putting /,, = 0, gives the 
distribution for p(m). Now let &, denote the mode of 
p(n). Then f is 2é,. With this substitution for # in 
Eq. (10), the power spectra of gusts are given in terms 
of the moments and modes of p(v), p(/), and p(m); 
When the variable of transformation is 


probability- 
Substi- 


Les, and V. 
w/ V’, then speed is out. 


SAMPLE PROBLEM 


An airplane in flight through random gusts for cer- 
tain wanted outputs is representable as a linear system. 
The gusts subject the airplane to random lifts and 





moments. If the wanted output is c.g. incremental 
load factor, then the airplane is adequately represented 
by the linear, rigid-body, longitudinal equations oj 
motion. The typical output spectrum of c.g. incre 
mental load factor in gusts is shown in Fig. 4. The 
spectrum is concentrated about the short-period fre. 
quency, showing that the motion is rigid body. 

In this section, a sample problem is worked demon 
strating the use of Eq. (10). Expected c.g. incremental 
load factors are predicted for a swept-wing fighter in 
flight through low-level gusts. The airplane flight 
condition is 14,000 Ibs. gross weight, 16,500 slugs-ft. 
pitch inertia, c.g. at 25 per cent M.A.C., Mach Num- 
ber 0.86, free in linear, rigid-body, longitudinal motion. 
Starting from the transformed equations of motion, 
the problem of predicting expected outputs leads nat- 
urally to Rice’s equation for the number of times per 
second an output equals or exceeds a given value. 


Eq. (10) is zero for frequencies > 71” /,,, rad., so the 


xX, = | T "W(w)dw 
0 
J, T|?W(w)w*de 
a 


converge and can be used to predict expected outputs. 


integrals 


Il 


Ix, 


9 
on 


Gust-spectrum equations a | w* on the interval (0, = 
do not give convergent integrals and are not useful for 
this purpose. 

Start the sample problem with the matrix equation, 


AX = B: (15 
where 

1 — (Vp/g) (1 /g)p 
A = 10 (V+ 2,)p (2, — 1)p + 2, 

0 (-p+M,)p Map + M, 
X = column (A,, 6, w) 
B = column (0, Z4) + Z,, Map + M,) 
4. = incremental load factor 
6 = pitch angle 
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POWER-SPECTRUM EQUATION FOR 
i’ plunge velocity 
transformed gust velocity, v(p) 
M 
M 
M ‘i a eee 
dimensional stability derivatives 


The transfer function, A,/V, is obtained from Eq. 
15). But first it is convenient to rewrite this equa- 
tion. Let 


repeated indices denote summation on the repeated 


Let Y, stand for any one of the unknowns. 


symbol. Now rewrite Eq. (15 
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But, A jjQir = |ald 
where (a) is determinant of A, and 6,, is the Kronecker 
delta, so that 
a|\X, = Aydt 
and T(X,/v) = A;id,/\a (16) 


The symbol 7 is placed before X, |” to denote that 
it isa transfer function for X, and v. 

Eqs. (10), (14), and (16) give information for deter- 
mining the expected number of times per second that 


an output equals or exceeds a given value, 


N(X, > Xn) = 
(1/2m) (1x,/X,2) 


exp) —X,,2, 2X7} (17) 
Eqs. (14) and (17) are given in reference 1. 
It remains to put values in Eq. (16) for 7(Y,/v) and 


obtain a spectrum equation from Eq. (10) for flight 


A ;i,X, = Ajjd;t through low-level gusts at 963 ft. sec. 
The expansion of Eq. (16) for 7(A, v) gives 
T(A./2) = —(1/g)p(p? — Map) (20 + 2%) — (p? g) (Mu + Map e4 (18) 
(p? — Mp) (Su — pee — p) + (Mab + Mi) (pz, + V, 
Putting 
V = 963 ft./sec., g = 32.2 ft./sec.*, 2 = —12.5 rad.—'|, 2, = —3.88 rad. 
Zs = —0.002 rad.-!, M, = —4.87 rad. M, = —0.014 rad.-!, MM, = —0.002 rad. 
and substituting in Eq. (18) gives the result, 
—0.00006p4 — 0.120p* — 0.580p° (19) 


a (An 
1b 
v p 


Eq. (19) is the gust velocity-c.g. incremental load- 
factor transfer function for a typical swept-wing fighter 
in flight through low-level gusts at 14,000 Ibs. gross 
weight, 16,500 slugs-ft.” pitch inertia, c.g. at 25 per 
cent M.A.C., at Mach Number 0.86, when free in linear, 
rigid-body, longitudinal motion. 

The next item required is a spectrum equation from 
Eq. (10). 

In most cases flight-test data for gusts are not com- 
plete and some compromises must be made to obtain 
the information required in Eq. (10). The variance is 
readily obtained from information on the number of 
gusts within set intervals of gust velocity. This infor- 
mation is usually given. The mean number of occur- 
ences can be obtained from the number of gusts, flight 
hours and mean speed, provided the percentage of 
flight time in gusts is known. Percentages of flight 
distance in turbulence versus altitude are given in ref- 
erence 11. This table can be used in case the percent- 
ages of flight time in gusts are not available from test 
Some duration data on gusts are given in 
For the problem at hand, an analy- 
low-level 


records. 


references 12—14. 


sis of unpublished flight test data over 
mountains puts 
n = 0.001 gust ‘ft., o,? = 292 (ft./sec.)* 


+ 10.666p* + 


32.203p 


£,, and /,,, are taken from measurements of duration of 


ground gusts reported in reference 12. Using the 
preceding values for # and o,*, and 
= @On., 4. = 18. 
from reference 12, and 
1” = 963 ft./sec. 


for airplane velocity, and putting these values in Eq. 


(10) gives 
0.296 (2.8 + 0.00215w" 


(] 


Ww = ae 
+ 0.00269w")? 


[1 — U(w — 96.3 r)] (20 
Eq. (20) is a power spectrum equation for low-level 
gust over continuous mountains when flight velocity is 
963 ft. sec. 
The operations required in Eqs. (14) and (17) 
by numerical integration, 


are 


conveniently handled re- 


placing the integrals in Eqs. (14) by 


) 2 a; T(w, )? W (a, } 


kA 


- a; T(a,)? 


h 


oe 


A? 


Ww, ey)” 


tan 


~ 
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Fic. 3. Smoothing of data-reduction procedures in sample of a 
time trace. 
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Fic. 4. Output spectrum of c.g. incremental gust load factors of 
an airplane in flight through random gusts. 
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Fic. 2. Probability-density function for the duration of gusts 
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and Eq. (2')) for low level gusts 
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FREQUENCY 
Fic, 8. Output spectrum of wing-tip normal acceleration of a0 
airplane in flight through random gusts. 
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POWER-SPECTRUM EQUATION 


TABLE 1 
Expect g. Incremental Gust Load Factors, V = 963 ft./sec 
( 1; Ani, g’s 
C N( An > An sec 
( 3; no./100 hours, 40% flight time in gusts 
( $: no./1,.000 miles, 40% flight time in gusts 
2 3 j 
3.831 551,657 20 993 
2 3.115 $48 528 17,069 
2.206 317,682 12,089 
1.361 196,008 7,459 
0.732 105,351 $009 
( 0.342 19 326 1,877 
7 0.140 20,119 765 
s 0.050 7,148 272 
4) Qa 0.015 ye 4b S4 
0.004 596 99 
0.001 140 5 
L.2 28 l 
1.3 5 
{ 0.8 


where the a,.s are suitable integrating elements. Using 
parabolic row integrating elements and frequency in- 
tervals of 0.1 rad. on the frequency interval (0,6), 
0.2 rad. on (6,50) and 1 rad. on (50,96) gives the nu- 


merical results 
A,? = 0.07248, I,4, = 48.207 


Substituting these values in Eq. (17) gives the results 
tabulated in Table 1. 

The operations described above were performed on 
high-speed digital computing equipment. 

In Table 1, column 1 is expected incremental load 
factor, A,,;, Column 2 is the number of times per second 
A, is > A 
when 40 per cent flight time is in low-level gusts, and 


. Column 3 is this information for 100 hours 


Column 4 is this information for 1,000 miles of flight. 


Flight speed is 963 ft. sec. 


OTHER FLIGHT CONDITIONS 


Statistical parameters for other flight conditions and 
altitudes are listed in Table 2. These values are from 
an analysis of unpublished data and references cited 


in the preceding section. 


EXTREME GUST SPECTRA 


Aside from the use of test records, it is easy to con- 
Struct probability-density distributions from gust 
velocities, select expected and extreme values for the 
other parameters for use in Eq. (10), and in this way 
obtain extreme-spectra equations for clear-sky and 
thunderstorm gusts. 

To do this for clear-sky gusts, construct a prob- 
ability-density function for a velocity distribution on 
Say that 
all velocities are equally likely to occur. Compute 


interval (—50,—2) and (2,60), e.g., Fig. 5. 

1,000 ft.°/sec.? from p(v). For low-level gusts 
use £ 60 ft. For gusts at other altitudes use £);, = 
300 ft. Take /,, = 10 ft. and # = 675,280, no. ‘ft. 
Put Eq. (10) in terms of the transformed variable 


“ = w I’ to get rid of velocity. Then, 
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1,105 (2.8 + 2,000 2° 


(1 + 2,500 Q?)- 


W(Q) = 


31,500 (2.965 + 81,200 Q? 


(1 + 84,100 Q*)- 


[ele] 


Eqs. (21) and (22) give extreme spectra for clear-sky 
gusts, Eq. (21) is for low-level gusts, and Eq. (22) is for 
other altitudes. 

For thunderstorm gusts, eliminate the small veloci 
ties in the preceding distribution. Make the velocities 
equally likely to occur on intervals (—50,—20) and 
(20,60)—e.s 1,500 ft.* sec. 
from this distribution. Use thunderstorm'* data for 
the mode of duration and put &,, = 130 ft. Take / 

10 ft. and 2 = 6.7/5,280, no./ft. Then, 


Compute o,° 


y Fi yr ) 
y., Fig. 0. 


(2.9166 + 13,200 2 
W(Q2) = 9,470 
(1 + 14,400 Q-)- 


Pee i)] & 


Spectra reduced from flight-test measurements sum 


marized in reference 9 and spectra from Eq. (20) are 
shown in Fig. 7. 


REMARKS 


It is common practice to design for gust load factors 
on the basis of a simplified representation of an aircraft 
that is free in plunge in flight through a discrete gust. 
Of course, it is recognized that in the past these simpli- 
fications were necessary because fast computing ma- 
chines were not available and advances in the mathe- 
matics of the analysis of random processes were not 
developed enough to use. But these conditions are 
not present now since fast computing machines are 


available and adequate methods of analysis of random 


) 


TABLE 2 
List of Statistical Parameters* 


Col. 1; Interval of altitude, multiply by 1,000 
Col. 2; 7, no. /5,280 ft 

Col. 3; &., mode of duration, ft 

Col. 4; o,2, variance of gust velocity, ft.2/sec.* 


(a) Clear-sky 
17 > } } 
Ground 5.3 60 292 
0-10 3.60 300 20) 
10-20 3. 4 300 107 
20-30 Re. 300 178 
30—40 9 5 300 57 
40-50) 1.4 300 113 
(b) Thunderstorm 
0-20 6.7 130 324 
*f and o,* for clear-sky gusts from analysis of unpublished 
data; £:, = 300, from reference 12, assuming thickness and width 
of turbulence areas distributed the same; £:, = 60, from reference 


13. Uselms = 10 ft. for both clear-sky and thunderstorm gusts 
7 In 1,000-ft. intervals 
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to the Buffeting Problem, Journal of the Aeronautical Sciences 


processes are known for a realistic analysis for the gust 
Vol. 19, No. 12, pp. 793-800, 822, December, 1952 


problem. Therefore it is no longer necessary to stick 


is as i ; : ; ae ®§ Liepmann, H. W., Extension of the Statistical Approach to t} 
to the previous simplifications. As a matter of fact, Buffeting and Gust Response of Wings of Finite Span, Journ 
the transition stage is already here since enlightened of the Aeronautical Sciences, Vol. 22, No. 3, pp. 197-200, March 
design decisions can be made on the basis of informa- 1955 
tion from both methods of analysis. ‘Miles, J. W., An Approach to the Buffeting of Aircraft Stru 


tures by Jets, Journal of the Aeronautical Sciences, Vol 21, N 


However use of the random-gust method of analysis cgi 
11, pp. 753-762, November, 1954 


for design purposes is limited because of the lack of 


. . 3 Luskin, Harold, and Lapin, Ellis, 4n Analytical Af Ich t 
critical design spectra. These criteria need to be the Fuel Sloshing and Buffeting Problems of Aircraft, Journal of 
established. Also it is necessary to set requirements the Aeronautical Sciences, Vol. 19, No. 4, pp. 217-228, April 
on system representation. System representation re- 1952 
quires the inclusion of airframe modes of vibration. ; 'Brem, Hersy, and Howkelt, — codrapeade PP : 
This need is demonstrated in Fig. 8, a finite output Come oenencetnana: se edheagtednggciny pu Aarne, Sm 

oe : ee ; : nal of the Aeronautical Sciences, Vol. 22, No. 1, pp. 17-2 
spectrum of wing-tip normal accelerations in flight January, 1955 
through gusts that caused a maximum c.g. incremental 1 yon Karman, T., and Howarth, L., On the Statistical Th 
load factor of 0.75g. Airframe-coupled modes through of Isotropic Turbulence, Proceedings of the Royal Society, Series 


A, Vol. 164, pp. 192-201, 1938 
11 Copp, M. R., and Coleman, T. L., Jnfluence of Flight Plan on 
Load Histories and Riding Comfort of Transport Airlines, Langley 
REFERENCES Aeronautical Laboratory, NACA Conference on Some Probl ms 
Bell of Aircraft Operations, Lewis Flight Propulsion Laboratory, pp 
8-8.10, 1954. 


fuselage first-bending are excited. 


1 Rice, S. O., Mathematical Analysis of Random Noise, 
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Oscillating Shock Boundary-Laver Interaction’ 


LEON TRILLING* 


Massachuselts Institute 


SUMMARY 


\ simple model of the interaction of an oscillating shock wave 


with a laminar boundary layer is proposed It involves a small 
perturbation of the steady interaction pattern found by Chap 
man, Larsen, and Kuehn! and Hakkinen, Greber, Trilling, and 
\barbanel.?. It is found that for any local Mach Number above 
1.5 and for any Reynolds Number there is a combination of fre- 
quency and shock strength for which shock oscillations are neu- 


trally stable In particular, normal shocks reach the stability 
threshold when the local upstream Mach Number is between 1.57 
ind 1.65, depending on the Reynolds Number. <A viscous inter 
action mechanism of this type may play a role in exciting tran- 


sonic buffeting and buzz 


SYMBOLS 


skin-friction coefficient 
Ae = change in skin-friction coefficient from flat plate value 
pressure coefficient 


( pressure-coefficient rise from uudisturbed to plateau 


value 
mY change in pressure coeflicient from flat plate value 
A pressure-coefficient rise across shock 
h normal distance from solid surface to inner edge of 
separated boundary layer 
h maximum mean height of separation region 
k reduced frequency based on free-stream velocity amd 


distance from origin of boundary layer to shock 
position 

K reduced frequency based on backflow velocity and 
length of separated region 

kK = eigenvalues of reduced frequency A 


length of separated region 


L distance from origin of boundary laver to mean 
shock position 
VW mean flow Mach Number 
N exponent of the viscosity-temperature relation 
recovery factor 
R Reynolds Number based on free-stream velocity and 


distance from boundary-layer origin to shock posi 
tion 
ratio of free-stream density to gas density at the wall 
= time variable 
] local temperature 
reference temperature 
velocity (variable 
free-stream velocity 
average backflow velocity in the separated region 
length coordinate measured along solid surface 
\ position of moving shock measured with respect to 
its mean position 
a inclination of incident shock 
ry) Mach parameter 6 = (./? — 1)!/? 
ratio of specific heats 
boundary-layer displacement thickness 
e increase in boundary-layer thickness due to pressure 


rise 
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of Technology 


0 angle between separated boundary layer and solid 
surface 
r ratio of backflow velocity in separated region to 


free-stream velocity 


rv = viscosity of fluid 

Mu = viscosity at the wall 

bu viscosity at a reference temperature 
p gas density 

p = gas density at the wall 

w = frequency of shock oscillation 


INTRODUCTION 


ik IS WELL KNOWN that most shock waves observed 
in wind tunnels and in flight are not steady, but 
oscillate slightly about a mean position. In the case 
of oblique waves whose position is anchored by the 
flow geometry (angle, generating wedge, etc.) these 
oscillations remain of small amplitude, but the normal 
shock waves which occur on transonic wings and control 
surfaces are observed to oscillate with large amplitude 
and their motion causes severe buffeting and control 
The mechanism by which 
off is not 


stood; it may involve the effect of pressure disturbances 


problems in transonice flight. 


these shock oscillations are set vet under 
traveling upstream from a trailing edge or other flow 
irregularities. In the present paper, it 1s shown that 
the interaction of a shock with a laminar boundary 
layer is also capable of originating self-supported (eigen 
value) oscillations and may therefore contribute to the 
initiation of transonic buffeting and buzz. 


DISCUSSION 


The mechanism of the interaction of a laminar bound- 
ary layer with a shock wave of finite strength has been 
described recently by Chapman, Larsen, and Kuehn! 
and by Hakkinen, Greber, Trilling, and Abarbanel.” 
It operates essentially as follows (see Fig. | If the 
shock is strong enough to separate the boundary layer, 
the high pressure downstream of the shock pushes some 
fluid 
which is lifted to form a free jet, but preserves its inner 
If the separation 


forward under the separated boundary layer 
structure almost without change. 
point is situated several boundary-layer thicknesses 
ahead of the shock impact point, it is found that the 
separation phenomenon is essentially independent of 
the shock structure, and is self-induced and self-similar. 

In fact, the separation conditions are obtained by 
stating that the pressure rise associated with separation 
is proportional to the rate of change of displacement 


thickness 


Ac, (2B) (de dx) (] 


/ 
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Fic. 1. Sketch of shock laminar boundary-layer interaction 


Furthermore, if the boundary layer thickens by back- 
flow of low momentum fluid near the wall, the momen- 
tum thickness does not change appreciably so that the 
momentum balance relation becomes 


— Ac, Ae(dc,/dx) (2) 
where Ac, is the difference between the skin-friction co- 
efficient at the wall and at the lower edge of the lifted 
boundary layer. 

Elimination of « from Eqs. (1) and (2) yields the fol- 
lowing result for the pressure rise associated with self- 
induced boundary-layer separation: 

Ac, = V 2Ac,/6 (3) 

This result is derived and discussed in detail in refer- 
ence 1, pp. 35-37 and compared to measurements in 
Fig. 34 of the same reference. The pressure rise up to 
separation involves a shear change Ac, cy where cy 
is the boundary-layer skin-friction coefficient observed 
in the absence of the shock. After separation the 
pressure continues to rise as the backflow builds up to 
a limiting intensity, and the skin friction reaches a 
value approximately equal to —c,; then the pressure 
ceases to rise and reaches a plateau, and the separated 
boundary layer forms a straight jet inclined to the 
stream by an angle 6, under which flows a uniform 
backflow of velocity 7. 

The height of the separated region depends on the 
mass of fluid forced into the separated backflow region 
by the overpressure due to the shock. Integrating the 
momentum equation from A to B (see Fig. 1) one finds 


2B eB 
| [puu, + pvu, + pr|dx = | (uu,),dx 
a A x 1 

But along AB, wu vanishes and experiments?” indicate that 
(uu,),issmall. Therefore, 


°B 
(u,)ave. | (— pv)dx _ Heseve. Mw | (1 2) p. a” x 


A 
Putty = (1/2) p.U~2Ac, 


Ps) 
or pu( ule Mw) _ Acs. (Cyr)ave. 


In the separated region, momentum is destroyed by the 


action of viscous forces alone, at the rate p,@*dh/dx = 
p»t"O per unit length, so that 
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ato = p.U 6; 


The parameters of the separated flow region are, there 


fore, 


P( 110 py) Ac Cy 5 


Y sh 


where Ac, | 


equals the shock strength less the pressure rises needed 


to separate and reattach the boundary layer. Som: 
manipulation and the introduction of the law which ri 
lates cy, to the distance Z from shock to leading edg 


then gives the approximate formula 


L/L = 2e,/B)*!4s* tO!” Ac. 6; 


sh 


where » is the exponent of the viscosity-temperatur: 
relation n/#. = (7 /T.)"*. 

Eqs. (4)-(6)* describe the steady shock boundary 
layer interaction mechanism and are described in detail 
and compared to measurements in reference 2. In the 
following pages, an attempt is made to examine the 
flow pattern which results from small harmonic oscilla- 
tions of the shock wave, both in strength and position 

Changes in the position of the shock induce changes 
in the position of the separation point and therefore in 
, a, and//L. Changes in the shock strength in 
duce changes in mass flow and, therefore, essentially in 
/ Land @. 
rated boundary layer which acts as an almost imper- 
meable perfectly flexible membrane and _ separates a 
region of slow viscous backflow near the wall from a1 
Any distor 


CH€ 


Py] 


Both of these modify the shape of the sepa 


outer region of supersonic potential flow. 
tion of this membrane induces pressure disturbances in 
the potential flow which affect the motion and strength 
of the shock. This model suggests the possibility that 
under some conditions the original shock perturbations 
and the shock perturbations induced by the wiggling oi 
the separated boundary layer so interact that a neu 
trally stable eigenvalue problem results. 

A simplified model of the oscillating perturbation 1s 
now constructed. The principal feature of the dis 
turbed flow is the distortion of the separated boundary 
layer under the action of variable backflow mass input 
It is assumed that backflow mass fluctuations travel 
upstream at the mean backflow velocity without dis- 
tortion so that the height of the separated boundary) 
layer above solid wall at any point and at any time 1s 


given by an equation of the form 


h(x, t) = 0(8/1) + het FF) = Oy Hb Ige’*®O 


where x is the distance from the mean position of the 
separation point (made dimensionless by division 
* The numerical values of the constants in Eqs. (4)—(6) ane 
the form of Eq. (6) are slightly different from those of references 
2, mainly to simplify the subsequent analysis at little cost ™ 


accuracy. 









is the pressure rise due to the shock, and 


or 
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through /) and A is the reduced frequency based on the where c,’(1+) denotes the pressure perturbation im 


re, there | backflow velocity. Note that the reduced frequency mediately behind the shock and c,’ the perturbation in 
| based on the overall parameters L, U.. is obtained mean skin friction over the boundary layer between 
from Eqs. (4)-(6) as shock and reattachment. The pressure behind the 
b oL/U K /2s"Ac. (R shock fluctuates for two reasons—the pressure level 
- sa | wh ahead of it fluctuates according to Eq. (9) and the in- 
When the boundary of the potential flow oscillates in coming Mach Number fluctuates since the shock moves 
) the manner suggested by Eqs. (7), it creates an oscil- at a velocity ¥ = iAXe'*. 
nck. and lating pressure field in the outer potential stream given Thus one has, for the pressure perturbation immedi 
by (see reference 3, p. 357): ately behind the shock, 
S needed j 
Some |  ¢,(x, f (2iKhy/B) A + ATRET” x c, (14 —c,'(1—) — (XN /u.) (dae, dM 14) 
vhich re- | { o— FAMERS (MP— DINK Mx/(M? — 1)] + ihn 
Ing edgy ax 
i1+N)K | etl \M?2/(M?—1)]}Ké x c,(1—) é [27K ho(1 x B lek {id K/(M?-1 y 
6 Ji (\K ME (Mf? — 4 dé} (9) Ji AMK (W2?—1)]/ + 21+A)A X 
=} 
perature where A ul’. so that AA is the reduced frequency eFC ME OURS FINK Mt/(M? — 1)) 2} 15a) 
based on free-stream velocity L’., and length of sepa- sth 
undary rated flow /. (xX /T (dAc, , dA) IAA X(d/dAM) } {A (y + I 4 
n detail The motion of the shock is given by X« ‘Kt ond the [cos? a — (1/M?)}} = [SAAN /(y + IM L5b) 
In the motion of the separation point is then the sum of the a Weare ; 
ine the shock travel and the change in length of the separated The skin-friction disturbance 1s due to the fact that the 
oscilla- region due to deformation of the separated boundary region behind the shock ee with the shock motion 
osition layer, |— (Ao aet**]. If the boundary layer reacts in- ane He panei 4 ung. — ™ seo “a rnd 
changes stantaneously to a change in conditions (and this is a tuation is needed, it Is expected that it will have a time 
fore in reasonable assumption since its characteristic time is lag with respect to X. The details of the flow ewe 
eth in s/U., as compared to 1/i for the separated flow, and / ture behind the shock have not yet been determined so 
— is 80 to 150 boundary-layer thicknesses while @ is that the tune ing cannot be estimated very accurately ; 
e sepa 1/5 to 1/3 of U.), then the pressure perturbation at sta order e’* and has so been assumed in the analysis. 
imper- the origin is obtained by perturbing Eq. (4) as follows rhus, 
ates a from now on the multiplier e*! will be dropped from ce /¢; = —(X/2L)e* (16) 
‘om an all perturbation terms) : ; Lo See - . 
Kaien Substitution of Eqs. (7), (10b), (12), (14), (15a 
eat Cy,/Cp, (1/2) (Ac,/cy) = —(ANnsep./4L) (15b), and (16) into boundary condition (15) gives a 
vent ~[(X = ho/)/4L] (10a) relation which the flow parameters .1/, Ac,, A, \ must 
v that satisfy if there are to be neutral self-sustained oscilla 
witliaaie or, solving for Y and substituting Eqs. (4), (6), and tions: 
Sateal )) to rewrite Cc, and ¢, . il = As*dc, +) + [5iK(L + 2)/2 1 ik 
A neu- 0X ho =1- [27A (1 + d) As"Ac, a! (10b) [SA M4 aa Ac, ’ 1 [—2A(1 4. > 
1001 1s The backflow velocity perturbation at v = 0 is also ob- s"Acp,,] —idy — [20KA*(1 + 2) ‘sc, 
e dis- tained from a perturbation of Eq. (4) as follows: { Jo[AMK /(M? — 1)] x &RE OMEN 4 
ndar\ : i i - 
ou u(0)/tiw = (1/2) (Cp /€p,, (11) iK(L + eT. — ih (1 
_ lf the backflow velocity perturbation is transported aeeeis 
a rate the particle velocity, then from Eq. (11) it follows Lil [ on HIF IAME/ (MPD) KE yg 
me is e ; 
6u/hot = [iK( + d)/2] ™ (12) Jo [PAKME (ME — 1) \d— 1S 
It remains to discuss the conditions at the shock (x = Eq. (17) defines a complex eigenvalue problem for the 
= |). If the shock motion is specified by a physical mech- frequency x oo te re wae & an steele 
defined by its Mach Number and Reynolds Number, 


anism such as a moving shock generator, then X is ; ay : ; 

, <7 - : . ‘ two equations obtained by equating the real anc 

f the specified and the problem is determined. However, a the : punt ins <i , ee : 
imaginary parts of Eq. (17) to zero give a series ol char- 


more interesting problem consists in requiring merely 


{sion : 
, . . ve - 6 acteristic frequencies and a critical shock ngth fo 
that the Rankine Hugoniot conditions be satisfied icteristic frequencies and a critical shock strength for 
nd TOS ; -ach frequency 

and across the shock, and that the rate of backflow be given cack irequency. , 
‘T nce a . ‘ = = y > c > > Oa > i T > < »<¢ ‘ > - 
pe by a perturbation of Eq. (5) Note that the Reynolds Number ippears in the prob 
wate lem only through the parameter \ which is related to 


u/u) + (h/6) = [c,’(1 +)/Ac, 


ray) — (es'/e) (18) R by 





= 0905 
[C72 — 1)R,]! ») 


Vv S(Cy B)! 4 
— 1)/2)rM} i”? 


A= u/U,= 


(1 |. [(y (19) 


The Reynolds Number dependence is therefore weak; 


in fact, the parameter \ goes from zero to 0.5 when R 
goes from infinity to 100,000 at 17 = 1.5. The Mach 
Number ahead of the shock is seldom below 1.5 and R 
is in most cases appreciably above 10°. The significant 
range of the parameter d is therefore between 0 and 0.5, 
and in the approximate solution of Eq. (17) several 
terms are expanded in powers of \ and only three terms 
are retained. 
In particular, the integrals /,, 7, may be evaluated 
as follows: 
IT, = {1/K(1 + A)] sin A + O(A?); ) 
I, = [1/K(1 + \)] (1 — cos K) + O(A2)f 


(20) 

Separation of real and imaginary parts in Eq. (17) then 

gives 

cos K(1 — As"Ac,.) — sin Kj (5/2)K(1 +2) + 
[4KA*(1 + d)/sAc, Jf = —[16K? X 


(1 + A)/(y + 1)Mis*dc, 2] +... (Zila 
sin K(1 — As"Ac,) + cos K}(5,2)K(1 + d) + 
[4KA*(1 + d)/sAcy ]f = (2KA/ Ac, ,) X 
{f2n(1 + A)/s] — [8/(y + IMI} +... (21b) 


Elimination of the trigonometric functions from Eqs. 
(21a) and (21b) then gives the following algebraic rela- 
tion for Ac, (A, Oe De 
[16A7(1 + A)/(y + 1 ) M's" Ac, *]? si 
(4K7A?/Ac, 7) {[2A(1 + A)/s] — 

+ 1)\M/3)}? = 


(1 — As"Ac,.)? 
(5K 2) (1 +2) + HK + A) sAcy, |}? (22) 


' 


[S (¥ 


In the special case when \ ~ 0, R > @& the solution 


of Eq. (22) is found as 


Ac, V 32K D(y + 1)1/*s" [1 + 0(1/12K")] (23) 


and the characteristic frequency & satisfies the equation 


AK, 1.79 K; = 
K» = 1.SO 


—(5/2)K 7.58] 


tan. A = 
‘ (24) 
a { 


In all cases, } [(5/2)A |? + 1) '/? may reasonably be re- 


( 
‘ 
placed by 5 2 A and in fact, Eq. (25) gives the pressure 


rise up to terms in A* with reasonable accuracy. The 
frequency for finite values of \ is given by 
tan A 

—(5/2)A[1 + ACL — 3V1.2s"/M*)] + O(A2) = (25) 


and the total pressure rise is found by adding to Ac, | 
the plateau pressure and reattachment pressure rise 


V 32K /5(y + 1) As" + 
[32/25(y + 1) AF] + (38A2/s) +... 


Ac, = 


} 
(26) 


Eqs. (25) and (26) give the shock strengths and fre- 


quencies at which neutral self-excited oscillations are 
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O 
Fic. 2. Intensity of the neutrally stable oscillating shock wav 
as a function of upstream Mach Number 
possible. The critical frequency is seen to be rather 


insensitive to either Mach Number or Reynolds Num 
ber changes over a wide range of values of these param 
eters. The critical shock strength is plotted as a func 
tion of Mach Number in Fig. 2 for several values of 
the Reynolds Number parameter X. 
The curves Ac, (.1/) intersect the Ac, (-V 


the curve which gives the pressure coefficient jump 


curve 


across a normal shock Ac, [4/(y + 1)] [l - 
(1/M?) ]—in M This 
means that at Mach Numbers below 1.6 the oscillation 


norm. 
the neighborhood of 1.6. 
is always stable; normal shocks are unstable near .\/ 
1.57 (A 0) to 1.65 (A 0.3, R ~ 1.2 X 10°). At 
higher Mach Numbers the critical shock is much weaker 
than a normal shock. 
critical shock strength corresponds to a flow deflection 
and at l/ 3 to a flow deflection of 4 


For example at 1/ ?, the 
of 6 


CONCLUSION 


It has been shown that the interaction of an oscillat- 
ing shock with a laminar boundary layer gives rise to 
an eigenvalue problem and that for certain values of 
shock strength and frequency, self-supported oscilla 
tions of the flow pattern occur. This mechanism or atl 
analogous mechanism involving a turbulent boundary 
layer may contribute to the initiation of transonic buzz 


and buffeting. 
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Transformation of the Compressible Turbulent 
Boundary Laver' 


ARTUR MAGER* 
Marquardt Aircraft Company 


SUMMARY 


nsformation of the compressible turbulent boundary 
demon 


The tra 
ver equations to their 


rhe transformation is essentially the same 


incompressible equivalent is 


strated analytically 
that for the laminar layer, first given by Stewartson, except 


as thi 
hat the explicit relation between the viscosity and temperature 


is not required \ key point in the analysis is the modification 


if the stream function to include a mean of the fluctuating com 


ind the postulate that the apparent turbulent shear, 


ponents 
ponent 


issociated with an elemental mass, remains invariant in the 


transformation 

The values of the incompressible friction coefficients and of 
pressure rise good 
ugreement with the experimentally measured and independently 
\n application of the transformation to the 


layers and to the 


causing separation thus transformed show 
reported results 
self-preserving boundary computations of 
general boundary-layer flow is shown 


SYMBOLS 


= speed of sound 
= specific heat 


= friction coefficient 


F, ( = functions of x 
= Maskell’s functions of '* and H* 
H = form parameter (= A/6 
= mean total enthalpy 
= modified total enthalpy 
R, Ri, k = constants 
5 = characteristic length 
VW = Mach Number at y = 6 
m = function of 1/ 
p = mean pressure 
= resultant velocity 
R = gas constant 
R = Reynolds Number at y = 6 and based on 
streamwise length 
= mean temperature 
= it at X¥ = 6 
u, = mean velocities 
X, ] = transformed coordinates 
x, ¥, = real coordinates 
a = friction coefficients ratio 
I = pressure gradient parameter 
= specific heats ratio 
A = displacement thickness 
6 = boundary-layer thickness 
0,0 = generalized momentum thicknesses 
6 = momentum thickness 
K = thermal conductivity 
mM = absolute viscosity 
; = kinematic viscosity 
p = mean density 
a = Prandtl Number 


shear 
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= stream function 
rature-viscosity relation 


= exponent in tem] 


Subscripts and Superscripts 
=aty = 6 
= final 
= initial 
incompressible at same Re 
stagnation 
( = virtual origin 
separation 
coordinate subscripts designate partial differentiation 
: = transformed 
= fluctuating value 
mean value 


INTRODUCTION 


KF R LAMINAR FLOW, the computation of the compres 
sible boundary-layer development is greatly facili 
tated by the use of a transformation which establishes 
This trans 
and later 


the corresponding incompressible flow. 
formation, first suggested by Dorodnitsyn 
modified by Howarth,’ Illingworth,* and Stewartson,* 
is extremely useful since it enables one to interpret all 
the hitherto available knowledge of the incompressible 
laminar boundary layer in terms of its compressible 
equivalent. The transformation from the compressible 
to the incompressible form is accomplished, in strict 
sense, only when the Prandtl Number is unity and there 
is no heat transfer between the wall and the fluid. 
However, even in cases where the equations do not trans 
form exactly, the use of the transformation eliminates 
the varying density from most of the terms and thus 
simplifies considerably the solution of the problem. 
In this way, for example, it was possible to apply the 
Thwaites method® to the compressible laminar bound 
ary-layer flow with heat transfer.' 

It has been shown, independently, by Dorodnitsyn’ 
and by Van Le’ that the transformation valid in the 
laminar boundary layer also eliminates the varying 
mean density from the momentum-integral equations 
of the turbulent boundary layer. This, of course, is 
not surprising, since in the integrated form, certain 
terms involving the mean of the product of the fluctu 
ating density and velocities are neglected, with the 
justification that, while these terms can possibly be 
large locally, when integrated over the whole boundary 
Con- 


sequently, the momentum-integral equations for the 


layer they will tend to cancel themselves out.’ 


compressible turbulent boundary layer are identical 
in form to those for the laminar layer, and one may 


expect the same transformation to apply. 
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In spite of this fact, however, one may not be con- 
fident that the laminar transformation is actually valid 
for the turbulent flow. 
be raised are as follows: 
verts the once-integrated form of the equations instead 


The main objections that can 
(a) the transformation con- 


of the differential ones and consequently does not 
ensure that the effect of the terms eliminated by aver- 
aging is not lost; (b) the transformation requires a cer- 
tain change in the wall shear, which in view of our 
limited knowledge of the turbulent shear, may or may 
not be correct; and (c) the transformation (at least in 
reference 8) requires a redefinition of the displacement 
thickness, which is not convenient conceptually.* 

The goal of this paper is to overcome these objections 
so that the application of the transformation to the 
turbulent boundary layer can be justifiable on more 
solid theoretical grounds and thus more acceptable. 
Furthermore, a comparison between the theoretical 
and experimental results for the skin friction of flat 
plates as well as for the separation pressure ratio of the 
turbulent layer will be made. Finally, a conversion 
of the quantities essential to the computation of the 
general compressible turbulent boundary-layer devel- 
opment will be demonstrated. 
TRANSFORMATION OF THE TURBULENT BOUNDARY 
LAYER 


The equations of the steady compressible turbulent 
boundary layer were derived by Young." This deri- 
vation requires, in addition to the conventional assump- 
tions, that the partial derivatives with respect to x 
of the mean stresses due to fluctuations be neglected, 
as small, when compared with the corresponding de- 
rivatives with respect to y. Furthermore, a number of 
terms involving triple correlations, the fluctuating 
thermal conductivity, and the dilatation were also 
neglected. When the mean values are taken, the 

* Since the presentation of this paper and as this manuscript 
goes to press, author’s attention has been called to three additional 
contributions to the problem of compressible turbulent boundary- 
layer transformation. One of these is a note by R. E. Meyer, 
Turbulent Boundary Layers on Nozzle Liners (Journal of the 
Aeronautical Sciences, Vol. 22, No. 8, pp. 572-573, August, 1955), 
in which the compressible flow transformation of the momentum 
integral equation is postulated and then the thus computed re- 
sults are shown to agree with those obtained from empirically 
derived formula. The other two papers are by G. W. Englert, 
Estimation of Compressible Boundary Layer Growth Over Insulated 
Surfaces With Pressure Gradient (NACA TN 4022, June, 1957), 
and by F. E. C. Culick and J. A. F. Hill, A Turbulent Analeg of 
the Stewartson-Illingworth Transformation (Journal of the Aero 
nautical Sciences, Vol. 25, No. 4, pp. 259-262, April, 1958). In 
both of these papers it is pointed out, as already noted by Van Le, 
that the possibility of transforming the momentum-integral 
equation hinges on the relation between the compressible and 
incompressible friction laws. The authors then substitute the ex 
perimentally obtained friction laws, involving a certain reference 
temperature, and show that the resulting transformation is very 
nearly equivalent to that of Stewartson. Although such a pro 
cedure weakens some of the objections raised above, still all of 
these treatments do not ensure that the effect of the terms 
eliminated by averaging in the original momentum-integral 


equation is not lost 
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equations become 


(pt). + (pv + pv’), = 0 ( 


puu, + (pv + p’v’)u, = —pr + (uu, — pu'v’),] 


puh, + (pv + p'v’)h, = 

) (Kk Ca) \h + 0.5(o — I )q?| — pu’h’| 5) 
In these equations the primes denote the fluctuating 
terms and /, the total enthalpy, includes the mean 
kinetic energy associated with the turbulent velocity 
components—e.g., h = (c,t + 0.5 q?), where g is the 
magnitude of the resultant velocity. It is quite ob 
vious from Eq. (3) that one solution of this equation 
when there is no heat transfer and when o = 1.0. js 
given by / = constant and v’h’ = 0. 

The essence of the laminar boundary-layer trans 
formation lies in the stretching of the physical coordi- 
nates into an imaginary incompressible system (with 
out varying density) in which the properties of the fluid 
are evaluated at certain reference conditions. Jn 
addition, in order to ensure that the two boundary 
layers are really corresponding, the stream function, 
which is the measure of mass flow contained in the 
boundary layer, remains unchanged by the transforma- 
tion. 

In the turbulent boundary layer, however, where the 
stream function does not describe the boundary layer 
completely, because of the presence of the turbulent 
stress, one may expect additional requirements which 
will ascertain that the transformed and the original 
layers really do correspond. In particular, if one as 
sumes that at any point in the flow field the turbulent 
shear is a property of the flow which is essentially con- 
stant over an infinitesimal mass of the fluid, then, in 
view of the invariance of the stream function, it seems 
logical to require that the total amount of shear asso- 
ciated with such elemental mass should also remain 
invariant. From a physical standpoint, this means 
that the mechanism of the turbulent process is un- 
changed by the compressibility effect ; however, the fluid 
is allowed to expand or contract in response to the in 
ternal heating.t In addition to this specific assump- 
tion on the mode of transformation of the turbulent 
stress, one also should choose the reference conditions 
at the stagnation value so that the total enthalpy (for 
o = |.0 and adiabatic flow) is the same. This choice 
then ensures that the internal heating affecting the 
turbulent stress transformation is properly scaled. It 
should be noted that this is somewhat different from the 
laminar transformation, where the choice of the refer- 
ence condition is essentially unimportant.} It is these 

+ This meaning has been pointed out to the author by Prof 
Hans W. Liepmann. 

t This has been called to the author’s attention by Prof. Nich 
las Rott during discussion at the 1957 Heat Transfer and Fluid 


Mechanics Institute 
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requirements on the turbulent stress transformation 
that will be used in the present paper. 

First, a transformation, similar to that used for the 
laminar boundary layer, may be postulated: 


Y= FG dx; 


0 dx = FG(O OX) + ¥.(0 OY (4 


Y=F (p/pu)dy; O/Oy = 


F(p/ po) (0 OY 


with the additional requirements that 


pdX dY¥ dz (5) 


. I *- 


y y (pu ‘y’ \ pdx dy dz = ( pult ‘y'* 


In these equations, X, }, and the quantities denoted 
by the star are understood to be in the transformed 
system, while the functions F(x) and G(x) are as yet 
undetermined. 

The original velocities may now be expressed in terms 
of the stream function, so that the continuity condition 


1) is identically satisfied ; 
pu = pol,; —(pv + pv’) = pov: (6) 


and the transformed velocities should be defined in 


the conventional incompressible manner; 
~~ * - 
— y X \é) 


stream function |Eq. 
relation between the 
Eq. (1) transforms to 


u* = yy*, —y* 


In view of the condition on the 
5)|, Eqs. (6) and (7) furnish a 
velocities in the two systems and 


its incompressible form: 


u = Fu*: pou t+ p’v’ = po(FGv* — Y,u* (S) 


(la) 


To proceed further, one first has to determine the 
explicit relation between the turbulent shears in both 
systems. Therefore, using Eqs. (4) and (5) one ob- 
tains the elemental area transformation 


dXdY = [O0(X, Y)/O(x, y) |]dx dy = (p/ po) F°G dx dy 


leading to 


pu'v’ = ina (9) 


F°G pou'v 


With this condition, Eq. (2) may be put, after some 


manipulation, into the form 


dU* [(1/p) + (1 


F) (dF dp)u*| 
u*uxy* + v*uy* = U* : = . 
dX  {(1/p,.)(1 F)(dF dp)U?| 


l 
I pu 
Tr ( F-uy* — F°Gpou'v’* (10) 
} 


poF*G \ po 


and G that 
will obtain the incompressible form when 
Realizing 


One must now so choose the functions F 
Eq. 10 
¢ = 1.0 and when the flow is adiabatic. 
that the total enthalpy is given approximately by 
; (c,t + 0.5 u*), and that this quantity must be 
1.0 and adiabatic flow, one requires 


hs 
a constant for ¢ = 
that F satisfy the differential equation 
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adF/F = (R/2c,) (dp/p) = [ly - 


for then the first term on the right-hand side becomes 
U*(dlU* dX). Incidentally, if the flow is not adiabatic 


then this term becomes L’*(dl'* dX 
The differential equa 


just as in 


the laminar boundary layer.' 


tion (11) must be solved, in view of relations (S), with 


the requirement that at the reference condition F 


1.0. This solution gives, just as in Stewartson’s 
transformation, 
f= (p pi Shi = a, a = (i ft 12) 


The second term on the right-hand side of Eq. (10 
can similarly be converted to the incompressible form 
when it is realized that the product of viscosity and 
density must behave very nearly like the pressure 
that is, it may be considered independent of y. Con 


sequently, it is permissible to set 


G = (pu/ poo) = (pett./ posto) = (u./po) F° ( 


without specifying the explicit dependence of viscosity 
on temperature. It should be noted that if the vis 
cosity is assumed to vary directly in proportion with 
the temperature, then the function G is identical with 
that given by Stewartson, 


(up = Rit t maid G = k(p Pp (l3a 
Using Eqs. (12) and (13), one is able to convert 
Eq. (10) to the desired incompressible form 
u*ux* + v*uy* = 
U*(dU*/dX) + wuyy* — (u’'v’*)y (2a) 


and demonstrate that the application of the thus postu- 
lated transformation to the turbulent boundary layer 
is indeed valid from the analytical standpoint. The 
slight modification of the Stewartson transformation, 
which was here introduced, should permit a closer 
description of the variation of viscosity with temper- 
ature and should, of course, be also applicable in the 


laminar boundary layer. 
EXPERIMENTAL VERIFICATION 


Flat Plate Skin Friction 


The usefulness of the transformation in the turbulent 
flow may be demonstrated by deriving the appropriate 
formulas for the local and average skin-friction coeffi- 
cients of flat plates. This derivation is especially im- 
portant since the comparison of the thus obtained fric 
tion coefficients with experimentally measured values, 
as well as with values derived by entirely different 
This comparison can 
the 


methods, can readily be made. 


thus serve as some measure of transformation's 
validity. 

Now in incompressible flow a modified Blasius for 
mula for the local coefficient is known to be applicable 


up to Xk, = 10° (see reference 11) 


14) 


0.5.c;* = 0.0296 (R,” 


) Cy = 


To convert this formula to the compressible flow, one 
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must first find the transformations between the friction 
coefficients and between the Reynolds Numbers. 
Using Eqs. (S) and (9) one finds that 


Cy* = (uo/medes Re* = (uc/po)? R, (15) 


With these relations the Blasius formula may be written 
as 
0.5 cr = 0.0296 (R,) —°** (a, /s49)°** (I4a) 


The variation of viscosity with temperature can be 
most accurately described by the Sutherland formula. 
Here, however, in order to avoid analytical complica- 
tions, one can take the simplified relation 


(e/mo) = (t./to)” (16) 


so that the ratio of the compressible to incompressible 
skin-friction coefficients, at the same Reynolds Number, 
becomes 


(C,/Cy, in) = (te/to (14b) 


This relation, with w = 0.76, is shown in Fig. 1, to- 
gether with an analogous relation obtained by Eckert’s 
“reference temperature’? method'' and a few experi- 
mental points of Coles'? and Korkegi.'* As may be 
seen from this Figure, despite the assumption of ¢ = 
1.0, the agreement is reasonably good. 

To illustrate that the good agreement is not due to 
the particular expression used for the incompressible 
coefficient [Eq. (14)], one may compute the average 
skin-friction coefficients by a different formula, equally 
acceptable. According to reference 11, the Karman- 
Schoenherr formula is valid up to kX, = 10°: 


0.242 (é,*)~°"” = logy (R,*é,*) (17) 
By the application of Eqs. (15), this expression may be 
restated for compressible flow as 


(1/2 


) (u, bo)! 2 a 
logw(éR.) — logio(uo/u,) (17a) 


0.242 (€;) 


Using Eq. (16) and setting the ratio of the compressible 
to incompressible skin-friction coefficients (at same R,) 
equal to (1a), Eq. (17a) can be written in a more con- 


venient form for computation as 


logio(€y, in R.) fl — V alt, fy)" |] = 
logio [a(to/t,)°| (7b) 


The solution of Eq. (17b) for R, = 10’ and w = 0.76 
is shown in Fig. 2. Also shown in this Figure are the 
analytical results of Van Driest,'* obtained by the simi- 
larity mixing-length method, and the directly measured 
experimental data of Chapman and Kester.'° Again, 
in spite of the assumed o = 1.0, the agreement with the 
experiment and with the Van Driest analysis is very 


good. 


Shock-Separated Turbulent Boundary Layer 


Since the transformation represents a certain corre- 
spondence between the compressible and incompressible 
boundary layers, if this correspondence is real, the trans- 
formation should be able to convert not only the vari- 
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ous analytical results, but also purely empirical rela- 
tions. Indeed, an experimental verification of the 
so-transformed empirical laws could be considered as 
one of the proofs of the transformation’s validity. Ap 
indication that the empirical laws transform properly 
has already been obtained in the previous section, where 
it was shown that the transformed skin-friction rela 
tions, which, after all, do involve some empiricism 
in their original form, agree reasonably well with ex 
periment and with other theories. An even more spec 
tacular illustration of the usefulness of the transforma 
tion and its validity is furnished by the estimation oj 
the pressure rise occurring with a shock-separated 
turbulent boundary layer. 

As is well known, the separation of the incompressible 
turbulent boundary layer, in the presence of very steep 
pressure gradients, has been observed to occur when 
the velocities outside the boundary layer satisfy 


U2* = kU (18 


with kj = 0.55. This relation between the velocities 
can also be argued on the basis of the semiempirical 
equation relating the shape of the velocity profile to the 
pressure gradient and to the shear at the wall. In 
compressible flow, applying Eq. (8), this purely empiri- 


cal relation becomes 
M ‘= kM; (1Sa 


and consequently one is in a position to evaluate the 
pressure rise connected with shock-caused separation. 
The details of such an analysis have been explained in 
reference 16; here we shall only show the final result. 
This is presented in Fig. 3, along with the experimen- 
tally measured separation pressure ratio for a number 
of shock-causing configurations plotted against the 
Mach number before the shock. As can be seen from 
this Figure, (and considering the difficulty with which 
the separation pressure is experimentally determined), 
the agreement between the transformed results and the 
experimental data is fair. 

As is well known, the pressure at which the separation 
occurs is not the highest pressure measured for such an 
interaction. This, of course, is the reason why it is so 
difficult to determine the separation pressure properly 
The highest pressure reached can be determined with 
considerably greater certainty. On the other hand, 
assuming that the additional pressure rise occurs due 
to the centrifugal effects and once more using the trans- 
formation, the ultimate pressure rise can also be found 
analytically and the details of such a computation are 
also given in reference 16. In Fig. 4 are shown the re 
sults of such computations, together with the experi- 
mental data obtained by number of investigators. 
As may be seen, the analytical results appear to be 
confirmed closely by experiment. 

It is important to note that, in any boundary layer, 
both the shear at the wall and the pressure gradient 
play an important role. In a laminar layer there 
seems to be a certain correlation between the two, as 
indicated by the success of the Thwaites method. In 
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a turbulent layer, however, the pressure gradient and 
the shear at the wall appear to play an essentially inde- 
pendent role. This is indicated, for example, by the 
fact that in the absence of pressure gradient, the so- 
called universal velocity distribution is obtained only 
when the effect of the shear at the wall is taken into 
account. 

Now it has been shown in the previous section that 
the compressible flow transformation gives correct 
results in the absence of the pressure gradient, while 
in this section, good results have been obtained for 
flows with strong pressure gradients but negligible wall 
shear. Thus, so to speak, two extremes of the turbu- 
lent boundary-layer problem have been tested. For a 
complete experimental verification, one should also 
check a case when both the wall shear and the pressure 
gradient play an equally important role. Unfortu- 
nately, however, the absence of experimental data for 
such a general case does not permit one to carry this 


out.T 
APPLICATIONS 


Self-Preserving Boundary Layer 


Townsend has recently shown” that when (py) «’v’*) 


is the only component of the turbulent shear, then a 
self-preserving (equilibrium) incompressible turbulent 
boundary-layer flow can occur if 


(U*/U,*) = [((X — Xor)/(X1 — Xor) |” (19) 


where subscripts 7 and or denote the initial point and 
the virtual origin of the flow, respectively, while con- 
stant k. signifies the appropriate pressure gradient. 
Of greatest interest are, of course, the cases when the 
pressure gradient is unfavorable—that is, when X > 
X,, and k; < 0. For such cases, Townsend has shown 
that the smallest possible value of k, is k; = —1 3. 
Applying the transformation and assuming for sim- 
plicity a linear dependence of viscosity on the tem- 
perature, Eq. (19) can be stated as 


M ig ae — ((8y¥—1)/2(7-1)] 4 

uw, LJ (1 +e ur) ax | x 
vi 7 —_ — [(8y— 1)/2(7-))) DX De 

if (: +? 9 : 1") ax | ' (19a) 


and the required distribution of the Mach Number is 


given by 


) 


: er > l [(3¥ 1)/2(4 1}] 
MV" = hs | (: + 7% ur) dx 


(19b) 


where k; is a constant permitting the correct match be- 
tween the pressure gradient and the boundary layer 
that is, ensuring that the flow comes from the appro- 
priate virtual origin. Eq. (19b) can be differentiated, 
rearranged and integrated for y = 1.4, with the initial 
+ The data quoted by Englert and by Culick and Hill give 
some evidence that the transformation is also valid when the wall 


shear and the pressure gradient are equally important. 
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conditions at x = x;, WZ = A,, to give 
mM) — m(M,) = k(x — x 


mM) = WY"}1 + [0.8/(1 + 2k.) ]M2 + os 
[(0.24/(1 + 4k,)] M4 + [0.032/(1 + Gke)| AM - 
[(0.0016/(1 + Ske) |As*I 


The function 31(.1/) is shown in Fig. 5, while in Fig, ¢ 


the distribution of 7/1; vs. (x — x;) ks is plotted. It 
may be seen that for a flow with the same virtual origin 
a much steeper decrease of Mach Number is possible in 
supersonic flow than in subsonic flow. 


General Boundary-Layer Flow 


The compressible flow transformation is particularly 
valuable in the computation of the general boundary 
layer flow. The methods for such a computation in 
incompressible flow are always somewhat empirical 
and an establishment of analogous methods for the 
compressible flow, in a conventional way, would require 
a tremendous amount of experimental information. 
Using the transformation, however, one can immedi- 
ately reinterpret the incompressible flow methods for 
the compressible flow. One such method has been 
given by Maskell'* and it appears to give very satis- 
factory results.'* This method allows to solve for the 
momentum thickness directly from 


(U*4-2 @* L) — (U*4-2@* L) + 


eX /1 
0.01173 | U*4-2 d(X/L 2] 
J Xi/L 


(9*/L) = (@*/L)9-822% (U*L /po) 0-175 


It should be noted that since 0* is only an auxiliary 
parameter, all one has to find is the conversion from 
0*(X) to O*(x). This can be accomplished by assum- 
ing a power law for viscosity 
(M*-*6*/L) = (M**6*/L), + 0.01173 X 
ax /L 
| M4231 + [(y — 1)/2]Ir EOF" d(x L) (2la 
xi/I 
with ky = 0.5(y + 1)/(y — 1). Also, since it can be 
shown that 


6= 41+ [(y — 1) /2]2!* 6* (22 
one has 


OL = 41 + [(y — 1)/2]2t** (0*/L)93225 x 
(Maoh /vo) 7" (2 


In addition to the momentum thickness, one also 
has to compute the form factor /7*. Now the Maskell 
equation for this involves a number of functions, some 
of which are in graphical form. Therefore, instead of 
elaborating on the details, one can indicate the general 
form of the equation and give the appropriate conver- 
sions which will permit to use the graphically plotted 
functions of Maskell. The equation can be written 

0,* (dH*/dX) = f(T*, H*) | 
0,* = 6* (6* U* yi )0. 268 (24 


,* = (0,* U*) (dU* dX) ( 1(¢))0-67817 | 
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\gain, it is only necessary to convert from /7*(X) to 


HT* (3 


und this can be done by writing 


dH* fi + [(y — 1)/2]Ar7- & 
(L/0,*)f(T*, H*) d(x/L) 
o* L 6* L)(W 414+ [(y — 1)/2)17I - 
(24: 
xX (0/L) (aol. /v ) nace 
r* ‘1 + [(4 — 1)/2] ayer? x 
0,*/L) (1/M) [dM /d(x/L)] (10)9-8#* 
Eqs. (24a) can be numerically integrated knowing 
W(x), O*(x), and @(x). Once H*(x) is known, the dis- 


placement thickness may be found by using the ex- 


pression suggested in reference 6, 
AL=6L (UI*}1 + [(y — 1) 2J4P7} + 
2 |Z?) (25) 


Finally, the separation of the boundary layer is assumed 
to occur Where the Ludwieg-Tillmann skin-friction co- 
efficient extrapolates to zero. In Maskell’s nomencla- 
ture this coefficient can be written as 

U*) (dU*/dX) 


cy” = (0.246/T*) (6* (26) 


so that in compressible flow it can be computed from 


cy = (0.246/T*) (0/L) (1/M) [dM /d(x/L (26a) 


and the f 
flow 


Using Eqs. (2la), (23), (24a), (25), (26a) 


general boundary-layer 
It should be noted that any 


functions of Maskell, 
can be easily computed. 
other method can similarly be converted to the com- 
pressible flow. The Maskell method was chosen only 
for the purpose of illustration and because it appears 
to give good results in incompressible flow. 
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Supersonic Flutter of a Cylindrical Shell—IT' 


JOHN W. MILES* 


University of California, Los Angeles 


SUMMARY 


The theory of Part I (see reference 1) is applied to pressurized 
shells filled with either gas or liquid. 
negative damping ratios are calculated; the latter also are pre 


A summary of the symbols in both 


The flutter speeds and 


sented in graphical form 
Parts I and II is given below. 


SYMBOLS 


in reference | 


A = space wave amplitude, see Eq. (3.5 

A, = surface wave amplitude, see Eq. (3.5) in reference 1 

yy = see Eq. (7.16 

. = Fig. 12 

D = bending stiffness |[Eh'/12(1 — vy?) for panel 

E = Young's modulus 

F = see Eq. (7.14) and Fig. 9 

G = see Eq. (7.18) 

G, = Fig. 10 

K, = modified Bessel function of second kind (Mae 
donald’s function 

L = see Eq. (3.2) et seq. in reference 1 

M = free-stream Mach Number (J = Ll’/a 

M, = see Eqs. (7.20), (7.21), and Fig. 11 

N; = axial force per unit width of shell 

N, = circumferential force per unit length of shell 

R = radius of circular cross section 

l = air speed 

l = critical air speed (at which Vy = U’ — aseca 

Uy* = flutter speed—viz., minimum value of / 

V = axial wave speed, see Eq. (3.14) in reference | 

J = JV for free shell 

i = see Eq. (6.4) in reference 1 

a = velocity of sound in air 

ay = velocity of sound—i.e., compression waves—in 
shell material 

r = wave speed of surface wave relative to coordinate 
system (x, y, 3, ¢) fixed in still air, see Eqs 
(3.6), (3.8), (3.14) in reference 1 

f = frequency 

g = structural damping factor, such that structural 
forces are multiplied by exp (/g sgn w) for the 
harmoni¢e motion exp (/wt) 

h = shell wall thickness 

h* = minimum value of / required to prevent flutter in 
absence of prestresses or damping 

1 = imaginary unit 

m = mass per unit area of shell wall 

n = displacement normal to wavefront, see Eq. (3.3) 
in reference 1; also denotes integer 

p = normal loading per unit area on shell wall 

Pa = aerodynamic pressure resulting from wave motion 

P; = internal pressure in shell 

sgn = ssnx = + 1ifx>Oor —1ifx<0 

t = time measured in (x, y, s) coordinate system fixed 
in still air 

t = time measured in (x, i, 21) coordinate system 


fixed in shell 


Received September 20, 1956. 
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“u,v, w@ = axial, circumferential, and normal (radi 


] 


ward) displacement of element of shell 


= mapping function, see Eq. (4.5) in reference 


‘. ¥, = Cartesian coordinates fixed in still air, see | 
3.1) in reference 1 

Xi, Vi» = Cartesian coordinates fixed in (and moving wit} 
shell, see Eq. (3.1 

= argument of A,, see Eq. (A-6) in reference 

a = azimuthal angle in polar coordinates (Fig. 2 
helix angle of surface wave on shell 

8 = see Eq. (4.6) in reference | 

= see Eq $.6) in reference 1; also specific heat ratio 

for air 

6 = damping ratio, see Eq. (4.12) in reference 

6 = —éat g = 0, see Eq. (4.13 in reference | 

6 = maximum negative damping ratio 

€ = see Eq. (7.17 

Ex, € = strain parameters, see Eq. (5.7) in reference 

¢ = transverse displacement of panel (¢ = —w for 
evlindrical shell of circular cross section 

n = see Eqs. (4.3b), (4.4) in reference 1 

0 = polar angle (Fig. 2); see Eqs. (3.9a), (3.9b) in 
reference | 

Kes Ke = stress parameters, see Eq. (7.1} 

N = wavelength of space wave in x, y, 2 coordinates 

d = wavelength of surface wave in x), ¥;, 2; coordinates 

d,* = see Eq. (6.5) in reference 1 

mn = see Eq. (3.19) in reference 1 

v = Poisson's ratio 

t = see Eq. (4.3b) in reference | 

Pa = air density 

p = density of internal fluid 

p = density of shell material 

Cx; Ge = axial and circumferential prestresses 

a = see Eqs. (7.5), (7.6 

t = see Eq. (4.3a) in reference 1 

o = velocity potential (see Appendix A of reference 1) 

y¥ = see Eq. (A-8) in reference 1 


INTRODUCTION 


The following analysis is continued from reference 1 and applies 
- 


the theory of that paper to the flutter analysis of 
An important example would be a liquid-fueled rocket 


prestre ssed 


shells 
The Sections, Equations, and Figures are numbered continu 


ously from Part I (see reference 1 


(7) PRESTRESSED SHELL 


W'° NOW CONSIDER a prestressed shell, for which 
the wave speed given by Eq. (5.5) when the 
inertia of the internal fluid is neglected (a liquid-filled 
shell will be considered in the following section) may be 


expressed ast 


Au* = 2x(12(1 — yp? V4(Rh 


cf. Eqs. (6.4) and (6.5). 
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120), a) = (1/2) V*24 [(y cos a) *]? + 


( | 


\.* (A, cos @)]?; + [(o, + ¢, tan? a) pi] (7.1 


The first term on the right-hand side of Eq. (7.1) 1s 
denote the axial 
then 


) 


identical with Eq. (6.3), and o, and a, 


and circumferential stresses, respectively. It 
follows that the minimum wave speed for the pre- 


stressed shell is given by 


Mf NV V/*2 TT \Oy Pi at Ay = Ai” 
and a = 0 t Be 

The corresponding flutter speed is given by 
U* =a+ VV V*™ + (o,;/p (7.3) 


It follows from Eq. (7.3) that axial prestress in- 
creases the theoretical flutter speed, but that circum- 
The 


prestress required to prevent flutter at any given speed, 


ferential prestress has no effect on this speed. 


as determined from Eqs. (6.4) and (7.5), 1s 


7>— [(h R) V3 — v*) Jai}, 
l>a 74a) 


= p(l' — a)*{1 — (hk h*)] 7.4b) 


where /:*, the wall thickness required to prevent flutter 
The 


is, of course, invalid for values of Ll’ too 


in the absence of prestress, is given by Eq. (6.7 
result (7.4) 
close to a, and the appropriate result at the sonic point 


[(h* —h 


R| atU=a (7.5) 


is given by Eq. (6.8). If h exceeds h* the 


give the maximum compressive 


where /)* 


results (7.4) and (7.5) 
stress that could be allowed without producing flutter. 
The stresses given by Eq. (7.4b) for duralumin and 


steel, evaluating a at sea level, are 


e 13 xX 104]1 — h*)|(\ — 1)? psi (duralumin) 


(7.6a) 
= 1.3 X 10°[1 — (A h*)](MW — 1)? psi (steel 
(7.6b) 


If h is small compared with h* it is evident that the pre- 
stresses required to prevent flutter may be of the same 
order as the ultimate allowable stresses. We also re- 
mark that if the axial prestress o,* is produced by 
internal pressure the resulting circumferential stress is 
265". 

We turn now to the calculation of the negative damp- 
ing ratio 6) for speeds in excess of the flutter speed. If 
.1) is substituted in the 


the wave speed given by Eq. (7 
constraint (4.11)—viz., 1) = U' — a sec a—we obtain 
(LU — asec a)? = (1 2)V**) [Ay cos a) \,*]? + 


\i* (A, cos a) ]?; + [(o, + o, tan? a) p;| (7.7) 


Solving Eq. (7.7) for A; cos @ yields two values having 
the geometric mean ),*, of which, according to Eq. 
(4.13), 


over, if, 


the larger gives the larger value of 6); more- 


as we shall assume, 


R 7.8 


CG, >» 0; > kh 


the effects of bending are negligible at the larger wave 


length, which then is given approximately by |sub 
stituting I1™* and \,* in Eq. (7.7) from Eqs. (6.4) and 
(6.5) and dropping the bending term | 

A, = 27R sec a X 

V | i” — asec a) a,|* — o, + o,tan* a) FE] 79 


Substituting this value of A, in Eq. (4.13b), we may ex 


press 6) as a function of a—viz., 


bo = (V 3/4 pal R akh x 


\J cos a } (cosa — | x 


| ( AM cos a — 1)- — (xk, cos" a + kK sin a (.10 
where the new stress parameters 

Kr = GO, py, a, aire i.lla 

Kj = 6 pid = Q, a)€ 7.1 lb 


We remark that the velocity of 


sound is relatively insensitive to changes of altitude; 


have been introduced. 


are, therefore, relatively direct measures of 
Mach Number at 


as obtained from Eqs 


ky and x, 
the static stresses. The minimum 
which flutter can occur, say .\/ 
(7.3) and (6.4), may be expressed as 

o-RV 311 — v*)]\ "7 7.12a 


My = 1+ 4,1 + [Eh 


7.12b 


where the approximation (7.12b) is based on the as 
sumption (7.8). 

The flutter motion at the initial Mach Number (.1/ 
\J)) is necessarily axially symmetric (@ = 0), but if 
MJ > AM (for fixed «,) flutter is possible over a range of 
values of a. We may anticipate, however, that maxi- 
mum instability will correspond to a 
range of Mach Numbers, say lJ) < JJ < J, 
\J > AJ, the value of @ will increase in consequence of 
In the former range 


0 in a finite 


while if 


the constraint I) = U — asec a. 
it is expedient to write the maximum value of 6, in the 
form [setting a = Oin Eg. (7.10) } 
é)* = [(a/a,)(R/h)(p,U*/2¢;) | FM. «x 
Mes Ms MW, (7.13 
FUM, «-) = 3°72 -**,, M-4 x 
(M — 1) (WM — 1)? — 4 ¢.14) 


F is plotted in Fig. 9 for values of J between 1/, and 
Ms, where J. is the value at which F achieves its maxi- 
mum value as a function of .1/; it may be proved that 
M. > M,, so that F is of no direct interest outside of this 


range. 
The determination of 1/,; and 6)* if \J > AJ; is ex- 
pedited by rewriting (7.10) in the form 
oo = l(a a,)(R h (p,U? 20, |? | er ee 7.1loa) 
= [(a/a,)(p,U? 2p;)}? Ge a, € if o, = 
p.R/h (7.15b) 
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where 
C= Mcosa (7.16) 
e = M-*[1 — (c,/¢,)] Céch7) 
CCS &, 6) me BF tC eC — 1)-8* & 


[((C — 1)? — «,(1 — €C*)]'/* (7.18) 


The form (7.15b) exhibits the dependence of 6) on the 
ratio of dynamic (p,l* 2) to internal (p;) pressure when 
a, is a result of pressurization. The function G(C) 
vanishes for some value of C > 1, say Cy(Cy = Jp at 
a = Q), and is real and positive for all C > Cy. The 


derivative G’(C)—namely 
, i — 1)-** x 
[((C — 1)? — «(1 — eC?) ]-?/8- | x, [(9C — 4) —- 
eC?(7C — 2)] — (7C — 4)(C — 137} (7.19) 


G'(C) = 37-122-4/3 2/8307 


has, by Descartes’ rule of signs, one and only one zero 
in C > 1, say C,,, corresponding to a maximum value 


of G, say G,,. It follows that the maximum instability 


occurs at a = Oif C, > M (since then VW cos a = C,, 
does not yield a real value of a) and at a = cos”! 
(Cy/M) if Cy, < MM. The Mach Number at which 
C,, = M is the aforementioned 1/,. 


The maximum value of 6) for 17 > J, may be ex- 
pressed as 


bo* = [(a/a,)(pgU?/2p;) ]?/* Gm(xy, €), M > M, 


(7.20) 


where G,, is obtained by equating the cubic in the braces 
of (7.19) to zero, solving for C,,, and substituting the 
result in (7.18). The value of ./; is determined im- 
plicitly by C,,(x,, €), since, for fixed (¢,/¢,), 


*{1 — (¢,/¢,)]} (7.21) 


AM, — Go ; Kr, AM, 


AJ, is plotted as a function of x, for various stress ratios 
in Fig. 10. A rough approximation is found to be 


M, = 1+ (5«,/3)!/? (7.22) 


and also is plotted in Fig. 10. G,, is plotted as a func- 


tion of e for various «x, in Fig. 11. C,,, which is required 
in determining the helix angle and frequency of maxi- 
mum instability is plotted in Fig. 12. 

The graphical presentation of 6,* as an explicit func- 
tion of J for J > 1/4, would require the introduction 
of both x, and «x, as parameters, and the results would 
be too cumbersome for arbitrary stress ratios. This 
presentation is made in Fig. 13 for the important, spe- 
cial case o, = (1/2)o,, however, in order to exhibit 
more clearly the variations of 6)* with 17. 

It is very likely that panel flutter would be most 
serious for Mach Numbers in the neighborhood of J/;, 
and, in any event, the value of 6)* calculated at JJ = 
Af, furnishes an approximate upper bound to the value 
at other Mach Numbers (provided, of course, that the 
actual dynamic pressure is used in the calculation). 
Using the approximation to J1/; — 1 given by Eq. 
(7.22), the corresponding approximation to the upper 
bound of 6y* given by Eqs. (7.13) and (7.14) is 
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b9* < 0.393 (0, pia?)'/®[1 + (50,/3pia?)'/?7]-4/35 &x 
(a ‘a,)(R/h)(p,U?/2¢,) }? 7.239 


or, for a pressurized cylinder, 


5)* < 0.393 (a ay)!(o, o,)'2(p:R Eh)"® X 
[1 + 1.29(a,/a)(o,/0,)'/*(p;R/Eh)'/?]-43 x 


(p,U?/2p;)° (.23b 


The frequencies corresponding to Eqs. (7.13) and 
(7.15) may be determined by substituting A; from Eq 
(7.9) in Eq. (6.12) and expressing the result in the 


dimensionless forms 


2rfR a, = cos? a(.V cos a — 1)[(.W cos a — | — 


(x; cos? a + x, sin? a) | (.24a 


= (39/4/4)M—*x,1/4 & 
[xy /*2(C — 1)G(C; «,, €-)]-*? (7.24 


= (M — 1)[(M — 1)? — «,]-"? at 


a=0 (7.24c 


where C and G are defined by Eqs. (7.16) and (7.18 
The frequency of maximum instability may be deter 
mined by substituting «,~'°(C, — 1) and G,, (from 
Figs. 11 and 12) in Eq. (7.24b) if J > 1, and directly 
from Eq. (7.24c) if JJ < J/,;. We note that the ap- 
proximation (7.22) yields 1.58 for the right-hand side of 
Eq. (7.24c), giving a frequency 12 per cent higher than 
Eq. (6.13) for an unpressurized shell at the minimum 


flutter speed. 
(S) Ligump-FILLED PRESSURIZED CYLINDER 


We consider in this section a pressurized, liquid-filled 
cylinder—typically, a fuel tank. We shall assume that 
the density of the liquid is such that 


pidj > 2arpih (S.1 


so that the inertia of the shell (;) may be neglected 
in calculating the wave speed; this assumption usually 
would be satisfied for practical fuel tanks, since p 
would be of the same order of magnitude as (although 
smaller than) p;, while \; must be large compared with 
the wall thickness. We also shall neglect bending, the 
implicit assumption being that of Eq. (7.8); other 
things being equal, bending is less important for a 
liquid-filled than for a gas-filled tank in virtue of the 
larger wavelength at which the wave speed exhibits 
its minimum value. 
The preloads are given by 

N, = (1/2)N, = p,R/2 8.2 
and the wave speed given by Eq. (5.5) may be placed 
in the form 
(A; cos? a) 27R] + 


Vo7(A1, a) = (Eh p;:R) 


(pi pi) [(A1 cos? a) 27R]—'[1 — (1/2) cos? a 8.9 








This wave speed has the minimum value 


(Vo)min = (2QEpjh/p7R)'* atdyy = A* = 
2tRV p:R 2Eh and a = 0 (S84 
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tively, should be used to determine the damping ratio for maxi 


mum instability 


half the circumferential stress 
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Fic. 9. The damping ratio for a pressurized cylindrical shell, 7 iia: ’ 
is given by Eqs. (7.13) and (7.14), on the assumption that both , 
bending of the shell and the inertia of the internal fluid are negli- ~ . > . 
gible and that JJ < AM, (see Fig. 10); if AV MJ, the damping 4 ies z ~ 
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Fic. 10. The parameter 14), required in determining whether Fic. 13. The damping ratio for a pressurized cylindrical shell, 
Eq. (7.13) or Eq. (7.20), corresponding to Figs. 9 and 11, respec assuming that both bending of the shell wall and the inertia of the 
iiternal fluid are negligible and that the (static) axial stress 1s one 
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and the corresponding flutter speed is 
U* =a + (2Epjh/ p;7R)'"' (8.5) 


It follows that the internal pressure required to prevent 
flutter is 


bi* = p?(U — a)tR 2Eh (8.6) 


The negative damping ratio when L’ exceeds L™*, ob- 
tained by neglecting the shell inertia in Eq. (4.13b), 1s 


59 = (V/3/4)(pe/pi)?7(M cos a — 1)75/8 = (8.7) 
subject to the constraint (Il) = U’ — asec a) 


(Eh/ p;R)[(A, cos? a) /227R]| + (pi/ pi) [(A1 cos? a) + 


27R]—' [1 — (1/2) cos? a] = (U — aseca)? (8.8) 


It is evident that the maximum value of 6) given by Eq. 
(8.7) corresponds to the smallest value of cos a per- 
mitted by Eq. (8.8); this, in turn, requires the left- 
hand side of Eq. (8.8) to be a minimum with respect to 
variations of \;. Imposing the latter condition, we 
may place the results in the form 


/ ° 9 
M1 = Ai*V 2 — cos? a sec? a (S.9a) 


(MJ cos a — 1) = (MW/* — 1)(2 — cos? a)'/4cos a (S8.9b) 


where ),* is given by Eq. (8.4), and J/* is the Mach 
Number specified by Eq. (8.5). 

The exact determination of cos a from Eq. (8.9b) 
requires the solution of a sextic equation, but the ap- 
proximation 

cosa = (M — M* + 1)" (8.10) 
obtained by neglecting the variation of cos a in the fac- 
tor (2 — cos? a)!/4, is exact for MJ = \/* and exhibits 
a maximum error of 0.19(.\7* — 1)(M7 — 1.1947* + 
1.19)~! relative to unity. Substituting Eq. (8.10) in 


ICAL SCIENCES MAY, 1958 

the right-hand side of Eq. (8.9b) to determine the cor. 
responding approximation to J cos a — 1, the damping 
ratio given by Eq. (8.7) becomes 


bo* = (V3/4)(pa/ pi)?/*[(M — M* + 1)/(M* — 1) 8 y 
[2 — (MW — M* + 1)-*)"? (8.11a 


= ().43(pi/ pa)'/*(po2R/ Ep h)®!"?(M — M* + 1)58 yx 
[2 — (M — M* + 1)-7]~*/?>) (8.11b 


where the specific heat ratio has been taken as 1.4 in 
evaluating a in Eq. (8.11b)}, and pp is the free-stream 
pressure. 

It might be thought that, in analogy with the action 
of oil on disturbed water, the friction between the in- 
ternal fluid and the shell could provide appreciable 
damping. This analogy is, in fact, exact for a wave oj 
given length (A,) and speed (I’ cos a) insofar as the 
wavelength is (as already assumed) short and the oil 
film is inextensible, and we find for the (positive 


damping ratio” 
6; = (1/2)V mv;/ VA, cos a (8.12 


where vy; is the kinematic viscosity of the internal fluid. 
It does not appear that 6; could exceed 10~* for practi- 
cal configurations and a fluid having dynamical proper- 
ties approximating those of water. We remark that 
we have not proved that 6; is directly additive to the 
damping ratio calculated in the absence of internal 
friction, but in view of the negligible magnitude of 4 
the question is of little importance. 
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On Unsteady Interaction Between a Weak 
Thermal Laver and a Strong Plane 
Oblique Shock’ 


C. T. CHANG* 
The Johns Hopkins Unwersity 


SUMMARY 


A problem of unsteady interaction between a weak upstream 


temperature disturbance and a plane oblique shock is investi 


gated rhe shock is assumed to be located at the concave corner 
of a wall. The temperature disturbance is taken as a layer of 
hot air of constant strength—namely, a step function Its wave 


front is inclined at a definite angle with respect to the upstream 
unperturbed flow (see Fig. 1). It is found for such a configu 
ration that, even for weak upstream disturbance, mixed flow 


, supersonic region outside and a subsonic region inside a sonic 


circle, respectively) usually occurs downstream of the shock 
Using a general concept of conical flow, the qualitative nature 
f the flow field downstream is predicted. Subsequently, a 


Based 
on the concept of constant states, the flow in the supersonic region 
the flow 


method for its quantitative investigation is presented. 
is determined by the method of characteristics; in the 
subsonic region is determined by following a procedure similar 
to the one used by Lighthill in his investigation about the dif 


fraction of blast waves 


SYMBOLS 


se = velocity of the unperturbed main flow 
upstream and downstream of the 
shock, respectively 

1,, A = speed of sound of the unperturbed main 
flow upstream and downstream of the 
shock, respectively 

ih, M = Mach Number of the unperturbed main 
flow upstream and downstream of the 
shock, respectively 

€ = inclination of the unperturbed shock 
plane with respect to the main flow 
direction upstream of the shock 

B = inclination of the unperturbed shock 
plane with respect to the main flow 
direction downstream of the shock 

N = component of the downstream Mach 
Number normal to the shock, NV = 
M sin 8 

v*oy* = coordinate system pertaining to the 


* 


shock, with y* axis taken in the plane 
of the shock, Fig. 1 

pertaining to the 
downstream flow 


x axis taken in the direction of the main 


voy = coordinate system 


unperturbed with 


flow velocity, U’, Fig. 1 
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Ul, Us 


Subscripts 


m1 


INVESTIGATIONS SHOWED 


— 
cally and experimentally—a 


garding 


layers and shock waves. 


(1) 


problems of 


pressure of the unperturbed main flow 


temperature of the unperturbed main 
flow 
specific heat of the medium at constant 


pressure 
ratio of specific heats of the medium 


perturbation of pressure, temperature, 


and velocity of the flow field, respec 


tively 


dimensionless quantities of the above 


perturbations as defined in Section 


Eq. (2.2 
rotational part and rotational part of 
the perturbed velocity field i, respec 


tively 


components of velocity perturbation 


normal to and parallel to the unper 
turbed shock plane, respectively 

inclination of the upstream temperature 
jump with respect to the main flow 
velocity U,, Fig. 1 

drift speed of the upstream temperature 
jump along the shock, Fig. 1 

shock displacement 

local shock angle and local shock propa 
gation speed, respectively 

conical coordinates corresponding to the 
physical coordinates x, y, and +, Eq 
(2.4 

derived conical coordinates appearing in 
the text 

polar coordinates corresponding to &’, n’ 

polar coordinates corresponding to &, n 

for the flow 


in the hyperbolic region, Eq 


a characteristic coordinate 

(3.2 

a radius vector, defined by Tschaplygin’s 
transformation, Eq. (3.11 


refers to the unperturbed main flow 

refers to upstream flow parameters in the 
vicinity of the shock 

refers to downstream flow parameters in 
the vicinity of the shock 


INTRODUCTION 


both theoreti- 


great interest re- 


interaction between thermal 


Most of these investiga- 


tions, however, are confined to interactions steady in 


time. 


Two types of interaction were considered; the 


upstream temperature disturbance was taken either as 


weak! * or strong.® 


For weak interaction, the per- 


turbed pressure field downstream of the shock, depend- 
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PLANE OF UNPERTURBED 
SHOCK — 





HOT AIR 





me 
* 
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Fic. 1. Model of interaction. 


ing on the shock angle concerned, is either entirely 


supersonic! or entirely subsonic.* For strong inter- 
action, mixed flow near the foot of the shock was re- 
ported.’ The object of this paper is to show that if 
the flow is unsteady, then even for weak upstream dis- 


turbance, mixed flow, usually, occurs downstream of 


the shock. 

As a model of interaction, the shock is taken to be 
originated at the concave corner of a wall, and the 
upstream disturbance is taken to be a jump disconti- 
nuity of constant strength, its discontinuity front in- 
clined at a definite angle to the upstream main flow 
velocity. The upstream disturbance thus drifts at a 
constant speed along the shock. From dimensional 
reasoning, the perturbed flow in the downstream must 
be of the conical nature. Based on the general concept 
of conical flow,’ the qualitative nature of the down- 
stream flow field can be predicted at ease. It is 
shown that the flow field downstream of the shock is 
divided into two regions; a supersonic region outside 
a sonic circle and a subsonic region inside a sonic circle, 

Finally, a method for the quantitative evaluation of 
the perturbed flow is presented. In the supersonic 
region, the solution is obtained by using the concept 
of constant states and the method of characteristics; 
in the subsonic region, the solution is obtained by fol- 
lowing a general procedure similar to the one used by 
Lighthill in his investigation concerning diffraction of 


blast waves.°® 


(2) FORMULATION OF THE PROBLEM 


The Model 


As a model of our analysis, the medium is taken to be 
an ideal gas, inviscid and obeying perfect gas laws with 
constant specific heats; the shock is assumed to be 
originated at a concave corner of a wall. In view of 
most cases occurred in practice, the downstream main 
flow is assumed to be supersonic. The upstream tem- 
perature disturbance is assumed to have a constant 
jump, A7’; the normal to its discontinuity front is in- 
clined at an angle 6 with respect to the upstream main 
flow Ul). The front of the discontinuity therefore 


advances at a constant speed 
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C, = [cos 6 cos (6 — e)]U, on 
along the shock (see Fig. 1 


Boundary Conditions 

For convenience of analysis, the flow parameters 
will be replaced by their corresponding nondimension 
alized ones. Thus, if Ap, A7, and Anu denote th per- 
turbations of pressure, temperature, and velocity 
their corresponding dimensionless parameters will be 


denoted by 


bp = Ap 'Ypbm T= ATT,, u = An'A (22 
In the above, p,, and 7), are the pressure and temper- 
ature of the unperturbed main flow, y is the ratio of 
specific heats of the gas, and A is the speed of sound of 
the medium. We shall also adopt the convention that 
whenever no subscript is attached, the flow parameters 
are those in the region downstream of the shock; for 
similar flow parameters upstream of the shock, a sub- 
script 1 will be used—for example, the pressure of the 
main flow upstream of the shock will be denoted by 
Pim, While the pressure of the main flow downstream of 
the shock will be denoted by ?,,. 

Referring to a set of coordinate axes ox*y* with the 
y* axis taken in the plane of the unperturbed shock, 
(Fig. 1), using the conservation laws of mass, momen- 
tum and energy, the flow conditions across the shock 


can be written as! 


r+ Ay 11. cos 2B 711 

p+ a Avy T 4 1 .\ cos ——s TT} } 
“.* Asi ' | a3. cos B Wy" 4 131 Ms 
v+ + | 10 | TT 4) | 0 


In the above, u* and v* are components of the per- 

turbed velocity taken normal to and tangential to the 
unperturbed shock plane respectively. Ay, As, 
T1, 11, ... etc. are parameters depending on the shock 
strength, defined as the ratio of main flow pressures 
across the shock—i.e., p,»,/ pin. Their explicit form 
was given in reference 6. 

Besides these shock conditions, the velocity com- 
ponent, v, normal to the wall and the shock displace- 
ment, ¥, at the corner of the wall must vanish at all 


time. 


Governing Equations and General Nature of the Solution 


As the effect of the upstream disturbance on the 
shock depends only on the drift speed C,, from dimen- 
sional reasoning, one would expect that the perturbed 
flow parameters downstream of the shock must also be 
functions of some physical parameters having the 
dimension of a velocity. Specifically, if one follows 
the idea of Busemann,‘ and introduces a set of coordi- 


nates defined by 
§E=x/7r, n= y/t (2.4 


then the physical parameters, e.g., p, T, u, v, etc. of the 


pe 
flo 
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2.1) | perturbed flow field must be functions of (£, 7). The 
flow field in the (&, 7) plane can be viewed as a repre- Ps 


















: sentation of the perturbed flow at the instant of 7 = 1. 
The flow field at any other instant can be obtained by 
a proper magnification or reduction in proportion to the 
meters instant to be considered—i.e., by taking a proper cut 
‘sion along the 7 axis (see Fig. 2). With respect to the pres- 
le per- ent frame of reference, the conservation laws of mass 
locity, and momentum can be written as & - a | 
=e mass: (M — &)h) — p, = —(%& + 2,) |} 
momentum: (M — &)u; — nu, = —pe> (2.5 
29) ) (M — &)x, — nv, = —?2D, \% 
mnper- For analytical works, it is more convenient to refer —_ 
tio of them to a set of polar coordinates defined by — 
ind of Te oe. fn /(t — ie ' 
ete (é — VW)? + 77], 6 = tan—'[n (E - ad Ade = Cs 7 
leters - we 
5+ for In terms of the polar coordinates, Eq. (2.4) becomes 
sub- rp cos #u, + sin @v, — 
1 the | (sin 6 r) ug + (cos 8 r)vs 
d by rl cos 8 p, — sin 6 pg r (2.7 
m of sin 0 p, + cos 0 py r ~—s a 
1 the It is easy to see that the particular solution of the s 
1ock. velocity field associated with the pressure field 1s irro- 
nen- tational, or it is the drift velocity associated with a "ag 
hock sound wave; while the homogeneous part of the solu- 
tion represents the rotational part of the velocity field aay cate ere 
: a a Fic. 3. Flow field in hyperbolic region 
i.e., the shear waves. From Eq. (2.7), three character- 
istic directions can be obtained; they are 
y | C,: 6 = constant whole region downstream of the shock, the last two 
eS Oe ae 28 characteristics exist only outside the sonic circle i.e., 
; forr > 1 only. The characteristics C+ and C— can be 
a C_: 6 — sec.~'r = constant ee é a 
2.3 shown to be simply lines tangent to the sonic circle; the 
Der- While the characteristic 6 = constant is valid for the circle itself isan envelope of the characteristics, hence, it 
the is sometimes known also as the Mach circle. Physically, 
) these characteristics represent discontinuity fronts of 
wn y sound waves generated downstream. Similarly, the 
e } characteristic C, corresponds to accompanying dis- 
- continuity fronts of entropy (or temperature disturb 
ance) and vorticity waves formed at the shock. 
on With the aid of these characteristic directions, a 
am { qualitative examination of the downstream flow field 
- can be made at once. A specific example is given in 
Fig. 3. Thus, if 0o represents the drift speed C, of the 
upstream temperature discontinuity, then 01 and Oe 
™ . will represent the sound wave and the vorticity wave 
(also the entropy wave) generated, respectively, at the 
he shock. 12 will represent the sound wave reflected at 
al the wall; 23, its reflection at the shock; 2e, the accom- 
ed panying vorticity and entropy wave formed at the 
be edge of shock. 
he incoming “Soon The sound field in regions (0), (1), amd (2) are hyper 
vs disturbance - bolic in nature—i.e., the pressure disturbances are 
i- waves with constant amplitudes. The sound field in 
region (3) is of a mixed type, hyperbolic outside and 
1) elliptic inside the sonic circle. If one wants to explore 


T=At further about the downstream flow field quantitatively, 


le Fic. 2. Conical coordinates of flow field. one finds that it is most convenient to analyze the 
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pressure field first. From Eq. (2.7), after the elimin- 

ation of « and v, a single governing equation for the 
perturbed pressure / results, thus, 

r*(r° — 1)D,, + r(2r* — 1)p, = Pon (2.9) 

Since Eq. (2.9) is hyperbolic or elliptic according to 

whether + is greater or less than unity, this confirms 


our previous statement that the downstream flow field 
is of the mixed type. 


(3) METHOD OF SOLUTION 


Flow Field in the Hyperbolic Region 
When r > 1, Eq. (2.9) can be reduced to a simple wave 
equation with respect to the characteristic coordinates 


éand o thus: 


Poo ai Pee = 0 for r > l (3.1) 
where ao = sec.—! (7) (3.2) 


Eq. (3.2) can be viewed as a transformation which 
maps a point (7, @) in the physical plane to a point 
(o, #) in an image plane. In particular, the Mach 
circle is mapped into the coordinate axis o = 0; the 
characteristics which are not fixed in the (7, 6) plane are 
mapped into two sets of lines inclined at 45° to the 
coordinate axis in the (o, @) plane (see Fig. 4). From 
Fig. 4, it is evident that, unless the shock degenerates 
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into a Mach wave, the family of characteristics coming 
from the wall will intersect with the shock; thus reflec. 
tion phenomena will continue (although the strength 
of the reflected wave might be zero in some special 
case),! until the last characteristic hits the coordinate 
axis ¢ = (—+1.e., the Mach circle. Clearly, the zone 
of influence for a fixed field point P in the hyperbolic 
region is the region on the left of the two characteristics 
passing through the point. Although the dominating 
effect at the point is coming from waves generated at 
the shock caused by the upstream disturbance and the 
first reflection at the wall, if the point under consider. 
ation is not too far away from the corner of the wall 
then clearly the effect of the corner will have a great 
significance. Consequently, in order to investigate 
the flow field in the hyperbolic region, first of all one 
needs the knowledge of the flow field near the corner 
Once the wave front of the upstream disturbance has 
passed the corner of the wall, the transient effect of the 
disturbance is gone. As a result, one might expect 
that the flow field near the corner would be the same as 
that produced by a steady disturbance; its region of 
validity is limited to that defined by the characteristics 

i.e., up to the line 01 for the sound field and up to 
the line Oe for the vorticity field as indicated in Fig. 3 
Using the fact that we are looking for a solution near 
the corner of the wall where the vertical component of 
the perturbed velocity, v = 0, the boundary conditions 
at the shock, Eq. (2.5), in the present case, reduces to 


P 4 = Av T; T+ 7 (.\J cos B Vs 
us sin B = As, Ty + m3, (VW cos B)y,» 3.3) 
us cos B = Tuy, 


where #4 is the component of the velocity perturbation 
along the wall immediately behind the shock. Since 
the region upstream of the shock consists of constant 
states’ separated by the discontinuity front of the given 
temperature jump, the hyperbolic region outside the 
Mach circle in the flow field downstream of the shock 
must also consist of regions of constant states separated 
by characteristics, and it follows that the above three 
equations will be just sufficient for the determination 
of the three unknown quantities p, u, and y,*. In par 
ticular, the sound field is given by 

(AoC — Az). cos? 8 + As; B sin B r 


Po = ‘ : 1 3.4) 
\J cos? 8 + B sin B 


» 


where c = F31/ F21, B — — 14 
and (to = —po/ M 


Since the shock displacement in terms of the conical 


coordinates 1s 
y(y*, r) = x(n*)r Dut 
it follows then at the instant 7 = 1, orn* = C, A, 


Xo = [Agi 7) cos B/ (a sin B — m3 J cos? B)| X 
(C./A) 3.6) 


In a similar manner, the vorticity and temperature 
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field near the corner can be determined. To explore 
further the flow field in the adjacent zones, it might be 
worthwhile to take another close look at the governing 
equation for the perturbed pressure in its original form, 
Eq. 2.9). If one compares our present problem, a 
two-dimensional unsteady flow with no restriction on 
rotationality, with the known case of three-dimensional 
irrotational steady conical flow about supersonic 
wings,’’ one notices that our governing equation for 
the perturbed pressure (more general, for any param- 
eters p, u; or v; of the sound field) has the identical 
structure as the governing equation for the velocity 
perturbation in the latter case, provided we identify 
our present pressure perturbation p with the axial 
component —w cot (sin~! | AJ) in the latter case. 
Using this analogy and the results of reference 8, it is 
easy to verify that the general solution of the flow 
field for r > 1 is given by 


%A+<a 
g(w) dw 7 | f(w) dw 
e 
%A—o 
u; = cos w [g(w)] dw + 
rA+e 

| cos w [f(w)] dw 
. 


2A—«@ 


sin w [g(w)| dw + (3. 


sin w [f(w)] dw 


%a 
u., = COS Ww [T(w Cw 


= sin w [J(w)] dw 
where [f(w), g(w), and /(w) are periodic functions of pe- 


riod 27 such that 


a» > 
PS a 
2\W) dw COS WwW [£\W ad@ 

0 0 


| sin w[g(w)] dw = 0 (3.8 
0 


e 


Similar relations hold for f(w) and /(w). In the above, 


u;, v;) is the irrotational part and (u,, v,) is the rota- 
tional part of the perturbed velocity field, respectively. 

Denoting the pressure jump across the characteristic 
01 by Ap and the velocity jump for the shear wave 
across the characteristic Ve by Ae according to Eq. (3.7) 
the jump of the flow parameters across these character- 


istics are given by 


Ap = pi — po 

Au; = (ui)) — (u = Ap cos 4, 

Ao, = (%) — (ede = Ap sin B, (3.9) 
Au, = (uy) — (Us) = Ae cos ay 

Az = (v,)r7 — (v.), = Ae sin Q) 


The flow parameters immediately behind the shock 
iow are seen to be expressible in terms of the two un- 
The local shock displace- 
Using 


known jumps, Ap and Ae. 


ment x(n*) furnishes an additional unknown. 
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the last three equations of the shock conditions of Eq. 
(2.3), with the upstream disturbance 7°_, deleted, these 
unknown quantities can be determined algebraically. 

In view of the fact that discontinuities of flow param 
eters can occur only across characteristics and the con- 
cept regarding constant states,’ the entire flow field in 
the hyperbolic region can be evaluated in a similar 
manner. 

In summary, the physical phenomena occurring in 
the hyperbolic region, generally, will consist of the fol- 
lowing: generation of waves of all three modes 
sound, vorticity, and temperature spottiness—at the 
shock: reflection of sound waves at the wall: and re- 
flection of sound waves at the shock with simultaneous 
production of waves of the other two modes. The 
process will continue or stop at any one of these stages 
whenever the last characteristic of the sound wave hits 
the Mach circle. An example illustrating all the above- 
mentioned phenomena has been given already in Fig. 3. 


Flow Field in the Elliptic Region 

Inside the Mach circle, 7 < 1, the flow field is elliptic. 
The governing equation for the perturbed pressure, 
Eq. (2.9), 
terms of a new set of polar coordinates (p, @) thus, 


can be reduced to a Laplace equation in 


(1p) (0 Op) [p(Op p)] + 
(1 p*) (02p 062) = 0 forp<1 (3.10) 


where the radius vector p is given by the Tschaplygin's 


transformation?® 
p=[l—-V1l-r]r (3.11 


The appropriate boundary condition now is to be 
specified along the boundaries of the region—namely, 
the shock condition on the segment ab, the “‘upstream”’ 
disturbance on the are bb,c obtained from known infor- 
mation in the hyperbolic region, the wall condition on 
the segment ced, and the ‘downstream’ disturbance on 
the arc da (see Fig. 5 Among these conditions, the 
one at the shock is the most complicated, for instead of 
dealing with one single flow parameter, one has to deal 
with a system of simultaneous equations stating the 
relationship between the flow parameters and the local 
shock displacement. After the elimination of the 
shock displacement x, one obtains from the shock 


condition the following system of simultaneous equa- 


tions: 


Cp, — i, = 0, Bp, — ww, = 0 (3.12 


where # and @ are components of the perturbed velocity 
along the coordinate axes E and %, respectively. With 
the aid of the governing equations themselves a single 
equation for the perturbed pressure now can be derived. 
With reference to Fig. 3, confining our attention to 
values of the drift speed, C, A below that of Ob, the 
boundary conditions for the pressure perturbation can 
be stated as the following. 
(1) On the segment ab, 
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pe (C — N)q + (BN)/9 
Ps 1— N° 
for —= —N at 
: ; P «(S.13) 
4 = — [McosB — (C,/A)]f 
(2) On the are bb;, p = pz (3.14) 
(3) On the are bic, p = pe (3.15) 
(4) On the segment ced, p,, = 0 (3.16) 
(5) On the areda, p= 0 (3.17) 


In addition, one has the compatibility conditions 
that the perturbed velocity @ along the shock at 7 and 
at 7. must conform to those values obtained previously 
in the hyperbolic region, or from Eq. (3.12) one has 


V1 — N? “na 
(B/n)p, da = | 0,d7 = —; (3.18) 
=-V/1 — N? at 


where @; is the component along the shock of the ve- 
locity perurbation in zone (3). 

According to Eq. (3.11), the domain of interest 
dabb,ced in the (7, 6) plane is mapped into a curvi- 
linear quadrangle d’a’b’b’\c’e’d’ in the (p, 6) plane. 
The segment ab corresponding to the shock is mapped 
into a circular are of a circle having radius equal to 
V(1/N2) — 1 and center at (—sin B/N, cos B/N). 
The Mach circle r = 1 is mapped into the unit circle 
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p = 1 with points on the circumference mapped into | 
themselves. The unit circle p = 1 intersects the circle 


corresponding to the shock, orthogonally at the image 
points a’ and b’ of the points a and b (see Fig. 5). To 
proceed further analytically, it would be much more 
convenient to deform such a complicated domain into 
a somewhat simpler one. By introducing a complex 
variable ¢ = uw + iv = pe”, it can be shown that the 
region inside the curvilinear quadrangle in the complex 
¢ plane can be mapped into the upper half plane, y, > 0 
in the complex 2; = x; + 7); plane through the following 


transformation! 


ee 3 K (1 — ae’®?) — (1 — ae’®) ce* | 
Sn* <2¢— in am Saree es 3.19 
lr (1 + ae”) — (1 + ae”) fe” f 
where 
@ = x — cos! (N) 


a = Vecos’? £8 — sin*® ¢ cos (¢ + 8B) 


| 
K = | dt (1 — #) (1 — k)]'? 


0 
k = \(zInG/z) = 
modulus of elliptic integral 


(45.20 


V cos 8 + sind + Vecos B — sin @ 
Vcos 8 + sin @ — Vcos B — sin ¢ 


Eq. (3.19) maps the boundary of interest in the ¢ 
plane into the real axis in the 2; plane. The are d’a’ is 
mapped into the negative half of the real axis; the arc 
a’b’ corresponding to the shock is mapped into the 
segment 0 < x; < 1; the are b’b’\c’ is mapped into the 
segment 1 <x, < k 
sponding to the wall is mapped into the half line x, > 
k~? (Fig. 5). For a fixed shock angle 8, the modulus 


> and the segment c’e’d’ corre- 


k can be shown as a function of the shock strength, 
Pm Pim. ‘Result of calculation shows that for 6 = 
30°, and for values of shock strength p,, 1, above 2.0, 
up to the accuracy of four figures behind the decimal 
the value of & is very close to unity.! 

Following a method used by Lighthill,°® the present 
problem can be formulated conveniently in terms of a 
function F(z,) defined by 


F(z) = py, + tpz (3.21 
The function F(z,) is regular in the upper half plane y; > 


0. In accordance with the boundary conditions stated 
previously in the (7, 6) plane, it is also subjected to the 


following conditions along the real axis y,; = 0. 

(1) When x; < 0, F(z;) is real (3.22 
(2) When 0 < x, < 1, 
: V tan? ¢ cot? 6 — 1 
arg [F(z,)] = tan— J = — 
( C; f 

{V tan? ¢ cot? é — 1I SS aie 
tan™* : (3.23 

{ c f 
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Vtan? @ cot?é—1= 
! a) tan | (m 2A )Sn 1% ¥1,R { » 4 
(3.2 


> o { é , . l ( 
1 (a?) tan? } (7/2K)Sn—! (Ym, R)} 


the constants C, and C. are to be determined from 


C+0=1/D, GAC, = (CD) -1 (3.25 
Cand D, in turn, are given by 
( 13) 7 D = —(m4/ 7) [NO — N?)] (3.26 
3) When 1 <x11<k 
F(z,) is real (3:37) 
1) When x, > R 
F(z,) is pure imaginary (3.28) 


In addition, the function F(z,) is to be supplemented 
by the following conditions. 

(5) Near the point 2; = x», the image point of b,’, at 
which a pressure jump occurs in the (7, 6) plane, must 


behave as 


F(z) ~ (ps — p2)/ (21 — Xo) (3.29) 
6) The compatibility condition, Eq. (3.18), further 
requires that 
B "11 + (1 a?) tan? } (4 2AK)Sn-! (Vix, R)f 
, ; Xx 
sn @Jo 1 — (1/a?) tan? }(r 2K)Sn—! (Vx, R)f 
[Im F(z) ],,=0 dx, (3.30) 


With some amount of manipulation, it can be shown 
that the first four conditions are satisfied by a function! 





HT(2,) = filti)fo(ai) (R72? — 24) 
where 
ae 7 
, (1 a) tan } (4 2K)Sn-! (V2, )f 
C, + 27 x 
L 1 — (1 a?) tan? )(@ 2K)Sn-! (V2, h)} 
‘ . 4 ‘ 
F (1/a) tan }(@ 2A)Sn-! (V 21, R)j 
L 1 — (1a?) tan? }(@ 2A)En1(V a, R)§ 
(3.31 
and 
j j i ee tan | —- (3, = 4) (s, + 7)] 
exp - : : x 
| a P pia tan x ee (s; — 7)/(s + 72)] 
YUN) 
g 
dxiy 
1+ Vy 
where g(x,) = arg [fi(t1)] on x, > 1 


In view of the remaining conditions, our required 


function, as can be verified, is given by 
F(z = D,|D2(z, — ia 1] F7(z,) (23 — Xo) (3.32 


The two constants, D,; and D2, are to be determined 
irom the remaining conditions stipulated by (5) and 


(6), respectively. 
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t) CONCLUSION 


By considering the shock as originated at the concave 
corner of a wall, and the discontinuity front of an up- 
stream temperature jump as inclined at an angle to the 
upstream main flow velocity, a problem of unsteady 
interaction between a thermal layer and a plane oblique 
shock is studied. It is found that under such a con 


dition, even for weak upstream disturbance, mixed 


flow generally occurs in the region downstream of the 
shock 
and an elliptic region inside the sonic circle. 


namely, a hyperbolic region outside a sonic circle 
Based 
on the general concept of conical flow, the downstream 
perturbed flow caused by the interaction is examined 
qualitatively at first. Subsequently, a method for 
quantitative investigation of the downstream flow field 
is formulated. Analytically, this is accomplished by 
formulating a single governing equation for the pressure 
perturbation alone. This equation, then, is shown to 
be reducible to a simple wave equation for the hyper 
bolic region, and to a Laplace equation for the elliptic 
region. 

For the flow in the hyperbolic region, the solution is 
started by assuming that the region near the foot of the 
shock is the same as that corresponding to a steady 
flow; using the concept of constant states, solutions for 
adjacent regions were obtained by the method of char 
acteristics. For the flow in the elliptic region, the 
solution is effected by using the method of conformal 
mappings after a Tschaplygin’s transformation. The 
domain of interest is finally mapped into a half plane. 
Due to the topological equivalence of the region to a 
rectangle, the mapping function involves an elliptic 
integral. The modulus occurring in the elliptic func 
tion is shown to be related to the shock strength and is 
found to be very close to unity for fairly strong shocks. 
Finally, following a procedure of Lighthill, the solution 


for the pressure perturbation is obtained. 
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Experimental Investigation of the 
Aerodynamics of a Wing in a Slipstream’ 


MARC E. BRENCKMANN* 
Institute of Aerophysics, University of Toronto 


SUMMARY 


An experimental study of a wing in a propeller slipstream was 


made in order to determine the spanwise distribution of the lift 


increase due to slipstream at different angles of attack of the wing 
The 


and at different free stream to slipstream velocity ratios 
results were intended in part as an evaluation basis for different 
theoretical treatments of this problem 

The comparative span loading curves, together with support- 
ing evidence, showed that a substantial part of the lift increment 
produced by the slipstream was due to a ‘‘destalling’’ or bound- 
ary-layer-control effect. The integrated remaining lift incre- 
ment, after subtracting this destalling lift, was found to agree 
well with a potential flow theory 

An empirical evaluation of the destalling effects was made for 


the specific configuration of the experiment 


SYMBOLS 


Ss = area of wing immersed in slipstream 
R = propeller radius 
( = wing chord 
y = spanwise coordinate, origin at midspan 
rg = total wing lift 
AL = total lift increase due to slipstream 
AL»y = direct lift increase due to slipstream 
C; = wing lift coefficient in uniform flow (nondimen- 
sionalized with free-stream g¢ 
ACy dest. = additional wing lift coefficient due to indirect 
slipstream influence (nondimensionalized with 
free-stream go) 
CI = local lift coefficient (nondimensionalized with 
free-stream go) 
Cy = local lift coefficient (nondimensionalized with 
“average” 7) 
= local lift per unit span 
= dynamic pressure in free stream 
lj = local dynamic pressure in slipstream 
qj = equivalent dynamic pressure of slipstream 
Aq = equivalent dynamic pressure increment of slip 
stream 
q = |(qj + qo)/2], ‘“average’’ dynamic pressure 
= slipstream thrust 
o = (qj — do)/(Gj + Q) Slipstream strength parameter 
\ = empirical factor 
a = angle of attack of the wing 
8 = angle of rotation in slipstream 
y = ViGj/9 


INTRODUCTION 


| anand in STOL aircraft configurations has put new 
emphasis on propeller slipstream deflection as a 

Received September 30, 1957. 

+ This research was carried out under the supervision of Dr 
G. K. Korbacher and Dr. H. S. Ribner, with the support of The 
DeHavilland Aircraft of Canada Limited, while the author was 
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* Now at The DeHavilland Aircraft of Canada Limited. 


Past 
as 


Iron 


means of creating lift at low forward velocities 
experiments on this subject are not systemati 
this new point of view, mainly because they were con 
ducted for other purposes. In fact, only two papers 
are usually cited when experimental values are com 
the 


are known to be inadequate for 


pared with theoretical results. On other hand 
the earlier theories! ° 
practical configurations and the validity of recent 


theories” * + is not definitely established. 

The investigation reported here was intended to 
supply experimental data in the form of the spanwise 
lift distribution of a straight wing in a propeller slip- 
stream for various angles of attack and for various 
free stream to slipstream velocity ratios. The present 
paper is an abridged version of UTIA Technical Note 
No. 11.7 


EXPERIMENT 


The measurements were conducted in the UTIA 32 
by 4S in. low-speed wind tunnel, using a free-stream 
speed variation from zero to approximately 100 ft. sec 

A wing of constant chord of S in. and of constant 
profile (NACA64A418) was designed to span the tunnel 
over 32 in. from wall to wall. It was composed of 19 
narrow sections independently springmounted on 
Each spring was straingaged so as t 
Aero- 


common spar. 
allow measurement of the lift at each section. 
dynamic continuity was obtained by means of a thin 
elastic skin wrapped around the wing, sealing the gaps 
between the sections but transmitting very little force 
from one section to another. The wallmountings oi 
the spar formed pivots for the variation of angle oi 


attack at the 1 4 chord line. 





1. Downstream view of the apparatus in the wind tunnel 
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The slipstream was produced by an S-in. diameter 
propeller driven at 15,000 r.p.m. by a 1 2-hp. pneu- 
matic motor which was suspended independently up- 
stream of the wing in a pusher configuration. 

The variation of the free stream to slipstream ve- 
locity ratio was achieved by varying the tunnel speed 
and keeping the air pressure for the pneumatic motor 
constant. The experimental layout is shown in Fig. 1 
and the construction of the wing in Fig. 2. 

Test runs were made for five tunnel speeds at three 
angles of attack with and without slipstream. Fig. 5 
shows the slipstream cross sections in dynamic pressure 
and helix angle for the different runs. 


REDUCTION OF RESULTS 


The curves of Fig. 3 showed that the slipstream 
velocity cross section was far from the ideal case of a 
step increment above freestream velocity. An equiva- 


lent step increment was therefore defined as: 
R 

Aq = (1 /R*) [ 2 y (qj — qo) dy 
/7 0 


In order to reduce the lift distribution to nondimen- 
sional form, a reference quantity was chosen which does 
not disappear for either of the extreme cases of uniform 
flow or of static jet; this was termed the “average” 
dynamic pressure g defined as 


(go + qj) 


7 with g; = go + Aq 


g = 


Thus, when g; > go (uniform flow), g — go, and when 
qo > O (static jet), gq > g; 2. 
The nondimensional quantities used in the presenta- 


tion of the spanwise lift distribution are then 


local lift coefficient cr = 1G. 


slipstream strength parameter 
o = (g; — qo)/ (g + go) = Aq 2g 


This particular parameter o arises naturally in the 
theoretical formulation of reference 4+. The results of 
slender body theory applied to this problem are also ex- 
pressed most simply as a function of o. The more 
common coefficient c, = //qec is easily calculated 


from those defined above as 
oo = 4/ — @) 
The total lift increase due to the slipstream was ob- 


tained as 


AL = 


. 
— ) 
| { by ith slipstream lwithout slipstream / dy 
¢/ span 


The slipstream thrust was obtained as 


> 


4dr i} (g — V 4,90) y dy 


Total lift increase values were rendered nondimen- 


T= 


sional in the forms 


AL/aT, AL/2maSAq, AL,/aT, AL,/2raSAq 
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as explained in the discussion of the results. 


RESULTS AND DISCUSSION 


(a) Lift Distribution 


Figs. 4(a) and 4(b) illustrate the spanwise distriby 


tion of lift as ¢, versus y R curves for different ¢. It 


can be seen that: 


(1) the rotation of the slipstream causes the lift in 
crease to be shifted towards the side of the wing which 
is in the propeller upwash. This influence diminishes 
with higher wing incidence and seems compensative 
in overall effect. 

(2) the limits of the direct slipstream influence are 
sharply defined and do not vary with a@ or o. The 
present configuration did not cause any appreciable 
spread of the slipstream. 

(3) the lift distribution outside of the slipstream 
shows a marked increase relative to the corresponding 
case without slipstream, especially for high wing inci- 
dence. This increase extends up to close to the tunnel 
walls, where the tunnel boundary layer renders the read- 
ings uncertain. This phenomenon suggests some flow 
modification over the whole wing which delays the stall; 
the evidence is developed in the following section. 


(b) Slipstream Destalling Effect 


Fig. 5(a) illustrates the lift characteristic of the wing 
at its midspan section in uniform flow at the various 
Reynolds Numbers encountered in the free stream dur 
ing the test runs. Fig. 5(b) shows the lift characteris- 
tic of the wing at a section outside of the slipstream in 
flow with and without slipstream. Fig. 6 shows the 
chordwise shifting of the center of pressure with angle 
of attack at sections inside and outside of the slipstream 
boundary for uniform flow and for flow with. slip- 
stream. 

These curves show the improved stalling behavior 
of the wing in the flow with slipstream. The progres- 
sive flow separation at the trailing edge which is nor- 
mally taking place at angles of attack above 10° is de- 
layed up to angles above 15°, the most probable cause 
being boundary-layer suction exerted both by the low- 
pressure region in the slipstream above the wing and by 
the jet pump action of the mixing region in the slip- 


stream boundary. 


(c) Total Lift Gain Due to Slipstream 

In order to evaluate the overall slipstream effect, the 
total lift increase AL is presented: 
This 


quantity can be termed slipstream turning effectiveness 


(1) as AL aT versus o (dotted lines in Fig. 7 


of the wing. The value unity signifies that AL has arisen 
from a turning of the thrust vector through an angle 
equal to the angle of attack. It can be seen, however, 
that curves for different a do not coincide, and AL is 
nonlinear with a and o. 

(2) as AL/2raSAq versus o (dotted lines in Fig. $). 
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Fic. 9. Graphical estimation of the destalling effect 


Here 27aSAg represents the lift increase according to 
The 


general shape of the curves is the same as for those 


simple strip theory without induction effects. 


under (1) above, the nonlinearity with a being some- 
what less pronounced. 

Comparison of both representations with a closed 
form solution given by slender body theory (adapted 
from reference 2) shows good agreement only for o near 


unity (static jet). 


(d) Direct Lift Gain Due to Slipstream 


In order to separate the lift due to boundary-layer 
modification from the total lift increase, the quantity 


AL, is introduced* as 


. 
AL, _ | (doen slipstream 
J span 


Thus, AL, is an approximation 


beithout sipstream ) dy 


‘“‘destalled flow 


as illustrated in Fig. 9. 
of the lift increase as it should be predicted by potential 
For comparison AL, is represented in the 
It is seen that 


flow theory. 
same forms as AL in Figs. 7 and 8. 
(1) AL, aT and AL, 2raSAq are independent of the 
angle of attack a, and (2) both sets of values agree very 


well with the solutions given by slender body theory. 


(e) Indirect Lift Gain Due to Slipstream 


Having correlated the direct influence of the slip- 
stream with potential flow theory, an attempt was 
made to formulate the indirect influence empirically. 
The c,; 
to be modified by the slipstream to yield higher ¢;, a; at 
distortion 


versus @ characteristic of the wing was known 


higher a. Therefore, \ was defined as a 


factor of the characteristic as 
[C L “+ (AC I )destalting | C i 
S | witout slipstream dy / Sf bwichout slipstream dy 


destalled flow’’ 


A= 


This factor \ plotted versus a for different o¢ is shown 


in Fig. 10. The relation 


* Dr. H. S. Ribner suggested the idea of subtracting the de- 
stalled lift and correlating AL, with slender body theory. 
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is shown to be a close approximation of the experimental 


results for the configuration investigated, indicating 


that d is independent of ¢ for ¢ > 0.25. 


CONCLUSIONS 


The slipstream velocity profile, rotation, and defor- 
mation are concluded to be of secondary importance. 

The total added lift due to the passage of a slipstream 
over a wing appears to consist of two parts: (a) a 
direct ‘“‘potential flow lift’ that is a function of slip- 
stream strength, and (b) an indirect ‘“‘destalling lift 
of comparable magnitude that is sensibly independent 
of slipstream strength above some threshold value. 

For the present configuration it was possible to corre- 
late the potential flow lift with an existing slender body 
theory and to approximate the destalling lift by an em- 
pirical function of angle of attack. 

The controlled use of the slipstream velocity to ener- 
gize the boundary layer would appear to be attractive. 
Since destalling effects will depend on the wing con- 
figuration, its stalling behavior, the Reynolds Number, 
etc., the indirect lift gain obtainable in STOL applica- 
tions will have to be studied for optimization on the 


specific configuration considered. 


REFERENCES 


! Durand, W. F., (Ed.), Aerodynamic Theory, Vol. IV, Koning 
C., “Influence of the Propeller on Other Parts of the Airplane 
Structure,’’ Julius Springer, Berlin, 1935 

2 Graham, E. W., Lagerstrom, P. A., Licher, R. M., and Beane, 
B. J., A Preliminary Theoretical Investigation of the Effects of Pro 
peller Slipstream on Wing Lift, Report No. SM-14991, Douglas 
Aircraft Co., Santa Monica Division, November, 1953 

3 Rethorst, S., Lift on a Wing in a Propeller Slipstream as Re 
lated to Low-Speed Flight, Aeronautical Engineering Review, Vol 
15, No. 10, pp. 42-48, October, 1956 

4 Ribner, H. S., Theory of Wings in Slipstreams, Unpublished 
Report, DeHavilland Aircraft of Canada Ltd., March, 1957 

5 Smelt, R., and Davies, H., Estimation of Increase in Lift Du 
to Slipstream, British ARC, R&M 1788, 1937. 

6 Stuper, J., Eznfluss des Schraubenstrahls 
Leitwerk, Luftfahrtforschung, Vol. 15, No. 4, 
as NACA TM 874. 

7 Brenckmann, M. E., Experimental Investigation of the Aer 
dynamics of a Wing in a Slipstream, UTIA Tech. Note 11, April, 


1957 


und 


Fliige 


translated 


aus 
1938; 


} 


S 


S 


) 


— 


[deg] 


mental 
icating 


defor- 
nce. 

stream 
(a) a 
f slip- 
x lift’ 
ndent 


corre- 
be dy 
n em- 


ener- 
ctive. 

con- 
nber, 
plica- 
1 the 


ming 
plane 


eane, 
Pro 


uglas 


s Re- 
Vol 


ished 
Due 


und 


ated 


lero 
pril, 


Readers’ Gorm 


RIEF REPORTS of investigations in the aeronautical sciences and discussions of papers pub- 


lished in the JOURNAL are presented in this special department 


Entries must be restricted 


to a maximum of 1,300 words or the equivalent of one JOURNAL page including formulas and 


illustrations 


Publication is completed as soon as possible after receipt of the material 


l. The Editorial 


Committee does not hold itseif responsible for the opinions expressed by the correspondents 





Contents 


A Vectorial Representation of Aerodynamic Forces Acting 
on a Thin Rectangular Wing Oscillating Harmonically in 
Supersonic Potential Flow. . JAGANNATH P. CHAWLA 


Note on Ballistic Trajectories and Orbits 
RICHARD F. HUGHES 


Flow Against a Vertical Plate With Small Suction 
R. Foore 


Hydromagnetic Effects on Stagnation-Point Heat Transfer 
JosEpH L. NEURINGER AND WILLIAM McILRoy 


Magnetohydrodynamic Analysis of Heat Transfer Near a 
Stagnation Point VERNON J. Rossow 


Space-Time Correlations of the Fluctuating Wall Pressure 
in a Turbulent Boundary Layer. . W. W. WILLMARTH 


Blunt Body Separation at Supersonic Speeds 
Pau F. BRINICH 
Some Solutions of the Navier-Stokes Equations With Time 
Dependent Density 
ROGER D. SULLIVAN AND COLEMAN DUP. DONALDSON 
On the Character of the Instability of the Laminar Boundary 
Layer Near the Nose of a Blunt Body A. M. KUETHE 


Relative Importance of Free-Stream Vorticity and Self- 
Induced Pressure Gradient on a Flat-Plate Boundary 
Layer HAROLD MIRELS 


One-Dimensional Transient Heat Flow in a Multilayer Slab 
W. F. CAMPBELL 


Asymmetric Starting for Hypersonic Wind Tunnels 
ROBERT H. JOHNSON 


Note on ‘‘Stability Equations for Conical Shells’’. . 
PAUL SEIDE 


Application of the Witte Rearranging Method to a Typical 
Structural Matrix BERTRAM KLEIN 


The Conservation Equations for Multicomponent Gas Mix- 
tures in Arbitrary Coordinate Systems. .F. A. WILLIAMS 


6) 
ord 


330 


w 





329 


A Vectorial Representation of Aerodynamic 
Forces Acting on a Thin Rectangular Wing 
Oscillating Harmonically in Supersonic 
Potential Flow* 


Jaganneth P. Chawla 

Aerodynamics Department, Systems Development Laboratories 
Hughes Aircraft Company, Culver City, Calif 

January 2, 1958 


A 


METHOD of representing aerodynamic lift and moment 


acting on a thin rectangular wing oscillating harmonically 


in supersonic potential flow is indicated. This may help the 
practicing flutter engineer gain a better understanding of the 
nature of aerodynamic forces exciting the structure. This rep 


resentation should also prove useful for (a) solution of flutter 


determinant by vector methods, (b) flutter studies by analog 


methods, (c) correlation with lift and moment of oscillating air 
foils measured in a wind tunnel, (d) dynamic testing of structures, 
and (e) flight testing for flutter. 

The section lift and moment acting on the airfoil may be ex 


pressed as follows: 


L = 2pcV?2 k? (Li + tLe) [h/(c/2)] + (Ls + iLy)a} l 
M = —pc?V? k? (M+ iM) [h/(c/2)] + (Mg + iMyal (2 
Defining a, = h/V = ikh/(c/2), the sectional coefficients become 
cq. = L/(1/2)pV2c = 4k(Le — ili ja, + 4RAL; + iLy)a (3 

Cm = M/(1/2)pV2c? = 
—2k( Me — iM,)o, — 2k2(.M; + iM,)a } 


* The author is grateful to the Hughes Aircraft Company for permission 
to publish this note, and should like to express his appreciation to his col 


league, Dr. William M. Bleakney, for heipful discussions 
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With (0c;/Oa) = 4/8, we can write 
-0.05 fo) 0.05 0.10 
Cte, /( Cia) = kB(L. — iL, (5 0.05 7 
Ctal (Cte = k?3(L; + 1L;4) (6) 
Cag lic = —(1/2) ka( VM tM, (7 
0 
Cme/(Ciqimo = —(1/2) RB (Ms + iM, (8 : j . 
The right-hand sides of these equations may be written as aa 
Ae'®. Hence, the amplitude and phase of the lift vector at any / = 
frequency are known in terms of their values at zero frequency; -0.05 + + or bd 
also, the magnitude and phase of the moment vector are known - ¢ 
in relation to the lift vector at a prescribed frequency of oscilla / . / ~ 
tion and Mach Number 25 A 
Plots of Eqs. (5) to (8), based on the results of linear theory ou 
(see reference 1) and for the axis of pitch at mid-chord, are 5 
shown in Figs. 1 to 4 
Figs. | and 2 show that the lift vector lies in the right half ell 
plane for all values of k and JJ. Its amplitude decreases mono 
tonically with k, for values of & of practical interest in flutter in 
work. The deviations in phase and amplitude from the steady — 
state value decrease as Mach Number is increased 
Fig. 3 shows that for mid-chord location of the axis of pitch, ee 
the moment vector due to vertical translation leads the vertical “0.20 ) 
velocity 4 by 90° for small values of & and lies in the upper half 
plane for values of k or j}@ = [2Rk\?/( 10? 1]! of practical 
interest 
4 a 
f 











Fic. 1. (0¢;/Oa),)/( Oc; / Oa) = kA(L tL. 














-a 
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Fic. 4 (OCm/ Oa) /( Or Ow ) = —(1/2)k*8( M/ T V 


Fig. 4 shows that for low supersonic Mach Numbers the m 


ment vector leads the angular displacement vector a by 90° for 


small values of k. As Mach Number is increased, lead angl 
decreases changing to a lag angle. The limiting value of the lag 
angle for JJ? > 1 is 90 
that for 1/*? > 1 amplitude of the moment vector increases 


Another noteworthy observation i 
linearly with frequency of oscillation. Position of the moment 
vector in the upper half plane, that is, with a lead angle of 0 
180°, is indicative of instability in pitch 

Curves of a similar nature can be drawn for other locations of 
the axis of pitch. For rectangular wings of finite aspect rati 
similar plots may be prepared with values of lift and moment 
integrated over the span, including the effect of pertinent mode 


shapes considered in the flutter analysis 
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— PATH of a particle moving in the earth’s gravitational 


field (see Fig. 1) is an ellipse which may be described as 


follows: 


Y 1— V1 — (2 — y)y sin? a sin 6 





YO™ for 
| angle 
the lag 
tion ts 
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rel where 
R ry /Ro 
tude of particle motion with respect to radial dire« 
it location r = 
dinates of path 
ticle velocity at location 7 , oO 
earth's radius 
vitational acceleration at the earth’s surface 


ristics of the particle motion are readily discernible 


pection of the form of relation (1 rhe maximum 
m distances of the path from earth’s center are given 
id —1 values of sin 6, respectively 
x sin? a@ 
1+ Vi 2 yy sin? « 
Relation (2) indicates that the conditions necessary for the 
t move in an orbital path about the earth are (a) y < 2, 


R 2Rv/ro, the escape velocity condition, and (b 


Condition (b) may be restated as a lower limit for the value of 


by inserting Ry for y minimum in relation (2) and solving for 


hus the inequality which must be satisfied for an earth 
tellite 1s 
2A(re/R ] 
2 9 : 3 
ry? sin? a/R)? l 
vniere sii? a g R r 


[The restriction on @ is required to satisfy condition (b) for 
values of y less than 2. The minimum particle velocity which 
ll produce an orbital path is given by @ equal 7/2 


Vy?/gR = 2/(ro/Ro + 1) (ro/R (4 


This orbit would measure a minimum distance of Ry and a 
from the earth’s center. If the path is 
R 2 1), the 


maximum distance of 7 
everywhere relatively near the earth’s surface (7 
minimum required orbital velocity becomes the familiar V gR 

Relation (1) may also be used to examine the ballistic trajectory 
f a particle moving between two points on the earth’s surface 
rhe trajectory range, distance traveled measured along the earth’s 


surface, may be determined from 


range = 2R,(@ @ rmar — 8 








Fic. 1. Orbital path of a particle moving in the earth’s gravita- 
tional field 
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here Ry, has been set equal to 7 é@r from relation (1) is 
7/2. @ is the coordinate angle for ry = i 
range = 2R)( 2/2 6 = 2R, cos sin 4 
i 1d 
l si : i 
range = 2R,) cos : 
Vi 2 1 ' 
For a fixed value of y, the value of a whic na izes range is 
a = sil ] V2 t 
max. range 2Ry cos 2v I 2 7 
id 
r Ry @ max. range = y/(2 + | Vi + 8 
For 0 ¥ 1, optimum a varies between 7/4 and 7/2 and 


maximum range varies between 0 and zR 
The maximum altitude attained by a particle directed straight 


outward from the earth may be obtained from relation (2) by 


allowing a —~ 0 


tmar/Ro (for a = 0 1/(1 7/2 a 
REFERENC!I 
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Flow Against a Vertical Plate with Small 
Suction* 


J. R. Foote 

Professor, Division of Engineering Science 
Purdue University, Lafayette, Ind 

January 14, 1958 


I’ IS WELL KNOWN that the complete Navier-Stokes equations 
for boundary-type flow can be solved for the case of a stagna 
tion point in plane steady flow. Here that problem is contained 
as a special case of flow having wall suction, typified by flow 
through porous cloth. The simplest effect of porosity is to alter 
a boundary condition without affecting the differential equation 

In the notation of reference 1 the equation and boundary con 


ditions are 


¢ + oy ¢ + | =) / | 
0) = w>0, ¢ (0) =0, ¢ = ]\| 
The independent variable is £, the wall is & = 0, and the suction 


parameter is ¢ This mathematical problem is related to the 


physical problem through 
“= Ag(t), v= Aaye(é 
t=ax, a= Va/», A= Voa 


ax, v = ay, which give tl 


At large distances, u = e meaning of 


the parametera. At the wall, 


u= -Mm =U, go = U/A 
and g can be quite large since A ~ pv This fact greatly affects 
the possibility of quickly truncating the series form of the solu 
tion 

A series method, set forth in detail by Fettis,? is applied, which 
leads to a system of linear differential equations with forcing 


terms dependent on previously calculated functions. Put 


* This note was written under the sponsorship of the National Science 


Foundation, Contract NSF-G3050, New York University Reproduction 
in whole or in part is permitted for any purpose of the United States Govern 


ment 


The author wishes to acknowledge the of Israel 


valuable assistance 
Krongold in programmer training and the UNIVAC staff at New York Uni 


versity in regard to machine time and computing aids 
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ot) = & + pelt) + pelt) +... (2) 


where p is a parameter to be suitably selected. The first term 
trivially satisfies the boundary condition ¢g’( 2) = 1. The re- 
maining procedure is to satisfy the wall conditions with ¢(£), 
while the higher-order terms satisfy zero conditions. These 
terms act to bring the series into agreement with the solution of 


the differential equation. In the present case it happens that p = 
1 may be chosen and no other parameter, in particular 1/¢», 
arises in a natural way to form either a power series or an asymp- 
totic series. 

It then is sufficient to use as boundary conditions 


)=0, gi+1(0) = 0, / 


e(0) = —g(0O) = 1/p,) gi+1(0) = ON 
Substituting the series into Eq. (1) and collecting terms by powers 


of p, for each 7 


; ae Pa ae = , 
¢i + §y;/’ — 2o;’ = file, o ..., Gi-1) = fil) (4) 
where 
fh(é) = 0, folb) = gr? — gi¢i 


and for each 7, fi(€) can be written explicitly as a homogeneous 





g je] ¢ 


r | 








G 
~ 
o 














ra 
/ 


o | 2 "We 
5 = 


'T KK 





























°o 


Solutions for ¢ at three values for ¢ 
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function. All the equations have the same form of comp} 


mentary function: 


3/92 , . £2/4 f 
A; + B,[(&/3) + &] + Cre */*D_4(é) 
where A;, B,;, C; are determined by boundary conditions, 
D_(é) is a parabolic cylinder function The method of yar 
ation of parameters gives as a particular solution, for 7 > 2, 
ot) + joolE (€3/3) + €] + iex()e*°/* D_a(t f 
where 
no: 
— 2 . } 
a(é) = —1/3 ¢ gi(n)dn, 
i ’ 
ol) = gi(n)D (n)dn, 7 
0 
9 
o3(é -{ e™ F (n? + 1)gi(n)dn, 
( 
gi(n) = e?/4 fin 
When the boundary conditions (3) are applied to Eq. (2), it js 
found that each constant A;, B;, C; contains the factor 1/p 
because the fi(£) are homogeneous in the functions ¢)(£ Thus 
it is sufficient to take p = 1, since neither the equations nor th 
boundary conditions place a necessary requirement on p. 
For 7 = 1 the particular solution is zero and boundary condi 
tions give a final form 
e(E) = [go — (2/8) V2/r] + 2V/2/r ee IDE) (8) | 
For 1 > 2, 
14, = —(2/3) V2/xBi, 
° 
B; = -{ gi(n)D-x(n)dn, C; = } 4 
Numerical calculations of the function { 
gE) => Et+ oiléE + o(E 
s/=—s $ s 
and its first derivative were made for g = 0, 1, 2, 3, 10, and 10/ 


using Simpson’s rule on the integrals, with an interval of 7/+/2 = 


0.005. 


The first three cases are shown in Figs. 1 and 2. For 
g) = 10, results seemed to be good to two digits, and for g = 


100 to only one digit. But for usual dynamic conditions of flow 


through cloth g = 250 seems appropriate. Hence, additional 
terms in the series, or a direct numerical solution of the differ 


ential equation, will be required for such a problem 
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Hydromagnetic Effects on Stagnation-Point 


Heat Transfer 


Joseph L. Neuringer and William Mcllroy 
Scientific Research Staff, Republic Aviation Corporation 
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December 9, 1957 


INTRODUCTION 


IT A PREVIOUS PAPER! by the authors, an investigation was 
initiated to study the 
electrically conducting viscous fluid with an applied magnet 


hydromagnetic interaction of au 


field in stagnation-point flow. In the present note (see als 

reference 2) this investigation is carried a step further, and the 

corresponding influence on aerodynamic heating is presented 
The notation employed will be that already clearly outlined ™ 


reference 1. However, a few more definitions will be necessary 
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READER 
mperature of the gas behind the detached nose shock 
of the surface at the stagnation point, assumed 


mperature 


ery much smaller than 7,* 


is velocity leaving the shock front 


eed of the body—i.e., the gas speed entering the shock 
/ ratio of the specific heats behind the shock 
} Prandt! Number of the gas behind the shock, assumed equal to 


kin-friction coefficient 27/pob?*x* 
Re , Reynolds Number, based on the local value of x* and the local 
external tangential flow velocity bx* 
Nusselt Number gx*/k(T-* T w* 
DISCUSSION 
[The theoretical model studied here is the same as that in 


reference 1. By virtue of the incompressibility assumption 
ind temperature-independent fluid properties, the momentum 
solved independently of the 


distribution, 


ind continuity equations can be 


ergy equation, and the solution for the velocity 


can then be used to solve for the temperature in the energy 


equation 
It is further possible to show that the solution reduces to the 


familiar solution obtained by Squire (see page 631 
subject to certain conditions 


reference 3, 
for incompressible stagnation flow, 
which will now be enumerated 
Making the usual boundary-layer assumptions in the energy 
3, page 612), but including the Joulean 
heat-flow 


equation (see reference 
dissipation term J**/o, the convective equation can 


be written in the form 


V*(o7*/oy*)| = 
+ p(OU*/dy*)? + (J*?/a 


pClo| U*(OT"/On*) + 


U*(Op*/Ox*) + k(027*/dy* 
This equation may be considerably simplified. In the first 
place, the small perturbations in the electromagnetic field 
strengths due to the magneto-hydrodynamic interaction may 


that the Joulean heating 


when the 


follows, therefore, 
as ol’ *2B,*?, 
following the 


be neglected It 


may be written simply definition of 


in reference 1 is used. Further, standard pro 


cedures of order analysis, it can be shown, (see reference 3 
page 613), that both viscous heating and Joulean heating are 
likely to be small compared to conduction heating. In other 


words, the heat-convection equation can be rewritten in the 


following simpler form: 











} U*(OT*/dx*) + V*(OT*/Oy*)} = «(0?T*/dy*? 
where «x = }k/c,po{, provided that 
{U,*2/c(T.* — T 
and }[h2(0 Uy*?/c,(T.* — 7 : 
TABLE 1 
{= 10° 

0 u(o t cyRey- 
to) MO, bg “ 6, (oF %, “Ke | Vo (oof so be t Nu | 
| 
© | 0.000/1.000 1.000 |1.000/1.000 1.000 |1.000/|1.000 4.96 | 





-100 | 0.100) 0.964/0,977 |1.005 0.982 0920 1.033 0952 | 5.09 
-300 0.900/0.726 | 0.852 | 1.037 |0.883 0.620 1.260 0.782 6.02 
-400 | 1.600 0.622 |0.789 1.055 0.834 0.492 | 1.430 | 0.703 6.72 
~500 | 2.500 | 0.535 |0.732 |1.077 |0.788 0.392 | 1.624 0.637 7.50 | 
~700 4.900) 0.412 |0.642 | 1.114 |0.715 0.267 | 2.052 | 0.548) 9.14 
— + = a 2 ee i = SC ——— 
B ee 2x10° 
‘ae a a | A a ae me ——n | 
~ 300 0.4850 0.830 0.910 |1.02! | 0.929 |0.760/ 1.138 | 0.865 | 5.53 


500 | 1.250 | 0.667/ 0.815 | 1.048 | 10,994]0.960/ 1.348 10.728 6.38 _| 
| 


‘ 
= ne [ome —_ T —_ 
800 |0.640 | 0.782 | 0.884| |.028 oon [1.19 [oeee | [5.74 
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Shear-stress and heat-transfer ratios versus interaction 
parameter h? (0)/4 


Fic. 1 


are both small compared to unity. For strong shocks the first 


condition is always satisfied, and for the values of the interaction 


parameter considered in this note, 47(0)/+7, the second require 
ment is also satisfied 

the above equation for temperature, use is made of 
The 
assumed 
of x* 


final form 


To solve 
the solutions for U* and V* obtained in reference 1 
effect of the Joulean heating on the 
negligible, so that a 


be used as boundary condition 


inviscid flow is 
independent 
The 


namely, 


constant temperature 


may at infinity 


of the result is identical to Squire’s solution 


Nu ca =a Pr Re tai 


2 , ) 
l/a(Pr) = J J exp _ Pr f u(y)dy~ dy 
o ¢ 0 ails, 


, Was calculated on IBM for various values 


where 


This integral, 1/a(Pr 
of h2(0)/>5 
The ratio of the heat transfers, with and without the magnetic 


field, is 

q/qo = asa b/by)''* 
Here } and by are the external, inviscid, velocity-gradient param 
eters with and without a magnetic field, respectively In 
reference 1 the ratio 6/b) has been shown to equal 
$1 + [h2(0)/y]} 2 


RESULTS OF THE ANALYSIS 


From the above expression for the heat-transfer ratio q/qo, 


it can be seen that it consists of two parts; which 


one part, a/ ao, 


is a function of the dimensionless velocity profile u(y) alone, and 


a second part which depends on the character of the external 


inviscid flow field. As a result of the analysis (see Table 1 

it is readily seen that the first part increases slightly with in 
creasing (#2(0)/7), while the second part decreases \ similar 
situation exists with respect to the shear-stress ratio. In ref- 


erence 1 it was shown that 


t/r) = )u"(O)/uo"(O)} }b/dof 
ie., a product of a nondimensional velocity-profile-dependent 


* ‘ ” ” . . . 
portion, }“"(O)/uo"(O)}, and an external inviscid-flow-dependent 
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Fic. 2. Typical velocity profiles. 


portion, (b/b))*?. In Fig. 1, the ratios g/qo and 7/7» are plotted, 


and it is shown that the heat transfer decreases less rapidly 
than the skin friction 

The effect of the magnetic field on the dimensionless velocity 
profiles, tangential and normal to the plate, is apparent from 
Fig. 2. The profile approaches more and more the shape com- 
monly associated with a turbulent boundary layer, and _ this 
would lead one to believe that the magnetic field has the added 
effect of increasing the stability of the flow. 

The value of the Reynolds analogy parameter, 


crRe,x/Nu = 2u"(0)/a(Pr) 


is tabulated in Table 1. It is seen that this is not a constant 
with varying h?(0) 

Finally, it must be emphasized that the present solution exists 
only for small values of distance x*, and for values of h?(0)/¥ 
consistent with the restrictions noted in the discussion. For 
sufficiently high values of this interaction parameter, the Joulean 
heating term assumes sufficient magnitude to be included in the 
energy equation. 


CONCLUSIONS 


As a result of the above analysis, it is concluded that it is 
possible to reduce the heat transfer at a stagnation point (when 
the gas there is ionized) by application of magnetic fields normal 
to the surface. The reduction can be of considerable magnitude, 


even in the presence of magnetic fields of reasonable strength 


and gases of relatively small electrical conductivity. As a 
typical example, from Fig. 2 it appears that for y = 10° and 
h?(0)/y = 4.9, the heat transfer is reduced 28 per cent with a 


corresponding decrease of 45 per cent in the skin friction 

With the typical electrical conductivity, density, external- 
velocity parameter ), etc., used in reference 1, this would 
require a magnetic field strength of 3,200 gauss 
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Magnetohydrodynamic Analysis of Heat 
Transfer Near a Stagnation Point 


Vernon J. Rossow 

Aeronautical Research Scientist, Ames Aeronautical Laborato; 
NACA, Moffett Field, Calif 

December 20, 1957 


SUMMARY 


Phe change in the skin friction brought about by a magnetic f er | 
dicular to a flat plate was studied theoretically in reference | for 
flow near a two-dimensional stagnation point In the analy t ¢ 
cally conducting fluid was assumed to be incompressible and t ive 
stant fluid properties The change in the imposed magnetic fic n 
by the electric currents in the fluid was also considered 

The purposes of this note are to present the results of a heat-trar 
inalysis for the problem of reference 1, and to compare the skin fri 
found in reference 1 with that obtained by the more approximate meth 
reference 2. The intent of the second objective is to find the degree of | 
proximation made by ignoring the secondary magnetic field which induce 


by the motion of the fluid through the basic imposed magnetic field Eq 
( 
i 


ANALYSIS 


Bye DIFFERENTIAL EQUATIONS Which describe the motion 
the fluid and the transport of thermal energy in the presenc¢ liff 


of the magnetic field are, from reference 2, 
u(Ou/Ox) + v(Ou/Oy¥) + (aBo?/p)u = 
(1/p)(Op/Ox) + v|(O07u/Ox?) + (O07u/04 val 
u(ou/Ox) + v(dv/Oy) = (1/p)(Op/dy) 
v|(077/Ox?) + (07 /O¥ lb ten 
(Ou/Ox) + (Odv/Oy) = V0 2 


C,u(0T/ox) + C,v(0T/dy) (oBy?/p)u? = 





(C,v/Pr)(0?T/dy?) + v(du/dy)? + (u/p)(Op/ox 3 


where the symbols a and By denote the electrical conductivity 
of the fluid and the imposed magnetic field strength. Eqs. (la 
and (1b) are the Navier-Stokes relations with the additional tern 
in Eq. (la) accounting for the magnetic field-fluid interaction 

A set of ordinary differential equations for the velocity and 
temperature distribution in the boundary layer near a stagnation 
point may be obtained by introducing the new variables, 


n= Va/vy fa 


u = ayxe'(n) th 


v= —V vaog(n) ( + 


> 


p Pp =(e 2 )ao?|x? 4+ (2Va v/ay?)d(n 4d 


T(n) = Ta — 


(U2 + V2)/2C,]e(n) — (Te — Tw0(n) (4 


where u and v are the velocity components in the boundary layer —_ 
and l’ = ax and V = ay are the velocity components outside 

the boundary layer. Here x is measured along the wall fron 
the stagnation point and y is normal to the surface into the fluid 
The stagnation and wall temperatures are, respectively, 7%; and 
7. When the magnetic parameter m is zero, UU) = aox and V = 
ayy. Since the boundary layer is assumed thin, the velocity 


its edge is 


Therefore Eq. (4e) becomes 
T(n) = Tse — (a?x?/2C,)0i(9) — (Tse T\)62(n Hf 
When the above variables are introduced into Eqs. (1), (2), and 
(3) and the Prandtl Number is taken equal to one, the following 
set of ordinary differential equations is found: 
¢’'’ + ge” — (¢’)? — me’ +1 =0 


(0,” + gf)’ — 2¢'8,)(a?/ao?) 


where m = oBy?/pa The boundary conditions are \ 
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TABLE | 
: Reference | 
aja = a°o 
mt ¢”(0 g’( = do” 4," TT q/Q h2(0)/- "(0 
0) 1. 2326 1. 0000 1.0000 0.5705 1.0000 1. 0000 0 1. 2326 
0 0953 1. 1844 0.9535 0.9091 0.5601 0. 9609 0. 9827 0.1 1. 1863 
() 3381 1.0756 0.8452 0.7160 0.5350 0.8726 0.9377 0.4 1.0735 
0. 6529 0. 9600 0, 5264 0. 5053 0.7789 0. S857 0.9 0.9611 
() 9923 0.8615 0.3847 0.4774 0.6989 0.8369 1.6 0.8635 
,= yo’ 0 
+ a / in 
6, = @ = () at y = (0 é 
6. = |] \ 
g’ = a/a) i 
at edge of boundary layer ae 
6 =0 |) Ie _— q 
= Qo T ee q 
Eq. (7) is the differential equation found by Squire and presented ~~ i 
~ 
in reference 3, but the boundary conditions are different The ‘ _ 
solution described by Eq. (6) with m set equal to zero or a con i ms 
stant not been considered previously Eq. (5) is the usual at 
lifferential equation found for the velocity in the boundary layer os r 
ir. a two-dimensional stagnation point with an additional term l be 
rising from the magnetic field fluid interaction, % : a 
Eqs. (5), (6), and (7) were integrated numerically for several m= CA 
- < <a . Po 
values of the magnetic parameter m on the IBM 650 Magnetic : 
Fic. 2. Reduction in viscous shear and heat transfer as a fune 


Drum Data-Processing Machine using the Runge-Kutta method 
The essential results are presented in Table 1. The velocity and 
temperature profiles in the boundary layer are shown in Fig. 1 
forx = VY v/ag and apv/2C, = Tr 1 
files would more nearly resemble those of the flat-plate boundary 
The reduction in the 


The temperature pro 


ver if a larger value of x had been chosen 
skin friction and heat transfer brought about by the magnetic 
field is Shown in Fig. 2, where the viscous shear and heat transfer 


re defined as 


Ou 


a 
= pyr ¢ (0 
=a | 0 \ v 


-+=) k(T Ty; “ 4.'(0 
- OV J/y=0 N: 


rhe ratios of these parameters in the magnetic case to the non 


magnetic case (subscript zero) are then, 


T/T) = |¢ (O ¢ (O)im=0 


0.'(0) 


m=0 


rhe values of the magnetic parameter were chosen so that a 
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FIG Velox ity and temperature profiles for sey eral values 
olm at x = V v/a Tt — Te = aqv/2C a) Velocity. 
b) Temperature. 


with 


tion of the magnetic parameter 


direct comparison could be made with the skin-friction results of 
It is to be Table 


parameter of one, an error of less than 1 part in 400 is caused by 


reference 1 noted in 1 that, up to a magnet 
ignoring the induced secondary magnetic field in the boundary 
layer solution. Hence, the simplification introduced in reference 
2 causes no appreciable error and is a good approximation for 
many problems 

The heat transfer is noted to decrease with increase in magnetic 
field strength independent of the distance from the stagnation 
point. Experimental results may depart from the present theory 


when density variations (which were neglected) arise either be 
cause of large temperature variations or because the results are 


applied too far from the stagnation point 
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ger PROBLEM of turbulent boundary-layer induced noise in 
side an airplane fuselage has been discussed and formulated 
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pressure fluctuations which excite the airplane skin make a sta- 
tistical description of the problem necessary. Thus, one asks 
for the statistical properties of the wall-pressure fluctuations in 
the turbulent boundary layer in order that the statistical proper- 
ties of the motion of the airplane skin and the internal sound field 
generated by this motion may be computed. 

It has been conjectured’: ? that the governing phenomenon for 
the generation of internal noise is a type of resonance in which 
traveling waves in the skin are excited by the pressure fluctua- 
tions in the turbulent boundary layer. The magnitude and other 
characteristics of this resonance are determined by the proper- 
ties of both the skin and the pressure fluctuations. Phillips‘ 
has recently discussed an analogous problem for the generation of 
ocean waves by turbulent wind. He also finds that a resonance 
phenomena dominates the situation and shows how convected 
normal pressure fluctuations can cause appreciable wave ampli 
tudes. 

In this note measurements of the streamwise space-time cor- 


relation 
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p(x, t — +) P(x + 4a,t 
are described. p is the pressure, ¢ is the time, x and Aw are dis 
tances measured in the stream direction, and 7 is the delay time 
The measurements were performed in a low-noise-level tunne} 
using barium-titanate pressure transducers® and Skinner's elec 


tronic time correlator.’ The time correlator samples one signg 


and stores it in a digital memory device for a time 7, whereupon 
a sample from the other signal is multiplied with the stored 
sample. The output of the device is averaged, giving the time 
correlation of the two signals : 

In Figs. la, 1b, and le the streamwise space-time correlation 
of the wall pressure in a turbulent boundary layer is presented 
for three different free-stream velocities. At each stream ye. 
locity the pressure transducers were separated by 0.104 ft., 0.25 
ft., and 0.50 ft 
delay 7, called the optimum time delay, at which the correlation 


If one forms the ratio 


For each separation distance there is a time 


between the pressures is a maximum 
for the nine cases in- 


Ax/ropt. it is found that Av/top:. = 0.85U 
vestigated. The boundary-layer displacement thickness was 
6* = 0.0062 ft. for all the measurements and the Reynolds 


Number was approximately 9 X 10° at the lowest speed ar 
15 X 10® at the highest speed 

One can interpret the peaks in the correlation of the pressures 
as the result of a convected pattern of turbulent flow in the 
boundary laver. It should be noted that a given pattern of tur 
bulence has effectively been destroyed, 


Pi(x, £ — r)pe(x + AX, 7)/V pr?V po? 0.1 
ina time 7 = 1006*/L 


A complete report of these measurements and of measurements 


of the spectrum of the wall pressure is being prepared 
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Blunt Body Separation at Supersonic Speeds 


Paul F. Brinich 
Lewis Flight Propulsion Laboratory, NACA, Cleveland, Ohio 
January 17, 1958 


Se TUBE EXPERIMENTS conducted at the Cornell Aero- 
nautical Laboratory showed the existence of a separated 
flow region ahead of a blunt body. This phenomenon was ob- 
served at a Mach Number of 5.0 and was found to be caused by 
the impingement and subsequent upstream reflection of dust 
particles present in the shock tube gas 

Blunt body experiments conducted at the NACA Lewis lab: 
ratory in the 1 by 1 ft. variable Reynolds Number wind tunnel 
This 


revealed the same phenomenon at a Mach Number of 3.1 
spark 


both with instantaneous 
with high-speed schlieren motion 


phenomenon was observed 
photography and 


Schlieren photographs of a blunt cone with and with- 


schlieren 
pictures 
out the separated conical flow region are shown in Figs. 1 and 2 
Out of approximately 20,000 frames taken with the high-speed 
motion picture camera, flow separation was observed only once 
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] Separated flow conical shock configuration Free 
stream Reynolds Number, 6.25 & 10°/in 





Conventional blunt tip shock configuration. Free- 
stream Reynolds Number, 6.25 & 10°/in 


On this occasion it appeared on three consecutive frames at a 
film speed of 1400 frames per second. The description of the 
separation phenomenon given in reference 1 agrees with that 
observed in the NACA motion pictures, inasmuch as the sepa- 
rated region suddenly erupted as a dust particle was reflected up 
stream at the blunt nose and gradually collapsed as the particle 
was again blown downstream. This relatively infrequent occur- 
rence of separation (three out of 20,000 frames) in the 1 by 1 ft 
variable Reynolds Number wind tunnel indicates that the phe 
nomenon probably has an insignificantly small effect on the heat 
transfer and boundary-layer measurements made in this facility 

It was suggested in reference 1 that the separation region at 


_the blunt nose might be used as a means of reducing the heat 


transfer in the stagnation region if a sufficient number of dust 
particles existed in the airstream. Whereas this method might 
work at the nose, Figs. 1 and 2 suggest that the heat transfer 
would be increased over the downstream portions of the body. 
\ comparison of the transition locations indicate that the 
extent of laminar flow is greatly diminished when the nose sepa- 
ration occurs, and that turbulent heat-transfer rates would be 
encountered over much of the afterbody 

A burst of turbulence which accompanied the unsteady flow 
separation is also indicated in Fig. 1 and seems to encircle the 
model. The presence of this burst suggests that early transition 
may be the result of increased spot production frequency and 
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} 


perhaps that the spot production in this case is somehow ri 


lated to the presence of dust particles in the airstream 
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Some Solutions of the Navier-Stokes Equations 
With Time Dependent Density 


Roger D. Sullivan and Coleman duP. Donaldson 
Aeronautical Research Associates of Princeton, Princeton, N.J 
January 20, 1958 


IT THE COURSE of an investigation of time-dependent boundary 
laver problems an apparently new family of exact solutions 
of the two-dimensional Navier-Stokes equations has been found 
They apply to cases in which the density of the fluid is a function 
of the time only, and the other variables are functions of one 
space coordinate and the time The viscosity is taken as 
constant 

These solutions represent the boundary layer on an infinite 
flat plate moving with variable velocity at low Mach Numbers 
in an environment of changing density, or a generalized Couette 


flow produced by parallel walls which move in and out as well as 


in their own planes 

Taking the density, p, as a function of the time, ¢, and the ve 
locity components in the x and y directions, u and 2, respectively, 
as well as the pressure, p, as functions of y and ¢, the Navier 


Stokes equations with constant viscosities, w and A, are 
p|(Ou/ Ot) + v(Ou/ Oy = p(O*u/Oy" l 


p|(Or/Ot) + v(Ov/O% = —(Op/Ov) + 


and the equation of continuity ts 
p’ + p(dv/dOv) = 0 Ss) 


where the prime denotes differentiation with respect to 
Eq. (3) can be integrated immediately to give 

} 

if we restrict ourselves to cases where y = Urepresents a W ill so 

that 7(0, 4) = 0. If there is another flat wall at y = A(t), from 

conservation of mass we must have 


ph = constant 5 
If the relative vertical motion of the plates is to be regarded as 


density variation, this motion must be slow 
enough so that the density is essentially a function of ¢ only In 


written v = (h’/h)y and in particular, 


the cause of the 


this case Eq. (4) can be 
v = h’ when y = h, as is to be expected 

Substitution of v from Eq. (4) into Eq. (2) and integrating gives 

p= 1/2)pl(p”/p) — 2 p’ p)*jve + pl, t 6 

The main interest lies in the determination of wu. Eq. (1), with 


7 determined by Eq. (4), becomes 
p(Ou/Ot) = pw(O*u/Ov*) + p'y Ou/O% d 


We note that a solution of this equation is given by 


“u = Cp) 
Considering the flow between two infinite plates, one at 5 
y = h(t), Eq. (8) becomes, using Eq. (5), u = 


and the other at 
coy/h, and we see that this solution represents the case in which 
the wall at y = 0 is stationary and the wall at y = A is moving 
in its own plane with the constant velocity c: (postulating no slip 
at the walls 

This solution suggests the substitutions 
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Fic. 1. Variation of boundary-layer thickness with time on 


a suddenly moved flat plate with exponential density variation 


P/ Pr) v 
T J (p/p dt lt 
where p, is a reference density Eq. (7) becomes 
Ou/ Or = v-(07u/On- 1] 
where 3 “ip This, of course, is the well known one-dimen 


sional heat equation, numerous solutions of which are recorded in 
the literature. Any of these can be used for arbitrary variations 
of density with time, provided only that these variations make 
sense physically. The mathematical requirements on the function 
p(t) are simply that it must be integrable and that it must remain 
positive (or negative) so that the transformation given by Eq 
1()) exists and is one to one. Since in all cases 7» is proportional 
to y at any time /, the instantaneous velocity profile—that is, 
the variation of « with y—always has the characteristics of the 
profile in the constant density case, where 7» can be identified 
with y 

For example, a solution given in reference 1 is written in our 


notation 


u = mu erlfe ¢ m(1 — erf ¢ 12 
where ¢ = 9/2V wr 13 


which gives a boundary-layer like flow for small 7. Eq. (13 
may represent the solution for a suddenly moved infinite flat 


plate in a region where the density may be supposed to be varied 


in an arbitrary manner by the vertical motion of a plane at 
distance from the suddenly moved plate which is large compared 
to the thickness of the boundary layer 

\s an interesting case of this solution, assume that the density 
varies in an exponential manner 


pe” 14 


p 
where m is a constant. The boundary-layer thickness, 6, can be 
taken as the value of y for which ¢ = 2 But, from Eqs. (9), 


10), (13), and (14), we find 


and, therefore, 

6 = 4V v-/m Ve" -— 1 
This relation is plotted in Fig. 1. It is seen that the thickness 
actually decreases at large values of ¢. A comparison can be 
made with what might be called the solution for quasi-constant 


density —that is, when ¢ is given by 


t= y/2V wt 
with vy = pri pr/p) = vre 
This gives 5 = 4V v,/m ew? VV mt 
which is also plotted in Fig. 1 This too shows a decrease for 


large t, but the effect is considerably delayed and weaker than in 


the exact solution 
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| 
These results may be interpreted as follows. Initially 
boundary layer on the suddenly moved plate grows rapidly q 
to the high shear and low kinematic viscosity After a tim 
however, the increasing kinematic viscosity causes a reduction } 
the boundary-layer thickness ‘his thinning is larger in t 
exact case as a result of the convection of the outer lay: oft 
boundary layer toward the plate, as required by continuit | 
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a a 
On the Character of the Instability of the 
Laminar Boundary Layer Near the Nose of a 
Blunt Body* 
A. M. Kuethe 
Professor, Department of Aeronautical Engineering 
Michigan, Ann Arbor, Michigar 
January 24, 1958 
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sistance of K. R. Raman and Ralph Deitrick in obtaining the experiment 
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veloc erical nose 11.5 im. in diameter, at wind speeds 
luti ft. sec. in the University of Michigan 5 7 ft. low 
unnel. Sample records at various angles from tl 
g the same meridian plane are shown in Fig. 1 rw 
Office p re shown for each condition \pproximate relative 
finn : plification factors for u’ are indicated Phe disturbances 
The "\ e probably produced by surface roughness since the, 
riment id not appear at all meridian planes Phe bursts of approxi 


usoidal fluctuations indicate that the instability is 

ab f the Tollmien-Schlichting type 
Comparison of the reduced frequency 8,v/u,? where 8 7 
frequene f the disturbances) and u, is the velocity at the edg¢ 


f the boundary layer versus 4,6,» with the flat-plate (Sch 





re stability theory! is shown in Fig. 2. Some of 


the points refer to records made with a rough ‘button it the 


Phe fact that the disturbances appear and grow at Reynolds 





rs far below the critical value given by theory indicates 





the existence of source of instability that is neglected in the 
theory I suggest that this hitherto unrecognized source is the 
stretching of the vortex filaments in the region of diverging flow 
lear the stagnation point 

\pproximate analysis? of the destabilizing influences of stretcl 
ing and surface cooling compared with the stabilizing influencc 
f curvature indicates that near the nose the net effect is de 
stabilizing Phe analysis uses the Taylor-Gértler vortices 
which stretch in the favorable pressure gradient) as a model, 


ited out? that analogous effects exist for those of the 





Schlichting type Qualitative agreement between 
theory and the present experiment is obtained for the conditio: 
Ino heat transfer to the body 

It is pointed out further? that, while the effect of cooling 
plate is known to be stabilizing to Tollmien-Schlichting 


listurbances, when convex curvature is introduced the centrifugal 





loree acting on a displaced particle causes a destabilizing in 
fluence vell It seems reasonable, therefore, that if the 
urvature is large enough the effect of cooling a surface may be 
stab g to Tollmien-Schlichting disturbances 
\ theoretical and experimental program is under way 
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Relative Importance of Free-Stream Vorticity 
and Self-Imduced Pressure Gradient on a 
Flat-Plate Boundary Layer 


Herold Mirels 
Lewis Flight Propulsion Laboratory, NACA eveland, Ohio 
February 3, 1958 


_— PROBLEM of the incompressible laminar boundary layer on 
a flat plate in a free stream containing relatively small 
amount of vorticity, wo, has recently been discussed by Glauert 
and Li.* In reference 1, Glauert refutes an earlier claim by Li 


that the free-stream vorticity introduces a pressure gradient 


which should be considered in the solution for the effect of free 
stream vorticity on the boundary-laver development 
Glauert did not make an order-of-magnitude estimate of tl 

self-induced pressure gradient and, consequently, Li remained 


unconvineed.? The purpose of this note, in su 








is to point out that if a self-induced pressure gradient must be 
considered, it is the pressure gradient associated with the displace 
ment thickness of the zero-order (« 0) solution rhe effect of 
the free-stream vorticity on the self-induced ] ure gradient 1 
of higher order and is negligible Criteria for evaluating the 
relative importance of free-stream vorticity as compared witl 
the self-induced pressure-gr idient effect ire Is deve loped 
herein for both incompressible and supersonic flows Unless 
otherwise specified, the notation of references 1 and 2 is followed. 
Consider first the laminar boundary layer on a flat plate in 
incompressible uniform strean Various attempts ive beer 
made to improve the zero-order (Blasius) solution by expanding 
in inverse powers of the local Reynolds Numbers Ku has 
shown that a first-order correction may be obtained by including, 
in the usual boundary-laver equation, the pressure gradient in 
duced by the displacement thickness of the zer rder solution 


rhe form of the stream function ts thet 


y V uorX |foln t L/V Re VY L/x fF ; t ] 


where Re tol /v, L plate length, and F,(x, is of order 
he form of Fi(x, 7) is obtained by computing the pressure dis 
tribution due to the zero-order displacement thickness at the 
plate and in the wake If a small amount of vorticity 1s in the 


free stream it will have a higher-order effect on the self-induced 


pressure distribution and therefore does not influence Fi(x, 


Now, consider the incompressible free stream to have a small 
unt of vorticity rhe stream function can be expressed as 





vorticity in free stream / vorticity in 


The relative importance of the self-induced pressure gradient 
and small free-stream vorticity can now be evaluated. Consider 


ing x/L, fi(n), and F,(x, 7) to be of order 1, the ratio of the free 


stream vorticity effect to the self-induced pressure effect is 


If R is large, only the free-stream vorticity effect need | 
sidered. If FR is small, only the self-induced pressure gradient 
effect need be considered If R is of the order 1, both effects are 
equally important. However, in this case the resulting solution 


rr 


is simply the linear superposition of the pressure-gradient effect 


is computed by Kuo) and the free-stream vorticity effect (as 


computed by Glauert or Li for the case of zero induced pressure 


\ similar discussion can be made for the case of supersonic flow 
\ssuming a linear viscosity-temperature law, u/u C(7 
7), the zero-order displacement thickness can be expressed is 
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5o* = A WV mxC/u 
where A is a constant For Prandtl Number 1, and uniform wall 


temperature 7),, an approximate expression for A is® 


A = 1.73(Tw/To) + 0.332(7 — 1)M? 


Using linearized supersonic theory, the pressure distribution in 


duced by the zero-order displacement thickness is 


b — py)/p W2?/V Mo? — 1) (din*/dx 
Wo2/2V My? — 1) (AV C/V Re, VL/x (5 


The boundary-layer perturbation, due to the self-induced pressure 
gradient, is of the order of the coefficient of YW L/x in Eq. (5) 
The effect of 


boundary-layer development can be found by using a Howarth 


free-stream vorticity on compressible laminar 


transformation. The perturbation parameter can be expressed 


as awl V C/(u V Rez) (see reference 6 The ratio of the free 
stream vorticity effect to the self-induced pressure gradient 
effect is, roughly, 

R 2V My? — 1/yA Mo?) (wol/u 6 


One-Dimensional Transient Heat Flow in a 


Multilayer Slab 


W. F. Campbell 


Gas Dynamics Section, Division of Mechanical Engineering 
National Research Council, Ottawa, Canada 


January 17, 1958 


SUMMARY 


In a recent contribution tc the Readers’ Forum! Wassermann gave analytic 


solutions for the temperature in a double layer slab, with a triangular heat 
rate input at one face, insulated at the other, and with no thermal resistance 


at the interface His solutions were for the three particular cases 
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t 
For hypersonic Mach Numbers, 1 V/,? and | 
R = wole/toM 
Thus, at hypersonic speeds, the parameter wl /u must be 
tremely large in order for the free-stream vorticity effect ¢ 
comparable with the self-induced pressure effect 
| 
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V a2/a 24 
M4, the number of different and distinct paths from the plar ( 
of the heat source 0, to the plane irriving 
in the direction of which 
possible order, ¢ reflections back into layer N« it the ( 
interface b, j passages through the interface in t ( 
direction of decreasing, and reflections back 
layer No. 2 at the interface hus, for example, VM nstit 
) 1 if 0 and O if x0, M for 
M,(1, 1, 1 2. M2(1, 1, 1 4; and M (3, 2,0 
10, ete 
V (a the number of different and distinct paths from the 
of the heat source, x 0, to the plane arriving at 


which involve, in any { 


into layer No. | att T 


the direction of x decreasin 


sible order reflections back 
ki/ Vai ko/V/a2, ki/y ai<— k v/a2, | Va >) \ interface, x 6, 7 passages through the interface int 
direction of decreasing, and reflections back L 7 
I propose here to give the general solution to this problem, to indicate laver No. 2 at the interface Thus. for example ” 
briefly how it is obtained using the method of reference 2, and to point out 0. 0) 0 and N.(0. 0. 0 1 ene 
that the solutions given by Wassermann are incomplete for times longer Febr 
than the duration of the heat input 
; DISCUSSION ' 
SYMBOLS* I 
R 1)/Gu2 + 1) R | ‘HE GENERAL SOLUTION, assuming, of course, that the therm I 
S 2 “ 1) = wo properties are independent of temperature, is: art 
eta 
ensi 
\ bh 
SOnar Vai 63? [x + 26 + fb] Va + 27 +k) yr rob 
Ti = —" a (Riz)! (SiSu)!) (RakEX 4 Mii, j, R)Perke 
ak; i=07=0 k=0 = V aa20 
: : —x +2(¢+ 7)b] Vas +207 + Rk) (CL Va, | ee 
Ni(1, j, k)ierfe : lishee 
2 V aad ! 
= F(@) for 0<0< ad 
TT," = F(@) — [1/11 — a)| F(6 — af) for ab SOSA 
T,'"’ = F(@) — [1/11 — a@)| F(6 — ad) + [a/(1 — a)|F(O — 6) for 0>86 
a 80... 1, 63/2 , ; 
7 ee > z :. (Riz)' (Si 1 (So1)’ (Ra)” X 
adk; i=0 j7=0 k=0 
j : , 1+27¢+))b Vaot+ [x —b +27 +R) (L b)| Va 
M.(1, j, k)Perfe - u 
{ 2 a)020 
“<= 
; . 2i+ 7b Vat+ (b-—-x +2Aj,+k4+ 1) C1 Vial 
N2(i, J, k)erfe 
2 V aad 
= G(@) for 0<5< af 
T.” = G(6) — [1/11 — a)|G(6 — a6) for ab SO< 8A 
7. = G(6) — {1/11 — a)lIG(@ — aé) “a la (1 — a)|G(0 — 6) for 4>4@ 


* These are in addition to those of reference 1. 
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1eTm 


READERS’ 


d used to write down the preceding solution is based on 





the well-known analogy of one-dimensional heat flow with the 


f electric current in a simple transmission line having 

ly series resistance and parallel capacitance Thus when a 
emperature W ive,’ propagated from the heat source, arrives 
| where the thermal properties change abruptly, part of 
wav 1S reflected at, and part passes through the discon 
nuit The respective wave strengths are determined by 
sily calculated reflection and transmission coefficients Phe 


uperature at a point is then made up of all the contributions 





from the possible wave paths. Due to dissipation effects, in 
general, the longer the path length, the smaller the contribution 
For the problems where this method is particularly useful, the 


nsient temperature for ‘“‘short’’ times only is required and all 
ut a few terms are negligible Phis method is fully described 
reference 2 for many different types of heat source and for 
different boundary conditions. A simpler description of 
the method is also given in reference 3 for a somewhat more re 
stricted range of boundary conditions 
[he solutions given by Wassermann in reference | lack terms 


f the form F(6 #) and are thus incomplete for @ > 4 
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Asymmetric Starting for Hypersonic Wind 
Tunnels 


Robert H. Johnson 
Genera! Electric Research Laboratory, Schenectady, N.Y 
ebruary 3, 1958 


— PROBLEM of flow initiation limits the maximum size of 
models one may study in supersonic wind tunnels. This is 
particularly true of hypersonic tunnels because the viscous effects 
retarding the flow establishment are exaggerated by the low 
density and large velocity gradients characteristic of such flows 
\ blunt, axially symmetric model further aggravates the starting 
problem. Here a strong detached bow shock wave interacts with 
he boundary layers and the resulting feedback may cause chok 


ing. Consequently, full tunnel flow may never become estab 
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Fic. 1. Schematic diagram of helium-tunnel starting system 
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Fic. 2. 3/8-in. diameter cone-sphere in helium at 22 


The size and bluntness of a model is thus limited in a giver 
tunnel. Since the shock wave from a blunt model in a hypersonic 
tunnel is the predominant part of the entire diffuser flow, it is 
possible to increase slightly the maximum possible size of the 
model by adapting the diffuser design to suit the shock-wave 
pattern. Thus, when the model shape or size is changed, it is 
frequently necessary to redesign the diffuser 

\ new method for starting supersonic wind tunnels has been 
developed for use on the hypersonic helium tunnel at the General 
Electric Research Laboratory rhis tunnel is currently operating 
in the Mach 21 to 28 range 

] 


number of small air jets distributed along the length of the nozzle 


rhe starting method employs 
wall, as illustrated in Fig. 1 Individual quick-acting valves are 
provided for each air line so that air jets in any combination may 
be used \ maximum of five 0.010-in. diameter jet orifices have 
been used. The starting sequence proceeds in the following man 
ner 1) open air jet valves, (2) start tunnel flow, (3) close ai 
jet valves. Once flow has been established by use of the air jets, 
it remains stable after closing the jet valves rhe hysteresis is 
sufficient 

The asymmetric flow disturbance in the nozzle due to air jets 
changes the effective area ratio of the nozzle This results in a 
lower Mach Number on one side of the test section The shock 
wave, therefore, does not intersect the diffuser walls in a single 
transverse plane—i.e., the line of intersection between the shock 
and the boundaries form an irregular curve instead of a circle 
This ‘spreading out’’ of the shock wave causes a weaker inter 
action with the boundary layer and reduces the tendency of the 
model shock wave to converge, and the shock-induced boundary 
laver separation, with the resulting feedback, is retarded 

Similar improvement in the starting characteristics of the 
tunnel have been made by using a series of solid flow obstructions 
short pieces of stiff wire) in place of the air jets. Other gases, 


such as helium, also have been used 





Fic. 3. 1-in. diameter cone-sphere in helium at J = 22. 
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The effectiveness of the asymmetric starting procedure is 
shown by the increase in the maximum model size which can 
now be used in our 3.6-in. diameter Mach Number 22 hypersonic 
helium tunnel. Without this new starting procedure, a 3/8-in. 
diaineter model is the largest for which flow can be established 
The new technique, however, permits establishing flow with 
(Mach Number 22 flow 


about a 3/8-in. and a 1-in. blunt model is shown in Figs. 2 and 3 


models larger than 1 in. in diameter. 
Larger models, resulting in stronger shock waves, require 
Thus, the 
optimum number and location of air jets must be determined 


greater asymmetry during the tunnel starting cycle 


for any given model configuration. It is conceivable that 
models even larger than 1-in. diameter may be used in the present 


equipment. 


Note on ‘‘Stability Equations for Conical 
Shells”’ 


Paul Seide 

Guided Missile Research Division, The Ramo- Wooldridge 
Corporation, Los Angeles, Calif 

February 11, 1958 


i. A RECENT NOTE,! Professor Kempner shows that the writer's 
equations for buckling of conical shells? are 


(although derived differently ) 


identical with 


those given by Vlasov and cor- 


rects an error in the writer's derivation of a single eight-order 


differential equation for the displacement normal to the cone 


surface, w. He then derives approximate uncoupled equations 
for the meridional and circumferential displacements, “ and 7, 
in terms of z 

The contributions by Professor Kempner concerning the prior 
derivation of the stability equations and the correction to the 
eighth-order equation were noted in the published version of the 
writer’s paper? where it was also noted that a still earlier deri 
vation of identical stability equations was apparently given by 
Mushtari.! 


to make any approximations in deriving uncoupled equations 


Some further analysis reveals that it is unnecessary 


for the displacements uw and v in terms of the normal displacement 
wand that a further refinement in the eighth-order equation is 
possible 

It can be shown that the operator L, defined by the writer? 
(or H in Kempner’s note) can be written as 

(1/s*)V* s? (1 

so that the final form of the eighth-order equation for w is 
l O7u' 


V 4s DV *w Ns Na X 
‘ Os? 


( 1 Ow l s) re) ( Oz ) Ou 
T ne 2N SA ms T o 
s Os s? O82? ds \s 06 Os 


1 Ow Et cot?a 0? ( O7w 
(a) - 4 Ss zi iQ) (2 
5s Ob 55 Os? Os? 


with 
v2 = (07/ds?) + (1/s) (0/0s) + (1/s?) (07/06?) 
vi = T°? 
s = the distance from the cone vertex 
6 = the angular measurement around the cir 
cumference 
a = the cone semivertex angle 
7] = @ sin a 
Vso Na. Nso _ middle-surface force components and sur- 
S,90{ face force components prior to buckling 
E = Young’s modulus of cone material 
t = wall thickness 
v = Poisson's ratio 


D = Et 


[12(1 — vp? 
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The uncoupled equations for « and v can be derived folloy 
Two of the equations of equilibrium of an element of th ne 


given in terms of displacements as 


}s?V2 — | (1 + v)/2] (02/062)} w 
(0/06) }[(1 + v)/2] (0/ds 2} 
v(O/Os l v ( 
(0/06) }1(1 + v)/2] (0/ds)s +1 — vlu 
»{(1 v)/2 v l 
(1 + v)/2] (07/06?) jv = (0/06)w « 
It can be shown that the operators (s°V? — 1) and (0/0 
commutative; that is 
(s*v* 1) (0/ds)s(_ ) = (0/0s)s(s?¥? l 


Then the uncoupled equations can be readily obtained 


tsev= 1) $(07/062)|u = s°V4su 
([(d/ds)s — (1 + v)] (9262 -— 1 
hi(1 + v)/2] (0/ds)s 3 + v)} (02/06*)) weota (5 
(s?v2 — 1)? 1(07/08)|v = s'Visv = (0/06 


2 + v)s*(0?7/0s? 3S(0/ Os r 1 (0*/08*)|w cot 
which replace Eqs. (11) and (12) of Kempner’s paper t 
interesting to note that the substitution 


introduced by Mushtari and Sachenkov® transforms Eqs. (5a 
(5b) into differential equations with constant coefficients whi 


may be of some aid in solving cone stability problems 
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Application of the Witte Rearranging 
Method to a Typical Structural Matrix 


Bertram Klein 

Head, Structure Theory and Methods Group, Convair Astronautic 
General Dynamics Corporation, San Diego, Calif 

February 10, 1958 


yee AUTHOR has shown!~* how the analysis of a complicat 
structure may be reduced to the solution of a large matt 
equation. The matrix in the equation is large because the aut! 
has sacrificed size at the expense of simplicity of setup Ther 
has been doubt expressed as to the practicality of this appro 
because one has to invert a large, unsymmetric matrix. 1 
purpose of the present note is to illustrate how this predicame! 
may be overcome 

Recently, Bruno Witte has devised a certain scheme’ particu 
larly suited to the type of matrix obtained in the matric met! 
of structural analysis advocated by the author—a matrix wit 
large number of zero elements and representing a structure wit 
weak coupling. (Coupling may be described as the interdepend 


ence of forces and/or displacements The method consists 





epend- 


OnSists 


ept Tor the 


utomatically by a 


ngularization and, hence, 





READERS’ 


itiall systematic rearrangement of the orginal matrix (in- 
luding the right-hand-side column if desired) so as to approxi 
ite ly as possible an optimized triangular matrix. The 
urpose of the present note is to illustrate the final form of re- 
rang nt when the Witte method is applied to a typical matrix 
reviou et up by the author 
The trix shown in Fig. 1 is the one appearing at the bottom 
fp.4 f reference 1 The right-hand-side column is included 
venteenth column. The rearranged form is that shown 
Fig. 2. The numbers to the left and above in Fig. 2 represent 


| locations of the rows and columns in the original ma 
of Fig. 1 Notice that the matrix is now triangularized ex 
last Therefore the new form is close to a 


te triangularization 


two rows 
Also notice the large number of zero 
ments still appearing above the main diagonal in Fig 

Remember that the Witte rearranging method does not require 


rithmetic operations such as connected with elimination 


method can be programed to be carried out completely 


high-speed digital calculating machine 


time necessary to perform the rearrangement is small 


mpared to the additional time necessary to complete the tri- 


to solve the matric equation. It ap 


rs that matric equations encountered in actual structural 


I 


oblems in which the method of references 1-3 are used can be 


ranged by the Witte method in forms containing a large per 


e of triangularized rows and a large percentage of zero 


ats above the main diagonal. Consequently, the time re 
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quired to solve such matric equations and the storage necessary 
for the solution should be considerably less than the corresponding 
Therefore, 


quantities for solid matrices of the same sizes one 


may conclude that the objection raised against the method of 
references 1-3 is not justified in light of the fact that the size of 
a matrix of the type encountered in itself is not a criterion of the 


time and storage required to solve the problem 
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The Conservation Equations for 
Multicomponent Gas Mixtures in 
Arbitrary Coordinate Systems 


F. A. Williams* 

Daniel and Florence Guggenheim Jet Propulsion Center 
alitornia Institute of Technology, Pasadena, Calif 

February 14, 1958 


’ / SHE CONSERVATION equations for multicomponent reacting 
gas mixtures are generally given only in Cartesian or orthog- 


onal curvilinear coordinate systems. Actually, the conservation 


equations are easily expressed in an arbitrary coordinate system 


he general equations in 


We present t tensor notation and then 
indicate the simplifications which arise for orthogonal curvilinear 
coordinates 

a multicomponent gas mixture consisting of 


Greek 


pis the de nsity, 7 


Let us consider 
different 


}, is the mass fraction of species a, 


species which we identify by subscripts. If 


is the mass 
average velocity of the mixture, and Vj; is the diffusion velocity 


of species a, then the general equations for conservation of mass 


of component a, total mass, momentum, and energy may be 
written, respectively, ast 

0/dt) (pVa pla + (pVaVa x 

(Qp/dt) + (p 0 2 

p( Ov; /Or p = I op» VaFa 3 

a l 
0/Ot) (p pe = O I p  ¥ al a Va } 
Here wag is the mass rate of production of species @ per unit vol 


ume by chemical reactions, ?;; represents the total stress tensor, 


and | is the body force per unit mass acting on the component 
a. In Eq. (4 
of the mixture per unit mass, and Q, denotes the total energy flux 


the quantity e represents the total internal energy 


covariant com and 


f the ten 


vector. English subscripts refer to ponents, 


superscripts identify contravariant components « sors 

In Eqs. (1 is made of the 
tion for repeated English indices Phe 
respect to 


through (4) use summation conven 


semicolon represents co 


differentiation with the coordinates whose 


The 


Variant 


indices follow it covariant derivative of a tensor A,,*, for 


National Science Foundation Predoctoral Fellow 


T For a discussion of the reasoning by which Eqs. (1) through (4) can be 
derived, see, for example, P. G. Bergmann, / he Ti Re 
Prentice-Hall, New York, 1942 
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example, is given by can be accomplished by dividing each covariant component j of 


1.) (OA;4/ox | the tensor by #4; and multiplying each contravariant Component 
a3” 4ij eee a a = ‘ 
iby h The new velocities u; and l'g;, body forces fy;, energy 
where the Christoffel symbol of the second kind is defined by the flux vector g,, and stress tensor p,, are therefore defined by the 
relation equations 
bits (l/ojg £i1/ OX tT (O¢;1/OX (Og;;/OxX (6 
In Eq. (6) the quantity g,;; is the metric tensor which defines the 
coordinate system by means of the expression for the invariant ‘al i (9) 
infinitesimal distance, ] 
re =h hyp 
ds? = g;dx'dx? (7 
I , ; The summation convention does not apply here nor in the follow. 
In an orthogonal curvilinear coordinate system the metric . : . ; ’ y , ; a 
; , jee ing discussion because of the change in notation Phe distinction 
tensor is diagonal and we can define the quantities 4; = Os/Ox ae a tes : ag pee 
between covariance and contravariance disappears in this formu 
by »equation : : ‘ RES: 
7 erage lation and we shall use subscripts for all indices 
By substituting Eq. (8) into Eq. (6) we find that, in orthogonal 
curvilinear coordinates, the relation 
(h;)? E : 
bef = (Six /hx) (OWK/OX 
It then becomes convenient to redefine a tensor in such a way that h:b 2) (Oh; /d (10) 
all of its components are of the same physical dimensions. This ; : f < amas 
is valid. Here 6;; denotes the Kronecker delta. If Eqs. (9) and 
(10) are used in Eqs. (1) through (4), then the following simpl- 


fied forms for the conservation equations are obtained 
(0/Ot) (p Va) + > {a h;) (0/0x:) [pYa (ui + Uai)] [pVal(ui + Uai » (1/h;h;) (Oh;/Ox;)f = & (11) 
F rare 


(Op/odt) + >» [(1/h,;) (0/Ox;) (pu;) + pu 2 (1/h;h;) (Oh;/Ox (12) 


1 


p(ou;/Ot) + 2. pu;\(1/h;) (Ou; /Ox;) + (1/hjh;) [u(Oh;/Ox;) — uj(Oh; 
7 


> {(1/hj) (Opi;/Ox;) + pi; (1/hih;) [(Oh;/d% 


(0/odt)( pe) + z. [(1/h;) (0/Ox;) (peu: + gi) + (peu; + gi) ph (1/hjh 


L 1 


- > '(p h;) (Ou;/Ox + [(pji/h; f u;/h;) (Oh;/Ox;\{ 4 > Uailei) (14) 


i.?7 a t 


In Eqs. (13) and (14), use has been made of the symmetry of the pressure tensor p The derivation of Eqs. (11) through (14) for 
orthogonal curvilinear coordinates illustrates the method of application of the general conservation equations given in Eggs. (1) 
through (4) 


Third U.S. National Congress of Applied Mechanics 
Brown Unwersity, June 11-14, 1958 


The congress program, being distributed in May, will feature over 100 papers in Applied Mechanics, including mechanics 
of rigid bodies, mechanics of deformable solids, mechanics of fluids and gases, thermodynamics, and heat transfer. In ad- 
dition, there will be four general lectures by leading authorities. 

To reduce the cost of attending the symposium to the utmost, the following all-inclusive living arrangement is offered. 
This covers rooms in the new West Quadrangle of Brown University and all meals from Breakfast on Wednesday, June 11, 
through Luncheon on Saturday, June 14, but ot the Banquet on Thursday night ($5.00 The rates will be $24.00 for each 
of two persons sharing a room, or $30.00 for a person occupying a single room. Rooms will be available without additional 
charge for Tuesday night and Saturday night. 

Registration Fees: Advance Registration on or before June 7, $4.00. Registration on June 8 or later, $6.00. Please 
make checks payable to: Applied Mechanics Congress, not to any person 

Requests for copies of the program and letters concerning advance registration should be addressed to Miss E. M. Addi- 
son, Box F, Brown University, Providence 12, R.I. In registering, please indicate number of persons in party and family 
relationship, desire for single rooms or willingness to share rooms, estimated arrival and departure times, and need for 
parking space. 
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